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THE SYSTEM AND HARDWARE DESIGN OF 
REAL-TIME FAN BEAM SCATTEROMETER DATA PROCESSORS 


1.0 SUMMARY 

The Remote Sensing Center at Texas A&M University (TAMU) has 
completed a major system design effort for the National Aeronautics 
and Space Administration, Johnson Space Center (NASA/ JSC) for a class 
of real-time radar scatterometer signal processors. A system design 
is presented which employs state-of-the-art integrated circuits and 
signal processing technology. This report presents the details of 
the design of a system specifically tailored for the C- and L-band radar 
scatterometer systems. The design is sufficiently flexible, however, 
to accommodate any scatterometer system with only minor modifications 
in sampling frequencies and filter bandwidths. 

The design is based upon a 512 stage CCD (change coupled device) 
analog transversal filter which is employed to perform a 512 point 
discrete Fourier transform (DFT) via the chirp Z-transform (CZT) algo- 
rithm. The CZT implementation provides for rapid power spectral density 
computations on the radar return signal from which the normalized radar 
cross-section a can be computed for a number of desired incident angles 
at a real-time rate. 

The scatterometer processing system design is presented in suffi- 
cient detail to allow for construction of an engineering test mode of 
the processor directly from the design data contained in this report. 
Breadboard testing of key elements of the system have proven the validity 
and applicability of the design. 


2.0 INTRODUCTION 


2.1 Historical Background 

The NASA Johnson Space Center has acquired a series of airborne 
fan beam scatterometers since the mid 1960's for support of various earth 
and space related programs. Data from an early system proved helpful in 
identifying an ocean wind measurement technique. This work eventually led 
to the scatterometer system aboard SEASAT A. Raw scatterometer data from 
these early experiments were post processed into quantitative engineering 
units using a general purpose computer. The time and expense required to 
process data caused the delay between data product delivery and data 
acquisition to be excessive. As a result, utilization of the scatterometers 
were eventually removed from service. 

More recently, radar data requests by the NASA soil moisture program 
provided enough justification to warrant reinstating the 0.4 GHz, 1.6 GHz, 
and 13.3 GHz scatterometers into service, and constructing a new 4.75 GHz 
scatterometer. The earlier experiences with scatterometer data processing 
led to a program of developing faster, cheaper methods of data handling 
and processing. The initial thrust of the program produced a demonstration 
processor for use with the 13.3 GHz scatterometer. The philosophy behind 
this processor was to provide a real-time quick look capability for verify- 
ing data characteristics {i.e., system operation) and 2) provide a method 
of identifying those data to be post processed on a large computer to the 
accuracies needed for analysis. This processor was developed and construct- 
ed at Texas A&M University (TAMU) under NASA contract [1] to process a 
single polarization channel over a limited set of viewing angles. In 
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addition to this hardware processor, TAMU developed software routines for 
utilization on general purpose coinputers to reduce the raw data needed 
for analysis to calibrated engineering units. 

Recent technological advances have suggested that by combining analog 
and digital processing methods into a single processor, real-time processing 
and real-time rate post time processing of scatterometer data to calibrated 
engineering units could be accomplished. Such a system could provide all 
of the capability in terms of viewing angles, resolution and adaptability 
that the post time software systems previously developed could provide, 
with a potential for more accurate results as a result of eliminating the 
analog recorder when operating in a real-time mode. Such a system would 
provide experimenters with calibrated data on a timely basis with many 
fewer manhours from data flight to delivery. This realization provided 
the basis for the current efforts reported in this document. 

2.2 Design Objective and Overview 

Designs for two airborne radar scatterometer processors for use with 
the NASA 1.6 GHz and 4.75 GHz scatterometers were identified and analyzed. 

A portion of the processor was implemented to evaluate a "state-of-the-art" 
component proposed for use in the processor. This component permitted a 
standardized design approach which is extendable to other NASA fan beam 
scatterometers. The current effort exploited design experiences from pre- 
vious hardware and software processors to minimize significant error 
contributions and to assure repeatability in performance. However, 
innovations were also introduced as a result of the hybrid sampled analog 
and digital approach to provide a flexible operator/experimenter oriented 
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system. As a result of these new insights, major improvements were also 
identified for use in the purely software approaches to processing scatter- 
ometer data. This benefitted another NASA sponsored program to develop a 
more efficient software routine to process scatterometer data on an interim 
basis while the hardware processors undergo development. This latter 
effort ran concurrently with this program and afforded an opportunity 
to also test, anticipate and prove the characteristics of the hardware 
design described within this report. 

The hardware design features a chirp Z-transform (CZT) approach to 
filtering the Doppler spread radar return. The CZT is implemented with 
multiplying digital to analog converters and a charge coupled transversal 
filter. The filtering operation reduces to that of performing a discrete 
Fourier transform (DFT) of the radar return when represented in complex 
valued form. As a consequence, no Hilbert transform operation (sign seha- 
ing) is required to separate fore and aft returns. Both are provided 
simultaneously with considerable reduction in complexity. The subsequent 
processing, is actually limited to the aft data; however, the fore data is 
available within the processor should future efforts require it. 

There are many advantages in the CZT approach. It permits high fre- 
quency resolution of the Doppler return. As a consequence, the return 
may be measured with good angular resolution. This also permits the pro- 
cessor to adapt with changes in aircraft velocity to tract the desired 
viewing angles by simply using a different set of spectral outputs. It 
will also permit arbitrary choices in viewing angles if desired. The 
CZT approach can be readily applied to scatterometers operating at other 
wavelengths primarily by altering the sampling frequency. 


The power spectral density (PSD) of the total return is formed from 
the chirp Z transformed data. The formation of the PSD requires that the 
spectral data be detected (squared) and accumulated (averaged) over a 
period of time. To achieve the accuracy and the dynamic range required in 
scatterometry the detection and accumulation is accomplished digitally. 

The detected and averaged data is converted to estimates of the 
scattering coefficients o° at eight viewing angles over the aft sector. 

The conversion is implemented in software and requires the application of 
radar range, pattern data, viewing angle, and transmitted power to yield 
a calibrated result. In addition, the software permits interactive control 
of the processor. Since the computations and control are provided by 
software, any portion of the operating system can be altered should the 
need arise. 

The design approach was partially evaluated by actually implementing 
a subsystem of the scatterometer processor. An evaluation to this detail 
was required to validate the performance and dynamic range of the charge 
coupled devices and associated circuitry since this is a "state-of-the-art" 
item. 

This report describes the system design theory, the system operating 
rationale and architecture, the hardware and software designs and an 
evaluation of the CZT approach for the scatterometer processors profiled 
above. In particular. Section 3.0 develops the system design theory back- 
ground. The characteristics of CW fan beam scatterometers are related to 
the scatterometer equation to identify the measurement theory. It is 
shown that the angular scattering characteristic can be resolved by 
estimating the PSD of the radar return. The precision by which the PSD 


estimated is related to the bandwidth - time product by analogy with clas- 
sical fading theory. The technique by which the fore and aft spectra are 
separated using a DFT method is then identified. The DFT is related to 
the CZT and the means by which the CZT may be implemented is then established. 

Section 4.0 is dedicated to establishing a suitable operating ra- 
tionale for the processor. Trade-offs between angular resolution, ground 
resolution, precision, and beam resolution are established and evaluated 
to identify a suitable operating mode to satisfy user requirements and 
system constraints. 

Section 5.0 describes in detail the system architecture and carefully 
distinguishes between the target system and the engineering development 
model proposed within this effort. An overview of the internal operation 
of the system is also presented. 

Section 6.0 discloses the hardware design rationale for the PSD 
estimation subsystem. The results of evaluating a breadboard model of the 
PSD subsystem is also included. 

The micro-processing subsystem and its interface units are treated in 
Section 7.0. A high speed floating point processor manufactured by 
Advanced Micro-Computer Corporation (AMC) is identified as a potential 
replacement for the 80/20 identified in work prior to this contract. The 
AMC mono-board computer is fully compatible with the 80/20 system but is 
much faster and permits efficient programming. The design and operation 
of three interface units which allow communication with the micro-processor 
are also described. 

Section 8.0 is devoted to the system controller and processing soft- 
ware. There the software design rationale, architecture and operating 
procedure are described. 
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The conclusions and recommendations appear in Section 9.0. The entire 
document is supported by a number of appendices. In addition to treating 
special topics, the appendices serve to preserve parts listings for each 
subassembly and vendors literature on special components. Schematics are 
provided in blue print form in a special container. Special 8080 software 
routines required to run and evaluate the CZT demonstration processor are 
also included in the appendices. 


3.0 DESIGN THEORY 


3.1 Introduction 

Airborne fan beam scatterometers permit simultaneous backscatter 
observations over a range of incident angles. By confining the antenna 
beam width in the crosstrack dimension and spreading the beam in the 
along track dimension, Doppler filtering may be employed in a CW system 
to resolve the average return power at various incident angles, each of 
which is spanned by a small angular window as illustrated by Figure 3.1 
Combinations of transmit and receive polarizations permit like and 
cross polarized scattering properties of distributed targets to be mea- 
sured. When the aircraft is flown over the same distributed target at 
different headings about the compass, the azimuthal as well as the inci 
dent angular behaviors may be documented. 


3.2 The Scatterometer Equation and Fan Beam Systems 

For a large class of distributed targets the returns from slightly 
different angular directions are essentially uncorrelated. Where a par- 
ticular direction is denoted by (G,^) within the coordinate system of 
Figure 3.2, the total return power may be described by summing returns 
from patches of the target located in various angular directions (6^,(()j). 
If the radar cross section in direction is denoted by a (6., <{)•)» 

1 J PH I J 

the total return power may be expressed as 


Wr* 


(4Tr) 


W. 


N 

E 

i=l 


M 


j=l 


(3.1) 
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where X = radar wavelength 

6tp “ transmit pattern directivity for polarization state p 
®rq " i^eceive pattern directivity for polarization state q 
= total transmitted power 
= h/cos0^. 

h = aircraft altitude 

p,q indices denoting the transmit and receive polarization 
states, respectively 

In order to primarily discriminate the backscatter within an angular 
sector of the incident angle 0, a Doppler filter, whose normalized 
transfer function is given by H(o)), is employed (see Figure 3.1). The 
portion of the return power appearing at the output of the filter is 
therefore given by 

XW. N M 


(3.2) 


where 


u). . = 4 tt V sin0. cos0 ./X 

I J * J 


(3.3) 


is the radian frequency associated with the patch in direction (0.,tf).). 
When the normalized scattering coefficient (0,(|)) is introduced, the 

rM 

double summation may be replace with a double integral given by 


X^W, 


w. 


Gj. g 

ff -Ajca 


^pq|H(^) 


dA 


' (471)=* R** 

In the interpretation of W it is important to realize that |h(u))|" 


(3.4) 
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participates within the integration since w is dependent on In an 

ideal fan beam scatterometer, (p ranges over a small interval about zero 
since the crosstrack antenna beam v/idth is small. As a consequence 

(jo K 4tt V sin 0/X (3.5) 

It is then observed that |H(w)|^ plays the role of a normalized antenna 
pattern having discriminatory power in the 0 dimension. 

When the bandwidth of the Doppler filter is sufficiently narrow, 
a° may be regarded as constant over the area A spanned by the Doppler 

r H 

bandwidth and the crosstrack beamwidth. The scatterometer equation (3.4) 
then reduces to the form 




H ^cos'edA 

t pq^ o' tp rq 


(47t)' 


m: 

h** 


(3.6) 


where 0^ is the incident angle corresponding to the center frequency 
of the Doppler filter, i.e., 

Wq = 4tt V sin 0 q/A (3.7) 

The recovery of a(0^) is, therefore, dependent on measurement of W^, 
and h and on an estimate of the double integral. In this regard the 
integral is often approximated by introducing an effective area Ag^^ 
so that 

ff G^G^jHl^cos** 0 dA = G^(0g,O)G^(0g.O) cos*^ 0^ A^^^ (3.8) 

The factor Ag^^ is often formulated on the basis of computed experience. 


3.3 Precision in Estimating 0 °^ 

In order to estimate a° at a set of incident angles 0 . , i=>l,2,..., 

PH ' 

n, a bank of filters is required. The center frequency of each filter 
is chosen in accord with equation (3.7). If S(w) is the spectral density 
of the return signal, then the output of the ith filter is given by 

W^( 0 oi) = 2 /Q|Hi(o3)|2 s(w)df (3.9) 


An alternative method could, instead, measure (estimate) the power spectral 
density (PSD) of the return signal and then form the return power through 
an integration, viz.. 


w = 2 ; 


ui 

^li 


S(w)df 


(3.10) 


where f^^^. and f^^. are the upper and lower corner frequencies associated 
with the ith angle. In the latter case the measurement of a° has been 

PH 

reduced to a problem in estimating the power spectral density. 

Regardless of the approach, it is necessary to reduce the variance 
in the estimate of the mean power return to assure a good estimate of 
the average scattering coefficient a° . As is well known, the radar return 

PH 

is characterized by heavy fading since the signal has a Rayleigh-like 
amplitude distribution [2]. As a consequence, it is difficult to estimate 
the mean squared statistic of such a signal. The theory for the precision 
in the estimating the mean squared statistic appears in references [2], [31 
and others. The extrapolation of this theory to the case where the PSD is 
to be estimated can be established on an intutive basis and is made pre- 
cise in reference [4]. The PSD is estimated from the periodogram and is 
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defined by 


Sm(I<) = M 


N-1 

E 

n=0 


s(n) e 


-j 2 tt kn/N 


(3.11) 


for the kth spectral line foir a signal represented in sampled form 
(s(o), s(l),... ,s(N-1)}. 

In the case where analog filtering, detection and integration is 
performed, it is well known that the standard deviation o in the estimate 
of is given by [3] 


a = W^//TT (3.12) 

where B is the pre-detection (effective) bandwidth and T is the integration 
period. The dependence of the variance ratio on the BT product is 

illustrated in Figure 3.3. 

Although the variance reduction is a good indication of the improve- 
ment in the estimate of W^, the precision of a system is better conveyed 
by a statistical 90% confidence interval. The 90% confidence interval for 
a BT product of 10 or better can be approximated by [5] 


- 1.645a < + 1.645a (3.13) 

where is the estimate of W^,. When the span of this confidence interval 
is expressed in dB, it can be written as 


p = m | ng 1 + 1.645/ /TT (3.14) 

1 - 1.645/ /IT 

The dependence of this precision factor on the time bandwidth product is il 
luGtrated in Figure 3.4. It is observed that a + 1 dB confidence interval 









requires a BT product of 50. Furthermore, it is noted that the improve- 
ment in precision becomes less rapid when BT >100. 

The above theory can also be applied to the case where the PSD is 
determined from the periodogram. To obtain precision in the estimate 
of the PSD, the periodogram must be averaged. This can be accomplished 
by averaging spectral estimates from a sequence of records and by smooth- 
ing adjacp’^t spectral lines from a single record. To make the association 
between the precision for an analog processor with that for a PSD proces- 
sor, it should be noted that the PSD estimates are based on the discrete 
Fourier transform (DFT) of the return signal s(t) (see equation (3.18)). 

The effective bandwidth associated with DFT processing (filtering) is given 
by 

Af = Y (3.15) 

where T is the duration of the signal record. It has been assumed that 
the record is unweighted. Therefore, the time bandwidth product associat- 
ed with a single line from a single record is simply 

BT = AfT (3.16) 

or 1. Since the spectral estimates from a single record are poorly 
correlated, the BT product can be enhanced by averaging over a window of 
adjacent spectral lines. The BT product can be further enhanced by 
averaging line estimates from a sequence of non-overlapping records. 

So, if adjacent lines from each record are smoothed and non-over- 
lapping records are processed, the resulting BT product for a filter of 
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(3.17) 


width NpAf is given by 

BT = (NpAf)(N^T) 

or simply NpN^^. This result indicates that the precision improvement is 
identical for analog and PSD processing. 


3.4 Power Spectral Estimation Using the Chirp Z-Transform 

As indicated above, the estimate of the PSD can be based on the DFT 
of a radar return record. If s(n), n=0,l ,2,. . . ,N-1 is an N point sequence 
representing the return signal s(t) over a time interval T, then the DFT 
of s(t) is defined as 

F(k) = V s(n) (3.18) 

n-o 

where k e{0,l ,2,. . . ,N-1}. F(k) is interpreted as the spectral amplitude of 
s(t) at a frequency of k/T when 0 ^ k ^ and at a frequency of -(N-k)/ T 
for ^ < k < N - 1. It has been assumed that N is even. The DFT formula- 
tion may be modified through the use of the identity 


2nk = n^ + k^ - (k - n)^ 


(3.19) 


to permit an implementation of the DFT by hardware. The identity results 
in a DFT given by 


F(k) = e 


^ «-JTrk^/N s(n)e’^^'^^/'^ e 


(3.20) 


n=o 


The implication of the above result is that s(n^) must be first down-chirped 
th convolved with an up-chirp and then post multi- 


wi 


plied by a down-chirp e The pre and post multiplications may be 


performed by analog methods and the convolution may be formed with trans- 
versal filters using charge coupled devices. The transversal filter 
requires 2N-1 stages to implement the DFT. When forming the PSD, the 
post chirp may be discarded since it does not affect the amplitude. 

A more efficient means for implementing the transversal filter, 
requiring only N stages in the transversal filter, is based on the slid- 
ing DFT. The sliding transform is defined as 

k+N-1 

F.(k) = Z s(n)e"^'^^"'^^'^ (3.21) 

^ n=k 

and differs from the non-sliding version in that the input sequence is 
shifted forward one sample for each new spectral estimate. The transform 
consequently operates continuously. Spectral estimates on the same line 
are updated every N samples since is modulo N in the parameter 

k. As a result of the sliding action, phase information is destroyed; 
nevertheless, the magnitude information important to the PSD estimation 
is preserved. 

It can be shown that through the use of the identity of equation 
(3.19), the sliding transform can be rewritten as 

F (k) = z s(k - m + Nje’J”"' ' ^ ^ 

^ m=l (3.22) 

This result implies that the input must be pre-chirped by a factor 
g-jTr(m-N) /N^ convolved with and then post multiplied with 

g-jTTk /N sliding DFT. When the PSD is required, the post 
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multiplication may be replaced with a squaring operation to yield 
jF^(k)l^. Since the pre-multiplication is periodic in m, the CZT filter 
can operate on the input signal continuously with the transversal filter 
only requiring N stages. 

3.5 The Application of the CZT to Doppler Filtering 

3.5.1 Signal Processing Theory 

A simplified block diagram of the CW fan beam scatterometer is 
illustrated in Figure 3.5. The transmitter illuminates the terrain at a 
radian frequency of The backscattered signal arriving at the receiver 
may be denoted as 

s(t) = a(t) cos [Wpt + c|)(t)] (3.23) 

where a(t) and i|)(t) may be regarded as random variables. The spectrum 
of s(t) is depicted symbolically in Figure 3.6a. The fore and aft spectra 
are distinguished from one another by "coloring" the fore spectrum as a 
rectangle and the aft spectrum as a triangle. As indicated in the receiv- 
er chain, the return signal is split equally into two channels. This sig- 
nal in the upper channel is coherently demodulated with cos w^t and low 
pass filtered to yield 

x(t) = Js a(t) cos(j)(t) (3.24) 

The upper channel is commonly called the cosine channel or the in-phase 
(I) channel. Demodulation and low pass filtering in the lower channel 
yields 


y(t) = a(t) sint|)(t) 


(3.25) 
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FIGURE 3.5 A SIMPLIFIED BLOCK DIAGRAM OF THE CW FAN BEAM 
SCATTEROMETER 







FIGURE 3.6 SPECTRA OF INTEREST 



where product modulation with the quadrature reference sin co^t has 
occurred. The lower channel is commonly called the sine channel or the 
quadrature (Q) channel. The spectra of x(t) and y(t) are illustrated 
in Figures 3.6b and 3.6c, respectively. 

A comparison of Figures 3.6b and 3.6c shows that the aft and fore 
spectra can be retrieved simultaneously if a complex signal 

z(t) = x(t) + j y(t) (3.26) 

is formed. An examination of the spectrum of z(t) indicates that the aft 
spectrum occurs for w > 0 and the fore spectrum for u < 0. 

Since the DFT can also be applied to complex signals as well as real 
signals, the above result shows that the fore and aft spectra can be 
simultaneously filtered using CZT techniques to implement the DFT. The 
configuration for implementing the CZT using charge coupled devices and 
analog multipliers is described below. 

3.5.2 The Implementation Technique 

It is advantageous to use the sliding version of the CZT, since 

many sequential measurements are required to estimate the PSD. The 

actual implementation method is best understood by separating the real 

and imaginary parts of the sliding CZT. When the sliding CZT of z(t) 

is taken, Zg(k) can be rewritten in the form 
N 

|Zg(k)| = 2 [cosTr(m-N)^/N + j sim (m-N)^/Nl 

m=l 

[x(k-m+N) + jy(k-m+N)]. 


(3.27) 


A careful interpretation of the arguments within the above magnitude 
suggests the implementation technique shown in Figure 3.7. Both x(t) 
and y(t) are pre-chi rped and appropriately summed to form the real and 
imaginary valued entries into the transversal filter bank. Four trans- 
versal filters are required to form the cross multiplication products of 
the complex valued signal entering the filters. The real and imaginary 
parts of the CZT without the post-multiplication are formed by differenc- 
ing the outputs of the transversal filters as illustrated in Figure 3.7. 



FIGURE 3.7 A TECHNIQUE FOR IMPLEMENTING THE CZT WITHOUT POST MULTIPLICATION 







4.0 SYSTEM DESIGN RATIONALE 


4. 1 Introduction 

As is the case with any system, the actual design is a compromise 
between user requirements and system constraints. The design of real 
time processors for the NASA scatterometers is no exception. Ideally 
the scatterometer should provide high precision estimates of a° with 
infinitesimal angular and ground resolutions. However, as will be shown 
below, the precision, angular resolution and ground resolution interact 
in such a way to prevent maximizing all three parameters simultaneously. 

In addition to maximizing the precision and resolution parameters, the 
user is also interested in achieving a reasonable accuracy to permit 
comparative analysis of the processed data at different view angles and 
polarizations and with scatterometer data from other sources (also 
presumably calibrated). 

A list of the factors which can potentially influence the performance 
of the scatterometer/processor system is shown in Table 1. The source 
of error is described in the left-hand column. The system performance 
factors most influenced by the error source is reflected in the middle 
column. The origin of the error within the system is identified in the 
right-hand columns. Most of the performance parameters are manageable 
by the processor provided that the scatterometer has been appropriately 
designed. Those that are manageable are treated below as well as in 
subsequent chapters to develop the processor design rationale. 
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TABLE 4.1 FACTORS INFLUENCING SYSTEM PERFORMANCE 


Error Source 

Affected 

Performance 

Factor 

Origin In Sj 

i 

astern 

U 

OJ 

i 

o 

u. 

a» 

<0 

o 

CO 

u 

o 

1/) 

0) 

u 

o 

L. 

Q. 

4^ 

fO 

£ 

•r- 

< 

4-> 

S, 

? 

1. Fading Signal 

Precision 




n 

2. Finite Doppler Bandwidth 

Precision & Angular 
Resolution 


X 


m 

3. Finite Record Length 

Ground Resol tul on 
& Precision 


X 


■ 

4. Filter Sidelobe Level 

Accuracy 


X 


in 

5. Bit Truncation 

Precision 


X 



6. Inversion Approximation 

Accuracy 


X 



7. Uncertainty in Altitude 

Accuracy 



X 

ni 

8. Aircraft Attitude 

a) Illuminated area 

b) Polarization 
decomposition 

Accuracy 



X 

X 

1 

9. Transmitted Power 

Accuracy 

X 



Hi 

10. Polarization 

Accuracy 

X 



Hi 

11. Non-Statlonary Return 

Accuracy 




B 

12. Beamwldth 

Accuracy & Angular 
Resolution 

X 



■ 

13. Pattern Sidelobes 

Accuracy 

X 



H|i 

14. Pattern Gain 

Accuracy 

X 
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4.2 Definition of System Design Parameters 

To identify the processor's mode of operation it is important to 
define various parameters associated with fan beam systems. In this 
regard such terms as angular resolution, ground resolution, scan length, 
beam resolution, ground track coverage, etc. must be clarified to arrive 
at the impact of these parameters on the system design. 


Angular Resolution 

As indicated in Section 3.2 the angular resolution of a fan beam 
scatterometer is dictated by the physical beamwidth in the crosstrack 
dimension and by the bandwidth of the Doppler filter in the intrack di- 
mension. If H(w) denotes the normalized voltage transfer function of 
the Doppler filter, then an affective bandwidth may be defined as 

00 

B=i; |HU)l'a) (4.1) 

—00 

where the normalization has been applied so that {| H(w) j) = 1 
when H(o)) represented in low pass form. An effective intrack beamwidth 
may be related to the effective bandwidth through the Doppler relationship 


A0 


when t|) 0. 


= AB 
2VCOS0 

A0 is defined to be the angular resolution. 


(4.2) 


Beam Resolution 

The beam resolution in the cross track dimension is given by 


= (h tan0) A(p 


(4.3) 


where A(p is the crosstrack angle subtended by the two way beam when pro- 
jected on the ground plane. The beam resolution in the intrack dimen- 
sion is 

”■ coiFe 

See Figure 4.1 to clarify these definitions. 

Scan Length 

Scan length is simply the ground track distance vT traversed by the 
aircraft during a single integration period T where v is the ground 
velocity. See Figure 4.2. 

Ground Track Coverage 

The ground track coverage is that entire length over which a 
radar return was observed. The coverage includes the scan length as 
well as the initial coverage within the beam, ie., 

I-C “ Pb 

See Figure 4.2. 

Ground Resolution 

The ground resolution may be defined in several ways. However, 
for the purposes of this design effort the ground resolution Is defined 
as an effective ground length over which radar returns have primarily 
contributed to the measurement as implied in Figure 4.2. It consequent- 
ly emphasizes that portion of the ground track which is repeatedly in 
view within the beam subtended by AG. 



1 


^0 


vT = scan length 

A0 = angular resolution 

p = effective ground 
resolution 

pg = beam resolution 



Lj, = length of coverage 



• V 



FIGURE 4.3 COMPARISON OF THE RESOLUTION PARAMETERS WHEN 
(a) vT^ < pg, (b) vT^ = pg and (c) vT > pg 
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Three cases may be identified as illustrated in Figure 4.3. The 
ground resolution is therefore defined as 


^ vT + p 

— ^ Pb 


p = < Pg if vT = pg 
vT if vT > Pg 

From the above definitions an important observation can be with- 
drawn. When the three ground resolution cases are ordered as implied 
In Figure 4.3 it may be shown that 


(4.5) 


and 



(4.6) 


L. = £2. 

TT P3 

where p^ is the ground resolution corresponding to T^. For a constant 
Doppler bandwidth these results imply the following conclusion: 

Case 2 maximizes the ground resolution consistent with maximizing 
the precision. It is also interesting to note that Case 2 maximizes the 
precision for a given coverage interval as demonstrated in Appendix A. 
In this case the design criterion requires that 

Pg = vT (4.8) 

at any incident angle. This conclusion will be used in identifying 
appropriate design theories below. 


4.3 The Theory for Constant Precision, Constant Angular Resolution 
and Constant Ground Resolution Designs 

Among the many system designs which could be considered it is 
helpful to limit considerations to three basic design approaches: 

(1) constant precision, (2) constant angular resolution and (3) constant 
ground resolution. The theory for each is presented below and their 
characteristics are compared in a final subsection. 

4.3.1 A Constant Precision Design 

A constant precision design approach requires a constant BT product 
at each incident angle to be processed. The implication deduced from 
Section 4.2 is helpful in assigning B and T so as to achieve an accept- 
able BT product. At the smallest incident angle 9^ the BT product may 
be maximized for a given coverage L. by requiring 


Pb' ' 

Then 


and 


T - h A9i 
" V cos^Gi 


(4.9) 


(4.10) 


B = 2v cosBi A9i/X 

The viewing window at 9j is therefore specified as 


AG 




3T \ cos^Gi 


2h 


(4.11) 


(4.12) 


where BT is chosen to achieve the desired precision. When A0^ is with- 
drawn from equation (4.12), T and B are uniquely assigned by equations 


(4.10) and (4.11), respectively. On pragmatic grounds the integration 
time must be constant at all viewing angles. Therefore B is also constant. 
As a result, the angular resolution, beam resolution, and ground resolution 
at the remaining view angles 9j^ become 


A0 


XB 


k " 2vcos0, 


(4.13) 


„ - XhB 

^Bk " ^vcos^0, 


(4.14) 


and 


n = XhB f-, . cos ^01 ■, 

^k “ 4vcos*0i ^ "coPiT” ^ 


(4.15) 


From the above results it is noted that the angular resolution is inverse- 
ly proportional to v cos0|^, the beam resolution is proportional to h/vcos^0|^ 
and the ground resolution is proportional to h {1 + 

COS w 

4.3.2 Constant Angular Resolution 

A constant angular resolution approach requires that A0j^ be constant 
at all viewing angles, say A0. Once A6 is assigned, the bandwidth at each 
0|^ is given by 




2v cos0|^A0 


(4.16) 


and the beam resolution by 


= . hA0 _ 
^Bk ” cos'^0, 


(4.17) 


The precision and ground resolution require a rationale to assign T. 
Once again it is convenient to maximize BT at the smallest incident 
angle for a given coverage. This requires that 


The remaining parameters then become 


(4.18) 


T = 


hA0 


V COS^Oi 


(4.19) 


2hcos0|^ A^0 

= xcos^e. 


(4.20) 


and 


Pk 


r , . COS^01 1 

2cos=^0i 


hA0 


(4.21) 


at each viewing angle 0|^ > 0^. 

4.3.3 Constant Ground Resolution 

A constant ground resolution approach requires that p be constant 
at each incideht angle. With p specified, the precision may be maximized 
at each incident angle by requiring p = pg = vT in accord with Section 4.2. 
As a result of this imposition 

A0|^ = pcos^0|^/h (4.22) 

T = p/v (4.23) 

and 
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(4.24) 


(BT),^ = 2p2 cos®e,^/A h 


4.4 A Comparison of the Design Approaches 

To evaluate the three design approaches, the nominal design guide- 
lines shown in Table 4.2 were employed for C band and L band systems. 

The guidelines were primarily applied at the smallest incident 
angle and were allowed to vary at the larger incident angles depending 
on the design approach. The results of this evolution are shown graph- 
ically in Figures 4.4 and 4.11. 

The first four graphs apply to the C band processor whereas the 
latter four apply to the L band processor., Each figure identifies and 
compares a single system performance parameter over the entire range of 
incident angles for the three design approaches. The graphs are para- 
metrically identified by the design approach: CP = constant precision, 
CGR = constant ground resolution and CAR = constant angular resolution. 

Table 4.2 Nominal Design Guidelines 


Parameter 

Value 

Units 

Aircraft velocity 

150 

knots 

Aircraft altitude 

1500 

feet 

Angular resolution (nominal) 

C band 

3 

degrees 

L band 

6 

degrees 

Ground resolution(nominal) 



C band 

25 

meters 

L band 

50 

meters 

Precision factor (nominal) 

C band 

50 


L band 

50 




FIGURE 4.4 THE DEPENDENCE OF PROCESSING PRECISION ON 
VIEW ANGLE FOR THE THREE DESIGN APPROACHES 



FIGURE 4.5 THE DEPENDENCE OF ANGULAR RESOLUTION ON 

VIEWING ANGLE FOR THE THREE DESIGN APPROACHES 







FIGURE 4.7 THE DEPENDENCE OF BEAM RESOLUTION ON VIEWING 
ANGLE FOR THE THREE DESIGN APPROACHES 






FIGURE 4.10 THE DEPENDENCE OF THE GROUND RESOLUTION ON 

VIEWING ANGLE FOR THE THREE DESIGN APPROACHES 



VIEWING ANGLE (degrees) 

FIGURE 4.11 THE DEPENDENCE OF BEAM RESOLUTION ON VIEWING 
ANGLE FOR THE THREE DESIGN APPROACHES 



From these graphs the following features are noted: 

1) . When constant precision is imposed, the angular, beam and 
ground resolutions degrade with incident angle; however, the degrada- 
tions, are only significant for angles greater than 45°. 

2) . When constant ground resolution is imposed, constant beam 
resolution is also realized. The precision, however, degrades at the 
larger incident angles but is useable to 50°. The £• , nar resolution 
increases rapidly at the large incident angles. The very high angular 
resolution at the large incident angles may preclude selecting a correct 
pattern gain at the large incident angles when converting to 0 °. 

3) . When constant angular resolution is imposed, the beam resolu- 
tion and ground resolution degrade at large incident angles. The per- 
formance, in general, lies between the constant ground resolution and 
constant precision design approaches. 

If a single design approach were to be selected among the three, 
it is apparent that the constant angular resolution approach represents 
a good compromise between constant precision and constant ground resol- 
ution. The constant precision and constant ground resolution designs 
may, however, suit some experiments better. The constant precision design 
is attractive in those cases where high precision is required on a single 
cell, particularly at the larger viewing angles. The constant precision 
approach may also be helpful in those cases where a large spatial average 
is required. The constant ground resolution approach is attractive for 
those applications where well metered - high resolution data are required 
along the intrack dimension. This design approach is consequently attrac- 
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tive for those targets which are highly nonhomogeneous . 

Since the scatterometer processor is under software control, it 
is conceivable to provide an experimenter with any design option. The 
constant angular resolution design represents a good compromise among 
the approaches, however, when a single approach must be taken. 




\ 



5.0 SYSTEM ARCHITECTURE AND OVERVIEW 


5 . 1 Target and Development System Architectures 

The objective of the scatterometer processing system is to 
provide real and post time conversion of two channels of scatterometer 
data, like and cross polarized signals, into a° estimates at eight (8) 
viewing angles: 5°, 10°, 15°, 20°, 30°, 40°, 50° and 60°. Processor designs 

were to be developed for the NASA 4.75 GHz and 1.6 GHz fan beam 
scatterometers. The efforts were 1) to emphasize a standardize design 
approach suitable for use with these scatterometers as well as future 
scatterometer systems and 2) to utilize design experience from the 
previous 13.3 GHz scatterometer processor project. An appropriate target 
system architecture using the CZT approach to Doppler processing and 
meeting the above stated objectives is illustrated in Figure 5.1. The 
system consists of two major subsystems, viz., the PSD estimation subsystem 
and the micro-processing subsystem. A number of interface units are 
also provided to permit control of the system, entry of data and storage 
of processed data. Also a special alignment generator not required 
during operation but helpful in aligning the CZT filtering unit prior 
operation is shown. 

The PSD estimation subsystem has two channels to handle like and 
cross polarized return signals. Each channel converts the quadrature 
signals into discrete PSD estimates over the fore and aft Doppler spectra 
simultaneously. The estimation technique is based on the CZT technique 
which may be regarded as an analog technique of implementing the discrete 
Fourier transform (DFT). The detection (squaring) and accumulation 
(averaging) of the spectral amplitudes is performed digitally by a 
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INTERFACE 



FIGURE 5.1 ARCHITECTURE OF THE TARGET SYSTEM 



















dedicated processor. 

The micro-processing subsystem is composed of a digital micro- 
computer, memory and associated software. The role of the micro-pro- 
cessing subsystem is 1) to control the scatterometer processor through 
a software operating system, 2) to select appropriate PSD estimates and 
to use them within the radar scatterometer equation to invert for 
values on the eight viewing angles and 3) to accept, store end transfer 
various data required by the processor or the experimenter. 

Communication between the micro-processor and various peripheral 
systems is provided by the interface units indicated in Figure 5.1. 

Among these interface units is the operator interface and display. All 
system functions are initiated through this interface by the operator. 

The display is also used by the micro-processor. 

The system of Figure 5.1 represents but a single architecture for 
the target system. It may not be the architecture of the final system 
but it will be close. Implied in Figure 5.1 is a single micro-processor 
to service both scatterometer channels. However, in view of the amounts 
of data and the number of computations involved, a single micro-processor 
may be unable to convert both channels of scatterometer to a° values at 
eight angles without large gaps between successive ground cells. A 
fast floating point processor such as the Advanced Micro Computer 
95/4000 may permit use of a single processor. However, it may be 
necessary to dedicate a micro-processor to each scatterometer channel. 

To make an assessment of the number of processors requires that the 
software design be reduced to machine coding and that the machine be 
specified. 
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In view of this uncertainty and with full realization that a complete 
design effort requires an iterative effort between paper design (called 
for in this contract) and laboratory evaluation, an engineering model 
which is expandable to the target system was actually designed. The 
architecture of the proposed design is illustrated by the block diagram 
of Figure 5.2. The development model is restricted to single PSD 
estimation channel and a single micro-processor. The software is suffi- 
ciently general to permit processing of like or cross polarized data for 
any transmit polarization. However, the processor must be cued 
externally as to which polarization channel the PSD estimation subsystem 
is connected. A single channel processor of this type is sufficiently 
simple to fully evaluate the total system design. It can be readily 
expanded to a two channel system, terminating in either a single micro- 
processor or two micro-processors. The expansion simply requires that the 
controller signals be routed to two channels and that simple modifications 
be incorporated into the software whether one or two micro-processors 
be required. 

5.2 An Overview of the Operation of the System 

As indicated above, the scatterometer processing system consists of 
two major subsystems together with the necessary interfaces to permit 
communication with various peripheral devices. The relationships among 
the major subsystems and interface units was shown in Figure 5.2. An 
alignment generator, although not required during operation, was also 
shown. 

When the system is powered, software control of the system is 
assumed by the micro-processor within the so-called RESET mode. Within 



figure 5.2 ARCHITECTURE OF AN ENGINEERING MODEL 
EXPANDABLE TO THE TARGET SYSTEM 









this mode the aircraft, operator and output interface units, and the 
CZT filtering unit within the PSD estimation subsystem are active (see 
Figure 5.2). When the system is within this mode, the operator may over 
ride certain aircraft parameters arriving at the aircraft interface 
should they be inaccurate or missing in the aircraft data stream. From 
the RESET mode, scatterometer data may be processed on a flight line by 
entering the RUN mode. 

Within the RUN mode, the software controller determines the number 
of 512 point sub-records to be processed to form an average return from 
the 512 beam resolution cells in view by the CZT filters. The number of 
sub-records N^ and the address to which the accumulated reserves are to 
be stored are transferred to the PSD controller (see Figure 5.2). Upon 
this initialization, the accumulator and detector are activated at the 
beginning of a CZT cycle. The PSD estimation subsystem then accumulates 
Np estimates in each spectral channel. At the end of the Np sub-records 
the detector and accumulator are halted and a DMA (direct memory access) 
transfer is made from the accumulator to the micro-processor memory. If 
the processor has not been halted by the operator, the detection and 
accumulation cycle is re-initiated by the micro-processor once the DMA 
has been completed. It should be noted that the filtering sub-section 
operates continuously between accumulations to avoid start up transients 

A closer look at the PSD estimation subsystem will indicate that 
this subsystem filters the fore and aft Doppler spectra of the radar 
return to form estimates of the return power is 512 parallel spectral 
channels. Of these channels, 256 are available to characterize the 
fore spectrum and 255 the aft spectrum. The remaining channel monitors 
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the return from the nadir point (this return is suppressed by the radar 
and the processor). A sliding CZT algorithm, as discussed in Section 3.5, 
is employed to form repeated estimates of the power return within each 
channel. Several estimates are summed to form an improved estimate of 
the average power. 

The format of the summed spectral estimates in the 512 channels is 
illustrated in Figure 5.3. When the forward CZT is formed on the complex 
signal x + jy, the aft spectrum appears in the first 257 accumulation bins 
(channels) and the fore spectrum in reverse order in the latter 255 
accumulation bins. Some of the spectral channels among the 512 channels 
are reserved for the calibration and polarization tones. The effective 
bandwidth of each channel is given by B = f5/512 where f^ is the sampling 
frequency. The frequency resolution, i.e., the separation between spectral 
estimates is also equal to B. The normalized frequency response of each 
filter in dB is illustrated in Figure 5.4. The filter efficiency to the 
first sidelobes is 90.3% and through the first sidelobes is 95.5%*. The 
effective bandwidth is also equivalent to the width of the mainlobe. 

The average spectral estimates together with calibration and polari- 
zation tone levels are transferred from the PSD estimation subsystem to 
the micro-processing subsystem through a DMA process initiated by the 

*It is possible to increase the mainlobe efficiency by weighting 
the tap points on the CZT transversal filter. Such a device is available. 
However, this design is not recommended since 1) the precision will be 
reduced by a factor of two (the adjacent channel outputs are highly 
correlated) and 2) the S/N ratio referred to the output of the weighted 
transversal filter degrades by a factor of two (the noise is primarily 
governed by the post amplifier and signal output is reduced when 
weighting is used). 
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PSD controller. Once the detector and accumulator are re-1nitiated, the 
micro-processor, using data from the aircraft interface together with 
invariant calibration constants stored in ROMs (Read Only Memories), 
converts the spectral and calibration data to a° estimates at the required 
eight viewing angles. Adjacent spectral estimates are summed about each 
viewing angle to form a total return within the desired angular resolution. 
The spectral lines are chosen dynamically to track the specified viewing 
angle, 

The processed scatterometer data and other parameters are stored in 
an array of the micro-processor memory for each of the viewing angles. 

The 0 ° estimates at the eight angles are stored in memory in an askewed 
fashion to provide near collocation of the returns on a single cell. 

When the re-ordered data is available on a single cell at all angles, 
it is written out to the output interface which in turn transfers it to 
magnetic tape. The transfer of the data continues until the operator 
halts the processor. When halted, the processor places an end of file 
indicator within the output array, transfers the partially filled array 
to tape and enters the RESET mode to await instructions from the operator. 
At the beginning of each accumulation cycle, appropriate parameters are 
withdrawn from the aircraft data channel for use in the computations 
and for transmittal to the output tape. 


6.0 THE PSD ESTIMATION SUBSYSTEM 


6.1 Introduction 

As indicated in Figure 5.2 the PSD estimation subsystem consists of 
three sections: 

1) the CZT (Doppler) filtering unit 

2) the detector and accumulator unit 

3) the PSD and DMA controllers. 

Filtering in parallel channels is accomplished by performing a sliding 
chirp Z transform on a complex valued input signal in accord with the 
theory discussed in Sections 3.4 and 3.5. The technique for implementing 
the CZT is based on a CCD transversal filter and requires a configuration 
as described in Section 3.6. The filtering section employs analog and 
sampled analog techniques. The remaining sections are implemented digi- 
tally in order to achieve the necessary dynamic range and accuracy. 

6.2 The CZT Filtering Unit 

The CZT filtering unit is composed of four elements: 1) the input 

conditioning circuit, 2) the pre-chirp multipliers and summers, 3) the 
CCD transversal filter and 4) the signal extraction circuitry. Each 
element is discussed below. Detailed block diagrams for the four elements 
are found in Figures 6.1 and 6.3. A timing diagram for the associated 
clocking/gating signals is shown in Figure 6.4. The associated design 
schematics are included in a separate blue print file and part listings 
may be found in Appendix D. 
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6.2.1 Input Conditioning Circuit 


The input conditioner consists of buffers, low pass active filters 
(LPF), sample and hold circuits, and a bank of variable gain amplifiers 
as illustrated in the block diagram of Figure 6.1. 

The signal in each input channel is buffered through an ac coupled 
amplifier to totally remove dc and drift components induced by the 
scatterometer or recording system. It is important to remove these 
components since, when large, they can severely limit the dynamic range 
over which a PSD can be executed. This requirement is imposed by the 
input dynamic range to the CZT transversal filter which is somewhat 
less than + 2.0 volts. The ac coupling causes a h'ighpass filter break 
point at 88 Hz for the C band processor and at 28 Hz for the L band pro- 
cessor. Although the amplifier is ac coupled, an dc test signal may be 
injected into the buffer amplifier through a separate jack to align or 
test the CZT filtering unit. 

After buffering the input signals in each channel, the upper bandwidth 
of the signal is established by a four pole maximally flat Butterworth 
filter. The design procedure for the filter is well known and may be found, 
for example, in Chapter 4 of reference [6]. The transfer function for 
the LPF in the C band system is given by 
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FIGURE 6.1 DETAILED BLOCK DIAGRAM OF CONDITIONING 
AND PRE-CHIRPING CIRCUITS 












where % = ^*0 KHz for the L band system. The low pass filters are 
required to establish the processing bandwidth prior to sampling to 
avoid aliasing of high frequency noise components into the desired spec- 
trum. 

An examination of the association design schematics will show that 
low drift and low offset operational amplifiers (LM308A) have been selected 
for the buffer and LPF stages. A low drift characteristic is assured in 
each stage by balancing the shunt resistances into the non-inverting and 
inverting terminals of the operational amplifier. This design approach is 
important to minimize the dc component before chirping. 

With the bandpass established by the buffer and LPF, the signal is 
prepared for analog sampling. A sample and hold (S/H) circuit converts 
the filtered input into a staircase approximation. The signal is sampled 
at a frequency 3 times the highest calibration/polarization tone. Analog 
sampling is required to stablize the amplitude for the subsequent pre- 
chirp multiplication and entry into the CCD transversal filter. 

An N channel MOSFET is employed as a series switch in the S/H circuit. 
Its substrate is biased at -8.2 VDC and its gate is protected with an 
external diode (see schematics ). The sampling switches are actuated by 

A high impedance FET input operational amplifier is employed to 
buffer the sampling capacitor. 

A pair of offset and gain adjustable amplifiers are employed to 
provide two parallel channels after each sampler. Two channels are 
required in each input channel to prepare the I and Q signals for complex 
valued multiplications and additions. The dc offsets in all four channels 


should be suppressed to the millivolt level prior the pre-chirp operation. 
The gains must be adjusted to nominally provide a 10 volt signal into 
the chirp multipliers for the maximum recorder output (1.5 volts). Final 
gain adjustment should account for the gain differences through each 
multiplier, however. It is important to operate the channels at the 
highest permissible amplitude (accounting for signal fading) to produce 
multiplications with good S/N characteristics. 

6.2.2 Pre-chirp Multiplication and Summing 

As indicated in Section 3.4, the CZT requires that the complex 

2 

valued input signal z be pre-chirped with where (k-n) is the 

sample index*. Therefore the pre-chirp and summing operations must form 
2 

^ ^^m ^^m^ (cosirm^/N - simim^/N) (6.1) 

where m plays the role k-n. When the real parts and imaginary parts are 
identified. 

. 2 

Re = x^cosirm^/N + y^^ sin-rrm^/N (6.2) 

and 

. 2 

Im {z^e = -X|^simrm^/N + y^^ cosirm^/N 

The above results indicate that the I and Q channel signals must be 
multiplied by cosine and sine chirps and then algebraically summed. 

*This may not be apparent to the casual reader. However, if the 
redundances of modulo Zt\ are removed, then the assertion stands. The 
sample index here is unlimited since a sliding transform is employed. 
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The stratagem by which this is accomplished is illustrated in the block 
diagram of Figure 6.1 . 

Multiplying digital to analog converters (MDACS) are employed to 
perform the pre-multiplication operation. The MDACS, in conjunction 
with the current to voltage converters and the read only memories (ROMS), 
are required to form the products cited in equation 6.1. An MDAC is 
composed of a resistor ladder network whose shunt currents may be diverted 
in either of two output ports. When these ports operate into virtual 
grounds as provided by an operational amplifier circuit, the current at 
output port #1 is given by 




(a^ + a22"^ + . .. + a^2"”) 


and at the second port by 

h ^ ^in ■ h 

where aj^e {0,1} and k e{n} and = 10"^ V^.^. The aj^s are externally 
selectable by appropriately addressing the ROM and storing the binary 
number (a^, ag, •••Sp) at that memory location. When these currents 
are sunk into the virtual grounds of the operational amplifier network 
shown in Figure 6.2, the output voltage is given by 


''o = >1n^2 


(1-2K) 


“1 “*2 n 

where K = a^2“ + a22~ + ... + a^2'". To produce a cosine chirp (i.e., 
the must be multiplied by a cosine chirp), taps on the ladder net- 
work must be selected so that 
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1-2K^ = cosirm^/N 
or 

2 

» _ l-cosTTin /N 
■ 2 

when the mth sample is inserted. Similarly to produce a -sine chirp, 
the MDACS must be addressed so that 

2 

„ _ 1 + simrm /N 

h ~ 2 

on the mth sample. Since 0 < K < 1, it is clear that the chirps can be 

achieved. The proper values of K and K are stored sequentially in a 

ROM and withdrawn synchronously with the input samples. A FORTRAN 

program which generates the closest binary values of K and K through- 

out a single chirp cycle is shown in Table 6.1. The values of K and K 

^ s 

are printed in hexi decimal and their BCD equivalent. 

As a note of caution, it is important that the operational amplifiers 
employed in the current to voltage converters have sufficient bandwidth 
and slew rate to support sudden reversals in polarity. The LF 356 
operational amplifier specified in the design schematics meet the band- 
width and slew rate requirements. Their low noise performance also 
makes them attractive in these circuits since the signals are limited by 
the current capability of the MDACS. 

After the I and Q signals are appropriately chirped, differencing 
and summing amplifiers are employed to form the real and imaginary parts 
required by equations 6.2, and 6.3, respectively. Prior to the summing and 
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TABLE 6.1 CHIRP 6ENERATI0N PROGRAM 


differencing operations, dc offsets are introduced into the real and 
imaginary channels since the CCD transversal filter requires that the 
chirped signals be offset on a positive voltage. 


6.2.3 The Transversal Filter and Signal Extraction Circuitry 

The real and imaginary parts of the chirped signals are injected 
into the imaginary and real channels of the CCD transversal filter, 
respectively*. Two transversal filters are required on each input to 
perform the complex valued convolution required by the CZT. 

With the connections as illustrated in Figure 6.3, the quad filter 
bank is matched to a down chirp. The reversal in sign in one of the 
sine channels is purposely included in the bank so that both of the 
output channels can be differentially connected. A differential connection 
is required because the convolution outputs must be extracted by pulsing 
the (j)| and (j)^ and the and outputs through coupling capacitors 
(see Figure 6.3). Once a differential output is pulsed, the capacitors 
act as integrating capacitors to receive the output charges from the trans- 
versal filters. Prior to receiving each new convolved output, the previous 
output sample is discharged using MOSFET switches within the CCD chip. 

The switches are activated by the and gating signals. Once the 

+ 

integrating capacitors have been discharged, the output on the 
channel is extracted with gating signal <}'■] and the output on the (j)^- 
channel is extracted with gating signal (t)^. The combination of signals 

*The CCD filter is masked to perform the inverse DFT. This requires 
that the inputs be t>^ansposed to produce the forward transform as required 
by this design. 


67 


T3R TIR 



FIGURE 6.3 DETAILED BLOCK DIAGRAM OF TRANSVERSAL FILTER 
AND SIGNAL EXTRACTION CIRCUITRY 









(j)^ , ‘J> 3 » ‘l >4 3*^® required to move the input signals down the 

transversal filters. The timing diagram for these gating signals is 
shown in Figure 6.4. 

An overview of the sampling and convolution operations can be 
established from the timing diagram. A basic sample period is broken 
into 64 sub-periods and the different segments of sub-periods are 
allocated to the timing signals. A new input sample is taken when 
(j)p 2 is high. This same signal discharges the previous (f)^ output signals 
on Cq. This input sample is pre-multi plied by the next chirp value 
clocked from the ROMS by CZTCLK. This pre-multi plied sample is loaded 
into the transversal filter with (Jig. Simultaneously 1^2 ®^^o latches the 
previous convolved output on the (|)2 channel into an analog output buffer. 
This output is extracted when goes high. At the end of the (l)^ 
signal, the previous cj)^ sample on the channels is discharged when 
(}ip^ goes high. Then (j)^ latches the newly convolved output on the (})i 
channel in its output buffer. The cj)^ output is extracted when (J)i goes 
high. The cycle then repeats. 

It is important to note that the transversal filter is actually a 
filter matched to a down chirp. Since all frequency components are 
down chirped in the pre-chirp operation, each frequency component is 
compressed by a factor 512 irt the matched filter. This implies that 
frequency components not "visible" on the input when chirped and com- 
pressed will be visible at the output of the transversal filter. As 
a result of this processing gain (27dB), the dynamic range of the 
transversal filter is quite impressive and is primarily limited by the 
noise of the post amplifier [7]. 
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FIGURE 6.4 TIMING DIAGRAM FOR CZT FILTERING UNIT 












The two output channels except for a constant phase factor contain 
the real and imaginary estimates of the DFT of the complex valued input 
signal. The output voltage is held and buffered by a sampling differen- 
tial amplifier. The integrating capacitors must be well balanced to re- 
tain the common mode rejection inherent in the differential amplifier. 
This assures that the extraction signal is rejected from the output of 
the amplifier. 

It is important to employ low noise operational amplifiers at the 
outputs of the transversal filter to realize the inherent dynamic range 
of the filter. Since low noise amplifiers typically have large band- 
widths associated with them, the differential stages are followed by 
low pass amplifiers to reduce the noise bandwidth, without affecting 
the signal bandwidth. The low pass amplifiers are followed by a series 
of S/H circuits. The cj)^ channel has an additional S/H to bring that 
channel into synchronism with the channel. The S/H stages also 
function to remove the residual common mode signal. This signal would 
appear at the output of the S/H chain as a dc offset if no dc compensa- 
tion were employed. However, to remain calibrated at the output a dc 
stabilization loop has been included to maintain the average output 
near zero volts. A null adjustment is available in each output channel 
to achieve this condition. The negative feedback is very effective in 
removing and stabilizing the dc component before digital conversion 
and allows the output channels to operate over a full dynamic range of 
+ 10 volts. 
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The output of each S/H chain is digitized to 12 bits by an analog 
to digital converter (ADC). The 12 bits allow a dynamic range of 41 dB 
with a resolution of + 0.5 dB at the lower end of the dynamic range based 
on a single record.* The output of the ADC appears as complementary 
offset binary but is converted to 2's complement by inverting the lower 
11 bits (see Figure 6.5). 

6.3 The Detection and Accumulation Unit 

The digitized output D03 and DOl , except for a post chirp phase 
factor, represent the real and imaginary parts of the DFT of the complex 
valued input signal. To convert these to a PSD, the real and imaginary 
parts must be squared and repeated estimates of the same spectral line 
averaged. These operations are equivalent to square law detection and 
accumulation of each spectral line. Within the PSD estimation subsystem, 
detection and accumulation are accomplished digitally in a manner 
illustrated by the block diagram of Figure 6.5. 

The real and imaginary parts are squared with 12 bit multipliers 
to yield a 23 bit result. The squares are summed and passed into an 
accumulation loop. At the beginning of an accumulation cycle, previously 
accumulated results are withdrawn from memory, ANDed with zeros, and added 
to new spectral estimates. The ANDing operation results in adding zeros 
to the first of spectral estimates and therefore initializes the accumula- 
tion cycle. The results of the addition are latched in a 32 bit write 
register, routed through a read/write switch and stored sequentially in 

*Note the total dynamic range is 66 dB; however, at low signal 
levels the experimenter cannot resolve + 0.5 dB over this range. 

However, some of this larger range is required for signal fluctuations. 
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FIGURE 6.5 DETAILED BLOCK DIAGRAM OF THE DETECTION AND ACCUMULATION CIRCUITRY 











memory. The next sample to be updated is then withdrawn from memory, 
directed through the read/write switch into the read latch. The output 
of the latch is ANDed with zeros until the initialization cycle has been 
completed. 

When additional sub-records are to be accumulated beyond the initial 
sub-record, the sequence of operations remains the same; however, the 
previous spectral estimate is ANDed with ones instead of zeros and then 
added to the new spectral estimate. This operation continues until the 
desired number of sub-records have been accumulated. The accumulation 
loop will permit accumulations up to 32 bits. For a homogeneous target 
this corresponds to averaging up to 256 sub-records. 

When the accumulation over the entire record has been completed, the 
contents of the 2x16x512 memory are unloaded sequentially under request 
by the DMA controller. The words from memory are transferred through the 
read/write switch to the read channel. Byte selection signals <fcA0 and 
(t)Al decompose the word into sequential bytes, which are then latched into 
the DMA buffer whenever the data bus is not busy. 

The sequence of operations for the detection/accumulation processor 
may be best understood from the timing diagram of Figure 6.6. During 
the 64 state cylce, the next spectral estimate is sampled when F/H goes 
low. An A to D conversion is initiated on the negative transition of the 
START CONVERT command. The digital sample is loaded into input register 
of a multiplier by the Command LDIR. The product is loaded into the 
output register of the multiplier by LDOR. DIEN directs the contents of 
memory location A into the read loop. The word is latched by RLAT so as 


START CONVERT 



TIMING DIAGRAM FOR THE DETECTION AND ACCUMULATION 





to permit ANDing with zeros or ones. The accumulation is performed 
between RLAT and WLAT and the result is latched by WLAT. The read/write 
switch is directed to write when DIEN goes low. The memory output 
function is disabled by OD and a write occurs when R/W goes low. 

6.4 The Processor Controller Subsystem 

The processor controller subsystem is shown in Figure 6.7. The 
subsystem can be further divided into a PSD controller and a DMA con- 
troller. The PSD controller provides all of the timing and control sig- 
nals for the analog chirp Z-transform circuit and the digital accumula- 
tion circuit. The DMA controller provides for rapid transfer of the 
PSD data from the digital accumulator memory to the main processor memory. 
Detailed descriptions of the operation of the PSD and DMA controller are 
presented in the following sections. 

6.4.1 The PSD Controller 

The PSD controller network performs two interrelated functions. The 
network generates a set of timing and control waveforms for the analog 
CZT board and a set of timing and control signals for the digital accumu- 
lation board. The PSD controller is implemented as a 64 state subsystem 
that is driven by the signal CCLK, the microprocessor system clock. 

The CCLK signal is divided by 16 to provide a signal called MCLK. 

The frequency of MCLK is 576 KHz which provides a minimum state time 
interval 1.74 ps. MCLK increments a 6 bit counter whose outputs provide 
addre! ng for the 64 state controller read only memory (ROM). There are 
16 state control) signals generated by the controller ROM. 
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PSO CONTROLLER 



FIGURE 6.7 PROCESSOR CONTROLLER SUBSYSTEM 


The state control signals are buffered by a data latch which is 


driven by MCLK. Buffering is provided to eliminate transition noise from 
the ROM outputs and to assure that all signal transitions are in exact 
time coincidence. Of the 16 state control signals, 6 are clock signals 
which operate the CCD transversal filter, 9 are required for A/D conver- 
sion, squaring, and digital accumulation, and one, CZT CLK, is utilized 
by both analog and digital circuits. 

The 6 analog control signals are identified in Figure 6.8 as TICL- 
T4CL, TIR and T3R. These signals are supplied to the CZT board, which 
in turn, generates the signals <})^ - <})pi and described in sec- 

tion 6.2. The signal CZT CLK increments a 9 bit chirp ROM address counter. 
The resulting address signals, RAO - RA8, select the appropriate chirp 
waveform values to be multiplied with the sampled input signals as the 
first step in performing the discrete Fourier transform. 

The digital control signals are also shown in Figure 6.8. The 
timing relationship for the signals are depicted in Figure 6.6. With 
the exception of S/TT, the digital control signals are shown to be gated 
off during the DMA transfer operation. 

There are 9 address lines required to address the digital accumula- 
tion network memory. These address lines are also generated as the output 
lines of a 9 bit counter as depicted in Figure 6.8. The accumulator 
address memory counter is increment by one of two possible clocking signals 
During data sampling and accumulation the memory address counters are 
clocked in synchronism with the CZT filter by the gated CZT CLK signal. 
During the data transfer operation the address counter is clocked at a 
much more rapid rate as determined by the signal INCP. 


■ .... ... . 
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state czt sample control 

CONTROL! 1 



FIGURE 6.8 THE PSD CONTROLLER 










During DMA only a single 8 bit byte can be transferred at a time so 
that each 32 bit data word must be transferred as four 8 bit data bytes. 

To identify the appropriate byte within each data word, two additional 
address, or byte select, lines are required. The byte select lines 
change with each cycle of INCP, whereas INCP must be divided by four to 
provide the clock signal drive for the memory address counter. 

The transfer of the control of the accumulator memory and address 
counter is implemented by means of a record counter followed by two flip 
flops. These circuits generate an interrupt signal, I NT 1 , that signals 
the end of a data accumulation cycle and that the DMA operation may begin. 

To begin an accumulation cycle, an OUT instruction from the micro- 
process preloads a counter with the number of records that are to be 
added together or averaged. Upon the LOADN command the accumulation 
process begins. Each rising edge of the most significant bit of the 
address memory counter, A8, decrements the record counter. When the count 
equals zero, the next rising edge of A8 causes the counter to generate 
a "borrow" signal. The borrow signal in turn sets the DMA request flip 
flop which the processor observes as an interrupt. The processor services 
the interrupt by means of an OUT instruction to address 09. The signal 
generated by the OUT instruction is called DMA. This signal activates 
the DMA operation. 

6.4.2 The DMA Controller 

The functions of the DMA Controller are to acquire control of the 
microprocessor multibus, transfer the accumulated PSD data to the proces- 
sor main memory and return control of the multibus to the microprocessor. 


Gaining and relinquishing control of the microprocessor multibus is 
accomplished by the circuit shown in Figure 6.9. The processor responds 
to a DMA request interrupt with the signal DMA. This action causes the 
"bus priority in" signal BPRN to go high. BPRN prevents the microprocessor 
CPU from accessing the multibus during the data transfer operation. 

In order for the DMA controller to transfer the data, the CPU must 
first relinquish control of the bus. This event is signaled to the DMA 
controller by allowing the open collector signal BUSY to go high. The 
DMA controller seizes control of the bus with the rising edge of BUSY, 
and indicates its control by asserting BUSY back low. At the conclusion 
of the DMA transfer, the signal labeled DMA COMPLETE initiates action by 
the controller to return control of the multibus to the CPU. 

In addition to initiating bus control transfer, the signal DMA in- 
hibits the accumulation control signals to the digital accumulator board 
and forces the signal, OD'GATE, low, OD GATE controls the digital accumula- 
tor memory data bus multiplexor such that data may be read from the accumula- 
tor memory and it also provides a means for synchronizing the accumulation 
cycle with the stai"t of the chirp filter cycle whenever a new block of data 
is to be taken. 

Once control of the multibus is established, the circuits in Figures 
6.10 and 6.11 provide for the actual data transfer. The circuit in Figure 
6.10 is driven by the computer clock signal CC'LK. During a DMA transfer 
cycle BPRN is high and the Q output of the transfer inhibit flip flop is 
also high. Under these conditions the clock signal, CCLK, propagates 
through the three input AND gate, is divided by 8, and subsequently clocks 
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FIGURE 6.11 DMA CONTROLLER ADDRESS GENERATION 





a D and J-K flip flop, both connected in a toggle mode. The D flip 
flop toggles on a rising edge of the divided clock signal whereas the 
J-K flip flop toggles on the falling edge. 

The divide by 8 circuit first generates a falling edge clock signal 
so that the J-K flip flop toggles before the D flip flop. The Q output 
of the J-K flip flop is designated as the signal INCP. INCP increments 
the byte select segment of the accumulator memory address register, and 
on every fourth transition, increments the memory address. 

After the proper memory address is provided, the D flip flop is 
toggled and a WRITE control signal is generated which commands the 
processor memory to store the PSD data present on the data bus at the 
memory location specified in the DMA controller address register. The 
complement of WRITE clocks the transfer inhibit flip flop to prevent 
the CCLK signal from incrementing to the next accumulator memory loca- 
tion until the write operation has been acknowledged by the microproces- 
sor main memory controller. 

The memory controller acknowledges the data transfer via the signal 
BXACK . Upon receiving the transfer acknowledge the controller is again 
activated and the process continues until all 512 PSD samples are trans- 
ferred. When the last sample is transferred, the signal DMA COMPLETE is 
generated and the DMA controller releases control of the multibus and is 
ready to begin a new DMA transfer cycle when a new record of data is 
received. 

The circuit shown in Figure 6.11 provides main memory addressing 
during the DMA operation. Two 8 bit data latches accept the upper and 


lower bytes of the 16 bit address of the destination memory locations 
for the first PSD data byte. The address held in the data latches is 
fully software selectable. 

The output of the data latches serve to preload a memory address 
counter which is incremented by the signal INCP. The counter outputs 
are buffered and placed onto the multibus when the signal TBUSY goes 
low. 


6.5 The Test and Evaluation of the PSD Estimation Subsystem 
6.5.1 Test System 

The PSD estimation subsystem was constructed, debugged and then 
evaluated. The schematic diagrams reflect the design revision identi- 
fied in the debugging process. The subsystem was evaluated using a 
test configuration shown in Figure 6.12. The PSD estimation subsystem 
was installed in and powered from the INTEL SBC-660 system chassis. 

Also included in the chassis was an INTEL SBC-80/20 single board computer 
and a SBC-116 combination memory and I/O board. 

As implied by the diagram of Figure 6.12 the micro -computer permit- 
ted communication with the PSD, estimation subsystem and other peripheral 
elements. An operating system was developed for the 80/20 microprocessor 
to permit interaction with 1) the Silent 700 terminal to include tape 
playback and record, 2) the PSD estimation subsystem and 3) the TI 
980A computer. Controlling software was developed and recorded on 
cassette tape. The Silent 700 terminal was employed to command the 
system operation and to enter the appropriate controlling software 
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FIGURE 6.12 PSD ESTIMATION SUBSYSTEM TEST SET-UP 






from tape. The software permitted interactive control from the Silent 
700 and provided sufficient flexibility to display spectral data in 
hexidecimal on the local Silent 700 terminal or to transmit it to a 
remote TI 980A Computer for conversion to floating point and graphical 
display. 

Various signals were injected into the test system to document the 
following parameters: 

1) dynamic range 

2) power linearity 

3) frequency envelope 

4) fore and aft signal separation 

Other test results of interest were also generated. Prior to each 
sequence of test runs the CZT filtering unit was aligned in accord with 
the procedure outlined in Appendix C. 

6.5.2 Test and other Illustrated Results 

To test for the dynamic range and linearity of the PSD estimation 
system, a 2000 Hz tone was injected into the x channel with the input to 
the y channel grounded. The amplitude of the tone was varied over a 45 
dB dynamic range. At each amplitude, 9 records were filtered, accumulated 
and displayed on the VERSATEC printer. Typical spectral plots for two 
different input amplitudes are illustrated in Figures 6.13 and 6.14. 

The peak spectial amplitudes were correlated with the input voltage 
represented in dB to demonstrate both the dynamic range and linearity 
of the subsystem. The result is shown in Figure 6.15. From this graph 
it is apparent that the subsystem exhibits a linear input power - out- 
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GURE 6.14 






put power relationship over a 45 dB dynamic range. Since the noise floor 
is more than 10 dB beneath the peak spectial line for an input of -30 dB 
(see Figure 6.14), the actual dynamic range is somewhat greater than 50 
dB with linearity undoubtedly occurring over this entire range. 

The actual frequency response of a single Doppler filter (one of 
512 parallel filters) was not measured because of the difficulty in 
maintaining a known frequency difference between the source and the 
Doppler filter center frequency. Visual experience from oscilloscope 
traces has, however, demonstrated that the unweighted RETICON filter 
has a sin X/X like response since when the source frequency matches a 
filter center frequency the side lobes essentially disappear. Similarly, 
where the Hanning weighted RETICON filter is employed, the first side 
lobes are 6 dB down when source and filter frequencies are aligned. The 
filter response envelope of the unweighted filter can be inferred approxi- 
mately from the test data taken in the dynamic range and linearity test. 
Since the source and filter center frequency were misaligned, the adjacent 
spectial data represents the approximate response envelope of a single 
filter. The actual response will be within a few dB of that actually 
observed. The results of this test are illustrated in Figure 6,16. The 
graph represents a high resolution look of spectial lines of the type 
shown in Figures 6.13 and 6.14. 

Since the NASA fan beam scatterometers illuminate the fore and 
aft track sectors (see Figure 3.1), both up and down Doppler spectra 
appear in the return. A well designed Doppler discrimination system must 
separate the fore and aft spectial data to assure isolation between 
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FIGURE 6.16 THE APPROXIMATE FREQUENCY ENVELOPE 


returns representing two azimuthal aspects of the distributive target. 

Test were conducted to determine whether the quadrature processing tech- 
nique employed here exhibits sufficient isolation between fore and aft 
spectra. The isolation is largely a matter of properly balancing the 
gains in the quadrature channels. 

To demonstrate the fore/aft isolation character of the CZT filter- 
ing unit, quadrature square wave signals were injected into the x and y 
input channels and the spectral information averaged over several accumu- 
lation cycles. As demonstrated in Appendix B, the spectrum of x+jy 
should be unsymmetrical with the 1st, 5th, 9th, etc. harmonics appearing 
on one side and the 3rd, 7th, 11th, etc. harmonics appearing on the 
opposite side of dc. The results of the these tests are illustrated by 
the graph of Figures 6.17 and 6.18 for the unweighted and Hanning weighted 
CCD filters, respectively. It is important to note that the spectra 
appear differently in the graphs for the two filters since their trans- 
versal taps are ordered differently. The dc spectral line for the un- 
weighted filter appears on the left (it is actually displaced two bins to 
the right because of the delay injected by the S/H amplifiers); on the 
other hand the dc spectral line appears in the center bin (256) for the 
weighted filter (again actually 258 for the same reason). The fore/aft 
isolation may be inferred by comparing the amplitude of the 1st harmonic 
with the location in the negative spectrum where it would have occurred 
if it were a double sided spectrum. An examination of Figure 6.17 shows 
that the weighted CCD filter rejected the first harmonic in the negative 
portion of the spectrum (the righthand side) by more than 30 dB. Since 


94 







the image of the first harmonic is not apparent on the righthand side, 
the amplitude of the test signal should have been incrieased to identify 
the actual isolation. But nevertheless 30 dB or more represents more 
than adequate isolation. An examination of the same spectrum produced 
by the weighted filter shows more than 25 dB isolation under the same 
input conditions. This result was anticipated since the gain through 
the weighted filter is 3 dB less than the unweighted. 

Other test data were also generated to characterize the performance 
of the CZT filter unit on noisy signals. The graphs of Figures 6.19 
through 6.21 represent PSD estimates of colored noise as generated by 
bandpass filtering broadband noi.se. Figures 6.19 and 6.20 show the 
affects of averaging 1 and 16 records, respectively, using the unweight- 
ed filter. Figure 6.21 represents a 16 record average through the 
weighted filter. Finally, the spectrum of actual C band scatterometer 
data was generated using the quadrature filtering technique. The result 
is shown in Figure 6.22, The 3000 Hz (approx.) calibration/polarization 
tone appears in the fore spectrum as it should. It is also noted that 
the aft and fore spectra are unsymmetric since the beam is primarily 
directed aft. 

6.5.3 Discussion of Abberations in the Test Results 

From this sequence of illustrations and actual experience several 
problems were identified in the PSD estimation system; however, all of the 
problems appear solvable by proper adjustments to the design, improving 
circuit fabrication practices for the analog board and reducing the 
background noise in the laboratory and its equipment. Identification 
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FIGURE 6.20 SPECTRUM OF COLORED NOISE CENTERED ON 1000 Hz (16 RECORD AVERAGE) 










and discussions if these problems are treated in the following paragraphs. 

An examination of the graph of Figure 6.19 suggests that the quantum 
jumps in the spectra between the major noise lobes is not associated with 
the CCD device or the input noise. These constant amplitude jumps may be 
attributable to a high noise environment for the A/D converter induced by 
using wire wrapping techniques in a high density digital circuit. The 
problem is thought to be a major contribution to the high noise floor 
level in all of the spectral illustrations produced by averaging several 
records. It is suspected that these quantum jumps may be eliminated by 
mounting the A/D converter on doubly clad circuit board to shield its 
terminals from the digital and laboratory noise. 

An examination of Figure 6,23 shows a spurious spectral response on 
the right of the 2000 Hz line. The source of the spurious response is 
unknown, but it is known to occur where high amplitude signals, not exceed- 
ing the dynamic range of the CCD device, are injected into the filtering 
unit. The spurious response primarily occurs on the positive frequency 
side of the spectrum and seldom occurs over the negative frequency spec- 
trum. This property is thought to induce the lack of symmetry in the 
positive and negative spectrums of Figures 6.20 and 6.21. The fact 
that it primarily occurs over the positive frequencies and not the 
negative frequencies strongly suggests that the problem is solvable. 



7.0 THE MICRO-PROCESSING SUBSYSTEM AND INTERFACING UNITS 


7.1 Microprocessor 

7.1.1 Multibus System Description 

The signal processing system consists of the following interconnect- 
ed units: 

1) Front panel and operator interface 

2) Aircraft interface 

3) Output interface 

4) CZT filtering unit 

5) PSD and DMA controller unit 

6) Detection and accumulation unit 

7) CPU 

8 ) Memory 

It is widely accepted that the most efficient way of establishing 
the necessary communication paths between these units is through the 
use of a system interconnection bus [8]. 

The Intel Multibus was selected for this application because of its 
industry-wide support and its flexibility. The Multibus includes the 
following signal lines: 

16 address lines 
16 bi-directional data lines 
8 multi-level interrupt lines 
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The address and data lines are driven by three-state devices, while 
the interrupt and some other control lines are open-collector driven. 

Modules using the Multibus have a master-slave relationship. A 
bus master can command the bus, while bus slaves cannot. The CPU is a 
bus master in this application, while the memory is a bus slave. The 
Multibus is expandable to a maximum of 16 masters. 

The system also 'has a DMA controller which acts as a bus master. 

Bus arbitration results when the DMA controller requests use of the bus 
while the CPU controls it. The bus clock, derived from the CPU, provides 
a synchronous reference for the resolution of bus arbitration. The 
Multibus provides for a maximum data transfer rate of SMByte/second. 

System use of the multibus is found in Section 6.4 and in vendor's 
literature describing the CPU and the memory [9]. 

7.1.2 System Components 

7. 1.2.1 CPU 

Early in this contract the 8080 family of microcomputer was selected 
as the unit best suited for this application [10]. It was recommended 
that the CPU, memory and available I/O be purchased rather than constructed 
to take advantage of vendor reliability and maintenance support. Since the 
Intel SBC-80/20 was introduced in 1976^ other vendors, including 
Christian Industries, National Semiconductor, Monolithic Systems, Ad- 
vanced Micro-Computers, MuPro, and even Intel (with their SBC-80/30) 
have introduced advanced products to dramatically extend the capabilities 
of the well designed Multibus system. For example, Advanced Micro Corn- 
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puters has introduced their Multibus compatible AMC 95/4000. Both the 
SBC-80/20 and the 95/4000 fulfill system requirements for 8080 compatibil- 
ity, on board RAM and EPROM, and programmable interrupt structure (see 
Figure 7.1). But the 95/4000 has two features which make it uniquely 
superior to the 80/ 20. 

The first is processor speed. The 95/4000 comes in two versions, 
the faster of which contains a 4 MHz 8080. This means that the 95/4000 
can execute instructions twice as fast as the SBC-80/20. The second 
superior feature of the 95/4000 is its on-board Ani9511 Arithmetic Process- 
ing Unit, It allows for the evaluation of ordinary and transcendental 
mathematical functions up to 100 times faster than software [11]. 

The advantages of the 95/4000 over the 80/20 may be quantified by 
examining the effect of the faster processor speed and arithmetic processor 
on program execution speed and memory requirements. 

Program execution time is divided between mathematical software and 
control software. Based on an analysis of the 13.3 GHz scatterometer 
processor software, which is assumed to be similar to that of the proposed 
C- and L-band processors, it was found that the processor spends about 
93 % of its time doing mathematical software and only 7 % doing control 
functions. If T is the total execution time for the program this may 
be written as 

T = 0.93T + 0.07T (7.1) 

where 0.93T is the time due to the execution of the mathematical software 
and 0.07T is the time due to control functions. If these times could 
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80/20 


95/4000 


Word size 

8 bits 

8 bits 

On board RAM 

2K bytes 

2K bytes 

On board EPROM 

0/8K bytes 

0/12K bytes 

Addressable memory 

64K bytes 

64K bytes 

Clock frequency 

2 MHz 

2/4 MHz 

I/O parallel 

48 lines 

48 lines 

I/O serial 

1 RS-232C 

1 RS-232C a: 

Arithmetic capability No* 

Yes 

DMA capability 

No** 

Yes 

Priority interrupt 

Yes 

Yes 

♦available with the 

addition of 

the SBC-310 

♦♦available with the 

addition of 

the SBC-501 


FIGURE 7.1 COMPARISON OF SBC-80/20 and AMC 95/4000 


be decreased independently by a and 3, respectively, then the new exe- 
cution time is 


Tj^Ew 0.93aT + 0.073T 


(7.2) 


where a 


time to do math at increased speed 

time to do math software without increased speed 


and 


_ time to do control software at increased speed 

time to do control software without increased speed. 


Using the 95/4000, a = 0.109 and 3 = 0.667 (see section 7. 1.2. 2). As 
a result, the new execution speed is 

Tf^Ew = (0.93)(0.109)T + (0.07) (0.667)T = 0.148T (7.3) 

New execution time is seen to be 85% less than the old. This means an 
increase in execution speed by a factor of almost 7. Concern has been 
expressed that the 4 MHz 95/4000 would be difficult to use with a micro- 
processor development system due to its increased clock rate. If we 
use the 2 MHz 95/4000, then 3 = 1 .0 and a remains at 0.109. The impact 
on execution speed of using the 2 MHz 95/4000 is then 

Tj^Ew = (0-93) (0.109)T + (0.07)(1.0)T = 0.171T (7.4) 

This corresponds to a decrease in execution time of 83% or an increase 
in execution speed by a factor of almost 6. The conclusion is that 
even with a "slow" 95/4000, execution speed is increased by almost 600% 
due to the large increase in speed of the mathematical functions. 


Similarly, memory usage may be divided between that used to store 
the mathematical software and that used to store the control software 
and data. For the 13.3 GHz processor, this may be written as 

M = 0.678M + 0.322M (7.5) 

where 0.678M is storage for non-mathematical software and 0.322M is 
memory for the mathematical subroutines. The only portion which we of 
a 95/4000 would decrease is the mathematical routines. The new equation 
is 

%W “ + 0.322yM (7.6) 

where y = yequiy’ed for math routines without 95/4000 

memory required for math routines with 95/4000 

The analysis show that y can be as small as 0.116. Thus, the new memory 
requirement is 

%W " (0.322)(0.116)M = 0.715M 

The use of the 95/4000 decreases memory requirements by 28.5%. The 
conclusion is that the arithmetic processor makes a significant differ- 
ence in both execution speed and memory requirement. 

When the above analysis are considered, it is concluded that 
because the arithmetic processor offers such a significant advantage, 
it would enhance the capabilities of any Multibus compatible processor 
module. In view of the experience gained with the SBC-80/20 and of the 
relative costs ($1430 for the 80/20 plus the AMC 95/6011 APU versus 
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$1390 for the 2 MHz 95/4000) it was decided that the purchase of the 
95/6011 APU and the SBC-80/20 would best satisfy present and anticipated 
system requirements. 

A complete description of the Intel SBC-80/20 and the AMC 95/6011 
is found in the vendor's literature for each part. 

7. 1.2. 2 Memory 

The software system required 50K bytes of RAM and 14K bytes of EPROM 
for its operation. Since the 80/20 provides 2K bytes of RAM and sockets 
for up to 4K of EPROM, an addition of 48K bytes of RAM and sockets for 
lOK of EPROM is required. Intel provides 48K of Multibus compatible RAM 
with their SBC-416. The memory map of the processing system as imple- 
mented is shown in Figure 7.2. A complete description of the SBC-048 
memory is found in the vendor's literature. 

7. 1.2. 3 Front Panel and Operator Interface 

The front panel is interfaced as an I/O device to the microprocessor 
subsystem. It is wire-wrapped and therefore requires two card slots. 

Its design and operation are described in section 7.2. 

7. 1.2. 4 Aircraft Interface 

The interface to the aircraft NERDAS computer is constructed on a 
wire-wrapped card and requires two slots. It is fully described in 
section 7.3. 

7. 1.2. 5 Output Interface 

The interface to the aircraft Bi (})-L tape recorder is constructed on a 
wire-wrapped card and requires two slots. It is documented in section 7.4. 
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A 

/> 


Address 


Location 



0000 

OFFF 

1000 


37FF 

3800 


3FFF 

4000 


FFFF 


EPROM on SBC-80/20 


EPROM on SBC-416 


RAM on SBC-80/20 


RAM on SBC-048 


FIGURE 7.2 SYSTEM MEMORY ORGANIZATION 
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7. 1.2, 6 PSD Estimation Sub-system 

The PSD Estimation Sub-system is constructed on three wire-wrapped 
circuit cards and requires six slots. It is described in section 6.4. 

7.1.3 System Chassis and Power Requirements 

7. 1.3.1 System Chassis 

The system as described above requires 13 card slots. The MDS-800 
chassis from Intel was first considered as a potential system chassis. 
However, although the MDS-800 has enough slots to contain the single- 
polarization development system, a dual -polarization system would require 
more than 16 slots. Thus, a MDS-800 system is not upward compatible. 

To satisfy potential system expansion requirements, an expandable chassis 
consisting of intially one SBC-604 and three SBC-614's (all from Intel) 
provides 16 slots of chassis space, expandable to 20 by the addition of 
a single additional SBC-614. A tentative card slot arrangement is given 
in Figure 7.3. Information on the SBC-604 and -614 is found in the 
vendor's literature. 

7. 1.3. 2 System Power 

The system power requirement is given in Figure 7.4 and is satisfied 
through the use of properly sized fixed voltage power supplies. However, 
greater flexibility allowing the adjustment of at.^log voltages to minimize 
noise is provided if the supplies for items (4) and (5) are made variable. 
Connection of the power supplies to the system bus is as follows: 








Voltage Current Requirement 


+5 

<30 A 

-5 

200 mA 

+12 

360 mA 

-12 

25 mA 

+16 

350 mA 

-16 

200 mA 

+20 

25 mA 


*Note: 30A is rated output of power supply; actual 

requirement is less than 30A. 


FIGURE 7.4 POWER REQUIREMENTS 
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Multibus 

Pin 

Definition 

5, 6, 81, 

82, 83, 84 

+5 

7, 8 


+12 

79, 80 


-12 

28 


+16 

30 


-16 

34 


+20 

11, 12, 27, 29, 32, 75, 

76, 85, 86 GND 


7.2 Front Panel and Operator Interface 
7.2.1 Design 

The contract requires an integral front panel allowing the operator 
to over-ride NERDAS aircraft parameters in the event of aircraft equip- 
ment failure during a mission. A design meeting this requirement is 
shown in Figure 7.5. The central part of the design is two Matrox 
MTX-Al programmable keyboard/display controller chips [12], The front 
panel (Figure 7.6) contains a 16 key keyboard for the entry of over-ride 
data and ^ keys for control of machine state (see section 8). The top 
11 display positions are for displaying the over-ridden parameters and 
their values. Over each of the bottom nine displays is a legend. The 
four on the left siiow an when the calibrate signal is absent, when 
the system's data buffer is filled, or when the integration time is 
inadequate. On the right most five displays, a blank is displayed if 
a parameter is not over-ridden; otherwise, an is displayed under 
the legend corresponding to the over-ridden parameter. Complete parts 
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lists, vendors, and schematics are given in the appendices. The design 
is otherwise self-explanatory. 

7.2.2 Overview of Usage 

The following addresses are used by the microprocessor for control 
of the CZT hardware and the front panel display. 

Address (hex) Function 


08 

LOADN/ 

09 

DMA/ 

OA 

LOADLSB/ 

OB 

LOADMSB/ 

OC 

DISPIA/ 

OD 

DISPIB/ 

OE 

DISP2A/ 

OF 

DISP2B/ 


These functions are described below in terms of their affect on the hard- 
ware. A later section details the software operation. 

DISPIA/ references the first of the two Matrox MTX-Al keyboard/ 
display controllers. Data may be written to positions 1-6 of the top 
display and positions 1-9 of the bottom display (see Figure 7.6) 

with an OUT instruction or read from the keyboard with an IN instruction. 
A read from DISPIB/ is used to load the contents of the status register 

of the first MTX-Al into the microprocessor accumulator. 


DISP2A/ references the second MTX-Al . Data may be written to 
positions 7-11 of the top display using an OUT instruction. An IN 
instruction from this address produces unpredictable results. 

DISP2B/ is used to load the contents of the status register of the 

second MTX-Al into the microprocessor accumulator. 

The operation of the panel has been optimized for user convenience. 
Suppose that the NERDAS altitude was to be over-ridden. The desired 
numeric value is entered via the keyboard, ENTER is pressed, and the 
alternate function of the 8 key is selected for ALTITUDE. The DISP 
function allows the operator to display the over-ride value for a partic 
ular parameter which has been over-ridden. The functions of the CLEAR, 
RUN, HALT, and RESET keys, as ^'oll as further details of front panel 
operation, are available in Section 8. 

7.3 Aircraft Interface 

The aircraft interface module provides the capability to extract 
the aircraft flight narameter data that is required for o° computation 
from the NERDAS Bi-tj) L bit stream. There are three basic functions 
required in the data extraction operation. The first step is to gener- 
ate a phase reference and master clock signal, from the NERDAS data 
signal. This operation is accomplished by applying the input signal 
to a full wave rectifier to remove the encoded phase information. The 
resulting signal serves as the input to a phase locked loop (PLL) cir- 
cuit which provides a stable phase reference signal at. the precise 
frequency of the input signal. 


The second step in the data extraction operation is to decode the 
incoming data stream. This decoding operation consists of detecting the 
encoded pattern of ones and zeros from the Bi-c() L signal and generating 
a corresponding TTL level data sequence. The TTL sequence is continual- 
ly examined to identify the presence of a 12 bit sync pattern that 
signifies the beginning of a frame of aircraft flight parameter data. 
When the beginning of a data frame is detected the serial bit pattern is 
synchronously accepted by one of the serial interface ports of the SBC 
80/20 microprocessor. The port converts the serial data to 8 bit 
parallel data such that each data byte corresponds to two 4 bit NERDAS 
BCD digits. The data are stored directly in the processor memory and 
are immediately available for use in the a° computation. 

By using one of the serial interface ports of the SBC 80/20 the 
design of the NERDAS interface is considerably simplified in terms of 
both hardware and software. The acquisition of NERDAS data is fully 
controlled by software and data transfer occurs only when expressly 
requested by the a° computational algorithm. 

7.4 The Bi-(j) L Interface 

The output interface of the real-time scatterometer processor is 
a 100 KHz Bi-(j) L serial data channel which will be recorded on a 14 
track analog magnetic tape unit. The output interface is very similar 
to that employed in the 13.3 GHz real-time processor [ i ]. 

The Bi-0 L output frame consists of 256 4-bit data words. The 
first three of these words are frame SYNC words. Following the frame 
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SYNC words is a frame identifier word which indicates the type of data 
contained in the remainder of the frame. A data frame may either 
contain calculated radar cross section data or it may be a fill frame 
which consists of alternating I's and O's. 

Figure 7.7 is a block diagram of the Bi-cf) L output interface. 

Three 8212 data latches contain the SYNC and frame identifier words 
which form the beginning of each data frame. Control logic is provided 
to select the proper frame identifier word at the beginning of each 
output frame. 

A universal synchronous receiver transmitter (USRT) performs the 
parallel to serial conversion required by the Bi-tj) L generator. A bit 
counter is provided to signal that 256 words have been transmitted and 
that a new frame must begin. 

Transfer of the o® data from the CPU to the Bi-cj) L interface is 
made via the SBC-80 parallel interface port. 
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'W 


... 


ENABLE SYNC 1&2 



FIGURE 7.7 BLOCK DIAGRAM OF BI-(f) L OUTPUT INTERFACE 












8.0 SYSTEM SOFTWARE 


8.1 Design Rationale 

The overall software system is designed around two major functions; 
i.e., 1) system operation and control, and 2) the scattering coefficient 

calculations. The first, system control, constitutes performance of those 
software functions related to the "HALT" and "RESET" modes. This soft- 
ware handles all operator interactions needed to specify conditions under 
which the system will "RUN". Specifically, this means identifying air- 
craft data values to be over-ridden, if any, and what override values to 
use. 

While in the "RESET" mode the software system will accept from the 
operator those specific over-ride values of the aircraft data system 
(NERDAS) to be used in the sigma-zero calculations. It will also display 
current contents of these over-ride values. 

Minimum operator interaction has been the over-all guiding criterion 
in the software design. For example, if NERDAS data does not need to be 
over-ridden the operator needs only to "RESET" the system, select the 
appropriate Receiver Polarization, and press^RUN". The software system 
will perform all other required functions necessary to generate the out- 
put data frames. 

In discussing the software functions for scattering coefficient 
calculations it is appropriate to first review the system geometric 
concepts. 
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8.1.1 System Geometry and Correction Factors 

The basic objective of the micro-processor software is to calculate, 
repeatedly, the estimated radar scattering coefficient, a , using timing 
and position data input from the aircraft data system and radar transmit/ 
receive power information from the chirp Z-transform hardware. This is 
done by implementing the radar equation in its general form 


''o _ (47r)' 
a = 


?ill _PR_ 

A • K 


( 8 . 1 ) 


where R = range to the target ground cell 
A = area of the target ground cell 
Z(f) = receiver roll-off filter correction 
PR = receiver/transmiiter power ratio 

K = system constants; e.g., ferite modulator, cable loss, etc. 

®r®t “ antenna receive/transmit gain correction. 

Several types of algorithms are available for performing the estimate 
defined in Equation 8.1, each differing generally in the way that the 
value of A, ground cell area, is defined. This system uses the Constant 
Angular Resolution ground cell as discussed in Section 4.3. As 
the name implies, this type of resolution provides for a ground cell in 
which the angle A0^^. , Figure 8.1, remains constant with viewing angle; 
i.e., A0 q^- = A0q^^-|. In this type of resolution cell the initial value 
of A0 q^- would normally be selected such that the first cell length near 
nadir, , Figure 8.1, is the minimum resolution required. Since and 
0|^^. are related by 

A0n. A0.. 

L. = H{tan(0p^. + — ^)-tan(0p^. —)} (8.2) 
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the L^. are specified for each 0 q^ selected, 
be the same for all 0^. ; specifically. 


The integration time would 


(8.3) 


The desired bandwidth for angle 0^^. is given by 
Bd^ = (2V L, cos^6p,)AH 


(8.4) 


where V = aircraft ground speed vector, 

L.j = instantaneous length of cell i, along ground track vector, 
0Qi = viewing angle, referenced to ground track vector, 

A = transmitter wave-length, and 
H = aircraft altitude. 

The desired doppler center frequency for angle 0^^ is given by 


fdc^- = (2 V sin9p^.)/A (8.5) 

Note that 9^^. is the viewing angle referenced to the ground track line; 
i.e., the intersection of fdc^ with the X-axis in Figure 8,2. Thus to 
evaluate Bd. and fdc^. the value of 9^.. must be first calculated. For a 
given look-angle, 9^^ , the values of X|^. and may be derived using the 
following coordinate transformations: 




— 

X" 


B 

Y" 

II 

0 

_ 


_-A. 


where A = H/cosi|), B = X" = X' 
r^^ = Z' tan 9|^^. . 


= (>^L^ Y' = -Htan i|) and 


( 8 . 6 ) 
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The two vectors 


and 


■x'1 


'l 0 

»1 


"x"" 


0 

r 

= 

0 cosij^ sinij^ 


Y" 

+ 

0 

Z' 


0 -sim|i cosi|< 


Z" 


H 

X 


|cos4> -sin<j) 

o" 


X' ' 


Y 

= 

sincj) 

COS(|) 

0 


Y' 


Z 


^0 

0 

1 

-J 


_Z'_ 



give the expressions 

Xii = (-H/cos(|j) tan 6^^ cos(j) + H tan sincj) 
Yii = (-H/cosi{j) tan sincfi-H tan \p cos<p 

p p ^ 

where tan 9'^^^ = -(tan tan'll!;) “cosij;, 
so that = Tan"^ /(H^ + 


(8.7a) 


(8.7b) 

( 8 . 8 ) 

(8.9) 

( 8 . 10 ) 

( 8 . 11 ) 


Notice in Equation 8.10 that when the roll angle exceeds the viewing 
angle, (}p > expression becomes indeterminate. Therefore, 

care must be taken when computing tan 0'[_i» machine results will 

be unpredictable. 

Because of this interaction between 0^^^ and ijj, minimum values are 
imposed upon the values of 0^^ and 0j^.. The minimum Doppler angle, 
®Dmirf derived by considering that the minimum allowable doppler 

center frequency, is BW^Z.; i.e., half the bandwidth of the 

viewing angle nearest nadir. Equating the two expressions (for the 
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near nadir case) 


BW^/2 = 


AH 


(8.11a) 


and 


fdc 


. 2Vs1nep„,„ 


trnn 


(8.11b) 


gives the expression 




(8.11c) 


L 1 

The solution for of course, depends on the parameter For 

example, at H == 457.3 meters and L = 20 meters ' 

Using rationale similar to the above, the minimum allowable viewing 
angle may also be derived. First, note in Figure 2 that geometric 
consideration Indicate that whenever fdc ^ BW/2 then X ^ L/2. Using 
vector expressions for the coordinate system in Figure 8.2 the value of 
X is given by 


X = X'cos(J)-Y'sin(|) 

Substituting for X' and Y' gives 
tane'L > simi; tancf) - X , 


(8. lid) 


(8. lie) 


but at X = L/2 and, using Equation 8.10, 

2 2 2 
tan 0, = K • + tan 

Lmin ^ 


(8. Ilf) 
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'*rr 




where 


K = -tant|< tan({) + (^) 


(S.llg) 


As in Equation (8.11c), the solution depends on the parameter L/2H. 

Analysis of the function in Equation (8. Ilf) indicates that if roll 
and drift are bounded by -10° £ ijj £ 10° and -10° £ (j) £ 10° then the lower 
bound of |0j^| is as shown in Figure 8.3. Since the value of the desired 
|0^| is normally £ 5° and -5° £ £ 5°, the limitations do not represent 

a large penalty in viewing angles, however, the system software must 
account for the above limitations in order to avoid Doppler bands that 
extend off the end of the filter bank. 

Note that the value of 0^ must also be properly limited. If the 
value of is always £ L/2, 0^ will automatically be correctly limited. 
These two limitations must be observed in the software if valid area and 
power ratios are to be calculated when using viewing angles such that 


> |0lI- 


The procedure recommended and used herein is to compare the absolute 
value of the expected or desired look-angle, 0|^^. , with the absolute value 
of roll, tjj, and whenever |i|i| > |, set the value of 0j^^. = for aft 

viewing angles, or generally, 0|_^. = J q -~ | ■ |ijj| . In addition, the value of X 
calculated in Equation 8.8 must be monitored and always set £ L/2. Test 
cases which have been run at TAMb' using this procedure indicate that satis- 
factory results are obtained for |ij;| £6°, |(|)| £ 8°, H = 457.3 meters, 

10 £ L £ 50 meters. Since the actual 0^^^ is given in the output, no loss 
of information occurs. 


L 
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Again, must be used to evaluate the desired doppler center fre- 
quency, fdc^ , and to evaluate the desired bandwidth, Bd^ . However, when 
evaluating G^G^, the receiver/transmitter gain, the proper angle to use 
is e'^^. given by Equation (8.10), but adjusted for actual Doppler center 
frequencies. This adjustment is needed because the actual Doppler center 
frequency is not always the expected center frequency. This shift in 
frequency occurs because a number of discrete spectral lines are used to 
make up the desired bandwidth, and the center of this set will not always 
correspond to the center of the desired bandwidth. For the same reason 
the desired and actual bandwidth are not the same, as illustrated in 
Figure 8.4. 

The number of filter indexes required to represent the desired 


bandwidth is given by the expression 

NF. = [Bd./Af + 0.5]j (8.13) 

The actual bandwidth is given by 

BA. = NF. Af (8.14) 

The center frequency index may then be calculated (for NF^ odd) by 

NC. = [fcd./Af + 0.5]j (8.14a) 

or (for NF^. even) by 

NC. = [fdc./Af]j (8.14b) 
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Discrete Filter elements 


1 Af. 


1 ’ 

1+1 


Expected Doppler, fdc 


f 




Actual Doppler, fac^. 



Af , 

Af 


i 

i+1 

i+2 




Expected Band, 

Bd. 


Actual Band, Ba^. 




FIGURE 8.4 DESIRED AND ACTUAL BANDWIDTH DIFFERENCES 
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The left most index pointer for band BA^ is given (for NF^. odd) by 


PRL. = NC. - (NF. - 1)/2 (8.14c) 

or (for NF^. even) by 

PRL. = NC. + 1 - NF./2 (8.14d) 

If NF^ is odd the actual center frequency is 

fac^. = [fdc^./Af + 0.5]j Af (8.15) 

When the value of NF^ is even the actual center frequency is 

fac. = {[fdc./Af ]j + 0.5 }Af (8.15a) 


From Figure 8.2 it can be seen that the actual center frequency 
places the actual viewing point (Xj^^. , , X|^^. ) on a different doppler 

line that expected, thus and 0^^ need to be updated before they 
can be used to evaluate exactly the cell area, range and gain. This 
may be done by first evaluating the new doppler angle along the ground 
track using actual center frequency, fac^. : 

®TDi ^ (fac.A/2V) (8.16) 

Next, the actual value of may be calculated by 

= {(H^ + Yj^.^) tan^ 0.pp.} ^ (8.17) 
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(8.18) 


so that the true viewing angle is 
e - T-an-’ ( 

9"P^i *“ Tsn X j 


and true Range is 


= H/cos0j|_^. 


(8.19) 


The true value of 0'|_^ for use in the gain tables is given by 

0'tl'. = Tan"^{-coS(j) (tan^0jl^^-tan^i|;)’'®} (8.20) 

The true viewing angle would be the appropriate value to be placed in 
the output data set with each estimated value of sigma-zero. 

The calibration tone, similar to the individual Doppler center 
frequencies, does not lie in the center of the filter element in which 
it is found in the PSD (Power Spectral Data) filter bank; therefore, 
appropriate corrections to the measured power level in that filter 
must be made. The measured power transmitted that appears in the 
calibration-tone filter must be corrected by 

PR = PRm (8.21) 

im)f 

where N = 512, number of filter elements, 

?TT 

= 2TTf /f ^ - M 4r> where 
c S c W 

M = [(f /f J512 + 0.5] + 1, 

C C S I 

f = sample frequency, 
s 

f = calibration tone frequency, and 
c 

PR = the measured power in the calibration tone filter. 

M 
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and 


W(J2) = sin ((2N-1) n/2)/sin (n/2) 

In the above discussion equations were presented which provide the 
exact solution to the radar equation for sigma-zero estimates. In an 
airborne, near real-time processor such as this implementation it is 
appropriate to consider more efficient and expeditious means of calculating 
these solutions. The only term of Equation (8.1) that promises any 
significant savings in computations is (R /A). 

Looking first at Range, it can be shown that computing Range, R, 
by the expression 

R^. = H/cos (8.22) 

as opposed to Equation (8.19), which used » the true viewing-angle, 
will result in a Range errors of up to 2.0% for l^Lil — 

_< ip £ 6° and -8° £ (|) £ 8°. These errors in the combined term (R^/A) 
can produce sigma-zero values off by as much as + 0.45db. 

Two alternative approximation equations for area were considered: 

A = LW/cos(t> cosijj (8.22a) 

and 


A = W(Y2-Yl)/cos(() 


(8.22b) 


Equation (8.22a) is a simple rectangular approximation where 

L = H {tan02 - tane^ }/cos(J) (8.22c) 

©2 = Doppler angle for fdu, upper frequency 
9^ = Doppler angle for fdL, lower frequency 

and 

W = 2R tan (B/2)/cos(i|i) (8.22d) 

B = Beamwidth referenced to the origin of the 
aircraft coordinate system 

Equation (8.22b), whose terms are defined in Reference [13], is a more 
rigorous approximation to the actual ground cell shape. Although com- 
putation time is somewhat longer for Equation (8.22b), it is considered 
to be the best alternative because it offers up to 10?^ more accuracy 
at viewing angle where |0Lil ^45° and with roll and drift excursions of 
up to + 6° and + 8° respectively. 

Recalling from Section 8. 1.1.1, Equation (8.2) must be solved eight 
times for eight values of cell length. By using small angle approxima- 
tions an equivalent expression is 

L. = z — (8.23) 

'' cos 9. .cos M 
LI 2 

which is accurate to + 2.0% for 0j^^. £45° and £ |i{;l . For ground 

cells defined by Constant Angular Resolution the value of A9, and conse- 
quently cos ^ , are constant; therefore, the cell lengths would involve 
only the calculations 


(8.24) 


H _ H A0' 
Li 


cos 9^^ 

2 

where A0' = A0/cos A0/2, and is a stored constant. 

For viewing angles several times grater than the aircraft roll 
angle (ip) other approximations can be made; e.g., 0^^ 0|_^ and for 

^ 0.0, 0|^^ in the Gain tables. However, scatterometer data 

processing done to date indicates that ip ^ 0.0 is not generally valid, 
therefore, the software algorithm presented in this design will account 
for aircraft roll angles. 


8.1.2 System Constants 


Equation (8.1) may be re-written in the form 


/s 

a 


0 

i 


,4^,3 Cl Z„(fdc,.) r''(6l,) (Pr(6o,) 

^•T -ggT^ • P4 " 


Where C|^ = cable loss 

K = calibration constant, plus other corrections as needed 
Z^(fdc^. ) = water roll-off filter correction for center frequency fdc^. 
GG(9'Li) “ antenna gain at antenna viewing angle 0^^- 
PR(0[ji) = receiver power at doppler angle 0^^ 
and 


Pj = transmitter power. 

The system constants will include the first three terms. The first two 
3 2 

constants, (4tt) /A and Cj^, remain fixed for a selected frequency. It is 
assumed that K may vary between horizontal, vertical, like and cross 
polarization combinations, therefore, four values should be reserved in 
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the software; i.e., Kyy» Kj^^, K^y. Similarly, there must be four 
segments in the gain tables, four segments of beamwidth and a roll-off 
table for a processor that will handle all polarization combinations. 

Two degree increments with linear interpolation in these tables has 
been adequate in prior usages. Normally, antenna data are not sufficiently 
accurate to warrant more refined tables. When indexing into the constant 
tables the software will use two flags, one indicating polarization of 
the transmitter and the other showing the polarization of the receiver. 
These flags will be combined to compute the address of the proper segment 
in the table. To make the address computation most efficient the tables 
should be ordered the same way for all constants; e.g.. Figure 8.5, gain 
values are ordered as GG^y, GGy^, GG^^y, and GGj^^. By setting the two 
polarization flags as indicated below it will be easy to compute the 
correct segment address into each table knowing only each starting 
address. 


Transmit/Receive Polarization 

Polarization Flag Combination 


00 

01 

10 

11 


VV 

VH 

HV 

HH 


If constant beamwidth values are used, as is currently being done with 
L-band data processing, no beamwidth tables will be required. Also, it 
is possible to use a function to compute roll-off, rather than a table, 
however, table look-up will be faster than function-evaluation, there- 
fore a table for roll-off should be included. 
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START 

ADDR. 


START 

ADDR. 


Sv 

BWvv 

Kvh 

BWvh 

><HV 

BW„v 



GGyy 

START 

GGvh 

ADDR. 

GGhv 

ROLL-OFF 




FIGURE 8.5 CONSTANT TABLE ARRANGEMENT 


140 


There are a number of other constants, such as A0', Tg, etc., 
associated with individual computation modules. These will be stored 
in the same region of core as the module itself. 

8.2 Software Architecture 

8.2.1 General Approach 

The software system design features a modular approach with the 
structure boundaries defined by specific functions that are performed. 
Simplicity of operation and minimum operator interaction have been 
the governing considerations in laying out the system. 

In general the system is defined by three specific "states" of 
operation. Fig. 8.6. These three states form the highest levels of 
the software architecture. The purpose of each operating state is to 
perform a specific task when directed by the system operator through 
entries on the system key board. 

The "RESET" state not only performs all initializations in prep- 
aration for the "RUN" state, but also performs all communications 
between the system and the operator. 

The "RUN" state performs all sigma-zero calculations using data 
acquired from the NERDAS frame and the radar signal integrator. The 
calculated data is aligned and formatted, along with appropriate 
NERDAS data, flags and alarms, and sent to the data output. 

The "HALT" state is simply a standby state for the processor. Once 
it enters this state and completes its close out functions, the system 
idles until a "RESET" command is received from the operator through the 
keyboard. 







8.2.2 Operator Interface 

All operator commands are entered through the keyboard and are 
received and processed by the "RESET" module as program interrupts. The 
software system Interrupt Processor, a part of the RESET module, is 
designed to respond only to those interrupts which are appropriate to 
the current "state" in which the system is operating, '^ir example, while 
in the RESET state, the integration timer (Interrupt Nb> . 1) is ignored, 
all others shown in Figu^ 8.7 are processed. 

The keyboard may be viewed as a set of three types of keys; i.e.. 


Key Type 

Key Name 

State 

: Reset, Run, Halt 

Function 

: Enter, Display, Clear 

Digit 

: Values 0, 1 , ... ,9. ' . ' 

CHS (change sign) and 
Vel , Alt, Drift, Roll , 


Pitch 


Any time a DIGIT key is pressed by the operator it creates an Interrupt- 
0, Figure 8.7. All such interrupts are processed by the RDKEY routine, 
a part of the RESET module. This type of interrupt will be processed 
only while in the RESET state, therefore, all aircraft parameter over- 
ride values must be entered while in this state. FUNCTION keys also 
create interrupts which are processed only during operation in the 
RESET state. These are Interrupts-4, -5, and -7 which are used in 
conjunction with the DIGIT keys to enter values of aircraft over-ride 
parameters and to display or clear specific values already entered. 

The operator procedure for entering over-ride values must have 
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FIGURE 8.7 VECTORED INTERRUPT CONTROL STRUCTURE 
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fjf; ... 



three properly sequenced phases; i.e., value, command, destination. 

For example, to over-ride the NERDAS values of altitude, the following 
sequence would be used: 

a. key in the string '1500', (value) 

b. key 'ENTER' (command) 

c. key 'ALT' (destination) 

This sequence results in the over-ride value of 1500 feet being used in 
all sigma-0 calculations rather than the NERDAS values. 

Depressing STATE keys such as 'RUN', 'HALT', or 'RESET' cause the 
software to process Interrupt-2, -3, or -6, respectively. As illustrated 
in Figure 8.7, the Interrupt-2 causes the processor to begin calculat- 
ing sigma-zero values. The system will remain in this state as long as 
power is supplied, and will process only Interrupt-1 or Interrupt-2 
(the HALT command). Thus any changes to the over-ride structure must 
be made by bringing the system first to HALT then RESET if already 
operating in the RUN state. 

Prior to putting the system into the "RUN" state the operator 
should select the appropriate position on the Receiver Polarization 
switch. 

8.2.3 Continuous Product Concept 

Once the system is placed in the 'RUN' state, values of sigma-zero, 
based upon input NERDAS data, are calculated continuously for each of 
a set of eight angles. As each set of eight sigma-zero values is calcu- 
lated, a set is transferred to the Bi-c() L output port for tape 
recording. As long as valid NERDAS and radar data are received, sigma- 
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A 


•RUN' INTERRUPT 



EOC INTERRUPT 

N/ 


SIGMA2 


WRITE 



SIGMAl : 

RESET INTERRUPTS (HALT & EOC) 
READ NERDAS, SET ALARMS AND 
OVER-RIDE VALUES 
SET AND START INTEGRATORS 
CALCULATE: 


0, _ (4tt)^ ^wi 

M A7- 


i = 1 , 8 


SIGMA2: 

READ NERDAS DATA, 

READ FILTERS, F . , i = 1, 12 
COMPUTE: ^ 

a”i = a°*. PR^i/GGLi, or 

a°i =0°'^ PR^^/GG^^, i = 1,8 

FINISH LOADING BUFFER LINE 
WRITE: 

OUTPUT A LINE TO Bi-tj, L 


FIGURE 8,8 4.75 GHz SCATTEROMETER SYSTEM SOFTWARE 




zero values will be produced at the data output port -- thus providing 
a continuous product for the user. 

The cycle time, or time delay between each record at the output, 
will be a function of the microprocessor overhead time and the integra- 
tion time. The integration time is radar frequency and aircraft ground 
velocity related and thus is somewhat independent of any software 
optimization scheme. However, the microprocessor duty cycle can be 
adjusted within certain limits to minimize the overhead. This should 
be done during program development testing by properly balancing tasks 
between the two sigma-zero calculation sections, SIGMA! and SIGMA2, of 
the RUN module. Figure 8.8. 

8.2.4 "RESET" Module A'l{(orithm 

This module provides all program Interrupt Service, as illustrated 
in Figure 8.7. From the operator's point of view, it is the module through 
which all communications are made with the software. In addition to 
operator inputs, the module handles all "pre-RUN" housekeeping functions 
such as Interrupt Structure set-up, Input/Output port initializations, 

Data Buffer initializations and resets all flags, alarms and the Front 
Panel . 

The RESET algorithm can be stated as below where the notation 
(ABC) 0 means that the value zero is assigned to the location ABC: 

1. Disable Interrupts 

2. Initialize the Interrupt Controller (8259) to accept 
all but Interrupt-1, initialize the CZT-Board with 
DMA address, and set-up the Bi-cp L board. 

3. Reset Data Buffers, and their Flags, and Pointers 
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4. Reset Front Panel 

5. Reset System Alarms and Flags 

6. Reset Interrupt Processor 

a. ^(STACK) ^ 0 

b. ‘^(DECFG) ^ 0 

c. ^(DISFLG) ^ 0 

d. ^(ENTER) 4- 0 

e. ^(X) 0 

f. ‘^(NNBR) ^ 1 

7. Enable Interrupts 

8. Wait for Interrupt i, i = 0, 2, 3, ...» 7 

a. If Received, transfer control to the 
appropriate processor, otherwise go 
to 8. 

The above algorithm is illustrated in flow chart form in Appendix F, 
Chart 1.0. Each step is represented by a box in the chart. Each box has 
a number whose second digit corresponds to the above step numbers. 

While operating in this module the Interrupt-1 (i-1). Figure 8.7, 
will not be accepted, i.e., it is masked out during the set-up. Step 2, 
above. Although the Interrupt-3 is acceptable, no purpose is served by 
invoking it in this mode of the operation of the system. 

The Interrupts-0, -2, -3,..., -7 are each created by the operator when 
a particular key on the keyboard is depressed. Table 8.1 summarizes the 
key types, names, and their respective assignments within the structure. 
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TABLE 8.1 INTERRUPT/KEY ASSIGNMENTS 


KEY TYPE 

KEY NAME (LABEL) 

INTERRUPT CREATED 

DIGIT 

0, 1, 2,. ...9, •, CHS, 
VEL, ALT, DRFT, ROLL, 
PITCH 

INTERRUPT -0 

STATE 

RUN 

INTERRUPT -2 


HALT 

INTERRUPT -3 


RESET 

INTERRUPT -6 

FUNCTION 

DISPLAY 

INTERRUPT -4 


ENTER 

INTERRUPT -5 


CLEAR 

INTERRUPT -7 


Whenever the opigrator, while in the RESET mode of operation, depresses 
a key, program control is transferred directly to one of the seven 
interrupt processing modules. Figure 8.7, tagged with i = 0, 2, 3,..., 7. 

Each of these modules has its own special purpose algorithm. Starting 
with i = 0, each will be outlined in enough detail to illustrate its 
principal functions and create an understanding of the Flow Charts in 
Appendix F. 

The RDKEY module handles all digit key activities. These activities 
are all related to entering the appropriate aircraft data over-ride values. 
Again, the algorithm is activated when a digit key. Table 8.1, is depressed. 
The RDKEY algorithm may be stated as follows (the notation ^(ABC) means the 
content of Location ABC): 

1. Read key value into a "save" location, INDIG. 

2. Interpret its value and 

a. If the value ^ II^q go to DECIM, Step 7. 

b. If the value IO-iq go to CHS, Step 9. 

c. If the 'ENTER' key has been depressed go 

to FARM, Step 11. 

d. If this is the 6th digit go to ERRl , Step 14. 

3. Save digit character and echo character to the 

Front Panel . 

4. Add 1 to the digit counter. 

5. Send End-of-Interrupt (EOI) 

6. Return to WAIT LOOP. 

7. If DECFG 0 go to RETURN, Step 5. 

8. Set DECFG, go to Step 3. 
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9. Multiply c(nnbr) by -1. 

a. If ‘^(NNBR)>0 ECHO ' + ' to Front Panel. 

b. If ‘^(NNBR)<0 ECHO '-'to Front Panel. 

10. Go to Return, Step 5. 


11. If the 'ENTER' key has not been depressed go to 
DISPL, Step 15. 

12. If the digit, ^(INDIG), is £ 3 go ERRl , Step 14. 

13. If ^(INDIG) < 8 go to ACPRM, Step 20. 

14. ERRl, (entry procedure error), 

a. Disable all interrupts except 'CLEAR' or 'RESET'. 

b. Flash Front Panel at 500 msec intervals. 

15. If ‘^(DISFLG) = 0 go to ERRl, Step 14. 

16. If '^(INDIG) > 9 go to ERRl, Step 14. 

17. Display parameter name and value on the Front Panel. 


18. Reset stack and flags: 

a. 

c (STACK) ^ 0 

b. 

C(ENTER) ^ 0 

c. 

c(DECFG) ^ 0 

d. 

c(DISFLG) ^ 0 

e. 

c(X) ^ 0 

f. 

c(NNBR) ^ 1 

19. Go Step 5. 


20. Save parameter string and sign 

21. Save parameter value 

22. Set bit in Over-ride Flag word 

23. Display '*' in appropriate LED on the Front Panel. 

24. Go to Step 18. 
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The above algorithm is illustrated in flow chart form in Appendix F, 
charts labeled Interrupt-0. 

Moving up the stack of modules in Figure 8.7, the CLEAR, Interrupt-7, 
processor has a relatively simple algorithm: 

1. Blank the digit LED's on the Front Panel. 

2. Reset the Interrupt Processor (Go to 
Step 6 of the RESET module). 

Next up the stack is the Interrupt-6, created by pressing the RESET 
key. This function does the same thing from the software point of view 
as turning power off and on again. It simply transfers program control 
back to Step 1 of the RESET algorithm. 

Interrupt-5, handled by the ENTR module is also relatively simple. 

Its function is to signal the software that the next key depressed will 
designate the destination of the parameter already keyed into the panel. 

It may be stated as follows: 

1. ‘^(ENTER) ^ 1 

2. Flash Front Panel LED's one time to signal acceptance. 

3. Go to Step 5 of RDKEY module, (Return to Wait Loop). 

The Interrupt-4 function interacts with the RDKEY module. By 
depressing the DISPL key the software is signaled that the next digit 
key pressed will designate the contents to be displayed. The algorithm 
follows: 

1. ^(DISFLG) ^ 1 

2. ^(ENTER) ^ 2 

3. Go to Step 2 of ENTR module, (flashes Front Panel 
and returns to Wait Loop). 
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Now, the operator, by pressing the appropriate digit key next will cause 
the control to pass through the RDKEY module to Step 18. 

The HALT module processes Interrupt-3, (i = 3, Figure 8.7). Its 
function is to perform an orderly close out of the data production 
operation, SIGMAl and SIGMA2 modules. It does the following: 

1. Disable Interrupts and reset mask for Interrupt-6 only. 

2. Build a 'Fill' frame for End-of-File identification. 

3. Write 50 'Fill' frames into the output buffer. 

4. Reset Front Panel 

5. Enables Interrupts 

6. Wait for Interrupt-6. 

The remaining two Interrupts in Figure 8.7 are associated with the 
'RUN' mode of operation. When the function key 'RUN' is depressed, 
control transfers directly to the SIGMAl module. Figure 8.8. While in 
this mode of operation the integration time on the CZT-board, will 
provide the Interrupt-1, (i = 1, each time it "times-out", or completes 
the cycle of integration provided to it by the SIGMAl module. 

8.2.5 "RUN" Module Algorithm 

As discussed in Section 8.1, the objective of the software is to 
repeatedly solve the estimated scattering coefficient equation 

:o _ 1%)! ^ ’’'’i 

1 '^s 

This is done for each of eight pre-selected viewing angles with either 
like or cross polarized input data. These viewing angles, 6|^^. , i=i, 8, 


are pre-loaded into the software object code and are not changed except 
by reprogratwning the entire software package. 

For each operating cycle (an output data set) the software will 
analyze the input data calibration tones to determine if horizontal (H-) 
or vertical (V-) polarized data is being transmitted. A front panel 
input will be sensed by the software to determine whether the receiver 
is in Horizontal or Vertical polarization mode. 

Having resolved these variables, the software system will calculate 
in sets of eight one of two forms of Equation (8.1); i.e., 


Like-Polarized: 


PR,/GGL^|_i 


(8.26) 


Cross-Polarized: 


“°xi = “°'i 


(8.26a) 


where the parameters are 


PR.., received/transmit power ratio in band i; 

p, 

K. 


cr .: , 

A 


for band i ; 


GG . , Antenna gain for viewing angle i, cross- 

A I 

polarized gain table. 

Aj^.j, ground cell area for viewing i using 
like-polarized beamwidth tables. 

» ground cell area for viewing angle i using cross- 
polarized beamwidth tables. 
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5/hen indexing into the beainwidth find gain tables the software will 
use the polarization tone and front panel input to determine whether 
to use the horizontal or vertical tables and whether to use the Like or 
Cross section in the table. Note that if a constant beamwidth value 
is used for all polarization combinations, the computation of will 
be done as part of the parameter. 

The system of calculation defined for a estimation in modules 
SIGMAl and SIGMA2 presume that the hardware is structured for single 
channel CZT operation, Figure 8.9; however, with appropriate changes 
in the set-up, control, and output portions of the SIGMAl and SIGMA2 
algorithms the software package could be used with a two-channel hard- 
ware system, Figure 8.10. The implications of using a system like that 
in Figure 8. ?0 are: 1) the processor must handle twice the data, either 

serially or parallel, thus the ground cell to ground cell separation 
will increase significantly, 2) the memory buffer size needed doubles 
for both the CZT and the microprocessor (yP), 3) the input/output 
times are doubled, and 4) the software must have a means of knowing 
which channel has the Like-Polarized data and which has the Cross- 
Polarized, (this could be handled by always assigning Cross-Polarized 
data to CZT #2 as suggested by Figure 8.10). 

8.2. 5.1 SIGMAl Sub-Module Algorithm, INTERRUPT-2 PROCESSING 

In the most basic form the algorithm for the SIGMAl module which 
calculates the a'° estimations may be stated as follows: 
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SIGMA! : 


1.1 Reset Interrupt mask to acknowledge the EOC and 'HALT' 
interrupts only. Read aircraft NERDAS and establish 
velocity, V; altitude, H; pitch, p; roll, drift, (}>; 
and clock time, T^, 

1.2 Evaluate L^, i=l, 8 using Equation 8.24. 

1.3 Compute Integration time using L-j from Step 2, 

Equation 8.3; start Integration on the CZT-board. 

1.4 Compute Output Buffer length, NL. 

1.5 Set values of load pointers and write pointers, up-date 

line pointers in the buffer. Move A/C data into 
buffer. 

1.6 For i=l, 8 using 0|^^. , do: 

a. Compute 0'^., X., Y., 0p. 

b. Evaluate fdc^. , using Equation 8.5 

c. Evaluate Bd^. , Equation 8.4 

d. Compute NF^. , Equation 8.13 

e. Calculate Ba^., Equation 8.14 

f. Evaluate fac^. , 0jQ.j j » Equations 8.15 

through 8.19. 

g. Using Ba^. and NF^. calculate index pointer, PRL^ 
(left most index), into the Power Spectral Density 
(PSD) data set, Equations 8.14a through 8.14d. 

h. Calibration and Polarization index pointers are 
stored constants 

i. Look up roll-off value, Z(fac^. ) 

j. If beamwidth is constant calculate area A^. . 

k. Evaluate the term , Figure 8.8. 

1.7 Wait for Interrupt-1, (i=l. Figure 8.5). 
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8. 2. 5. 2 SI6MA2 Sub-Module Algorithm, INTERRUPT-1 PROCESSING 

Upon receipt of Interrupt-1 program control is passed directly to 
the SIGMA2 module which does the following: 

SIGMA2 (AND WRITE ) : 

1.1 Transfer CZT data to yP memory, after interrupt 
Flag checks and resets 

1.2 Using pointers from Step 6.h, evaluate calibration and 
polarization power. Set polarization Flag Bits and 
calibration alarm (if appropriate). 

1.3 For i=l , 8 do: 

a. Evaluate PR^. , power received in band i 

b. Look up G Gx, gain for antenna, (and if appropriate 

I w 

BMW^. , beam-width for Ql-,-)" 

c. Compute A^. , area for ground cell i, unless already 
done in Step 6.k of SIGMAl. 

^0 ^0 

d. Evaluate o or a equation 8.26a or 8.26b. 

1.4 Store viewing angle values, a or a values. Alarms 
and Flags in the output buffer, formatted as shown in the 
Flow Charts. 

1.5 Pack and write the appropriate output buffer line to the 
Bi-(J) L interface. 

1.6 Set the Pass Number, ^(NPASS) 1 

1.7 Read aircraft NERDAS and establish V,H,P,i|;,(i), and Tj^, 
the current clock time. 

1.8 Using previous time, T^, (SIGMAl, Step 1) compute last 
cycle time, T^. 

1.9 Move Tj^ to T^. 

1.10 Go to SIGMAl, Step 2, after clearing the Time Alarm. 


The system will remain in the computational cycle until the 
operator depresses the 'HALT' key. 

8. 2. 5. 3 Additional Considerations for SIGMA! and SIGMA2 

In conjunction with preparing the software code which will execute 
the SIGMAl and SIGMA2 algorithms there are additional considerations and 
details related to some of the above steps. These are summarized as 
fol 1 ows : 

SIGMAl : (Additional considerations): 

Step 1 - The code for this step must include logic for 
a) selecting over-ride values in place of NERDAS data 
received, b) evaluating the values of and cj) 

against some standard (if not over-ridden) and setting 
the computational value to the input or to the standard, 
c) setting a bit in the output data frame to indicate the 
type data used, and d) writing an '*' to the Front Panel 
LED when failure occurs. NERDAS failure must be defined 
prior to coding this logic. Buffering data into Output 
Data Line cannot be done here. Where appropriate, over- 
ride or default values should be substituted into the line. 

Step 2 - In calculating L. , i=l , 8, the expression in 
Equation 8.24, L- = (hA0')/cos 0, . , the parameter A0' is 

p AQ 

a stored constant whose value is calculated by A0/cos (-^) 
where A0 = 3° (0.05236 radians), so that A0' = 0.05239 
radians . 

Step 3 - Integration time for the kth cycle is calculated 
using Equation 8.3 where L-j is the length of the ground 
cell associated with viewing angle 0 ^ 1 ; i.e., the smallest 
0j^^. . Integration is started by sending the CZT Board the 
number of sub-records desired, NS, where 
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(8.27) 


NS = [(Tj/512)fg + 0.5]j , 
or NS = [Tj/Tg + 0.5]j 

where = 51 2/ f^, and as before = sample frequency. 

and f„ will be stored constants, 
s s 

Step 4 - Output buffer length is computed by 

NL = [(H taneLg)/VT^ + 0.5]j (8.28) 

where 0j^g is the largest viewing angle and is the cycle 
time. When ^(NPASS) = 0, the first pass through SIGMAl , 
will not have been evaluated, therefore the software 
should use a dummy value, = 2Tj. The value of the tan 0j^g 
should be a stored constant. If the value of NL is greater 
than the reserved buffer length (200 lines) the software 
must set an alarm; i.e., byte 79H (xx H means Hexadecimal 
notation) in the output data frame, and the Buffer-Over- 
Flow LED on the Front Panel gets and written into it. 
Finally, provisions must be made to set NL to the maximum 
buffer length value when the Over-flow occurs. 

Step 5 - On the first pass through SIGMAl, ^(NPASS) = 0, 
the buffer bias value is set to NL + 10; i.e., 

MXo = NLq + 10 (8.29) 

represents the line number at which the output buffer 
loading begins. On each successive pass the value of MX 
is incremented so that in general 

MXk = MXjj + k (8.30) 
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where k = 0,1,2,..., 255 and NL^ is the initial value of 
NL. The value of k is simply the pass number. Care must 
be taken when the register accumulating MX over-flows, 
else a frame of data will be lost. This is averted by 
simply setting MX = 1 wherever over flow occurs. The 
input or load-line pointer is calculated by 

= [MXi^/Mo]*^ (8.31) 

D 

where [ ] signifies the remainder of the integer division, 
and Mo is the modulus for the buffer. Its value is 

Mo = NLq + 20 (8.32) 

so that the value of LLj^ cycles the load point through the 
buffer. Line pointers, LP^ , for saving each of the sigma- 
zero values and their associated viewing angle is evaluated 
by 


'-'’l = (8.33) 

Lo 

If LP. < 0, set LP. = NL + LP.. 

Column pointers into the buffer for each value will be 
determined by their respective byte position assigned 
in the Bi-c|) L output line, see Flow Charts for the 
•WRITE' module. 

The index of the output line which is sent to the Bi-<j) L 
interface is calculated by 

KO,^ = [(MX^ - NLq)/Mo]'^ (8.34) 

If KO < 0, set KO = KO + NL^. 
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step 6 - This step must first take care of those roll angles, 
i|;, which exceed the desired viewing angle . Thus, the 
appropriate order of computation is as follows: 

a. 1) if !♦! > 6^, = -I'd 

2 2 

2) tan 0'^^^ = -cosijj (tan 0j^^ - tan tp) ‘ 

3) XT = (-H/cosi|)) tan0'|_^- 

4) YT = H tanijj 

5) X. = XT cos<p + YT sin4>, but if X. < L./Z, set X. = L./2 

6) Y. = XT C 0 S( 1 ) - YT cos<p 

7) XD = + Y.^ 

8) 00^ = Tan-'' (X.^/XD)^^ 

b. Evaluate 6.b through 6.g as indicated. 

c. Calibration and polafi2;ation index pointers are evaluated 
before program loading and are stored constants. Consid- 
eration must be given, however, to the type of scattero- 
meter when these constants are calculated. 

For C-band scatterometers the left most pointers into 
the PSD (Power Spectral Data) for calibration tone power 
are found for horizontal and vertical polarization by 

PCH^ = [2860 (^) + 0.5]j (8.35) 

PCVl = [3280 (^) + 0.5] j (8.36) 

when the tones are inserted on the lower sideband. 

However, if the insertion is on the upper sideband 
then 


(8.37) 


PCHj^ = 512 - PCHl 
and 

PCV^j = 512 - PCVl (8.38) 

Note that all pointer expressions assume that the filter 
spectral lines have been ordered with respect to zero 
frequency; i.e., = 0, w, = f^/512, = 2f /512,..., 

for i £ 255. However, zero-frequency (dc) 
actual!^ appears in the CZT spectrum at wg* thus all 
pointers must be biased by +2. Also, this assumes that 
the discrete Fourier transform of X + jY has been 
implemented by the CZT board where X is the true cosine 
channel and Y is the true sine channel. If significant 
power is present at either PCH or PCV with respect to 
the noise level then either horizontal polarization or 
vertical polarization, respectively has been transmitted. 

For L-band scatterometers the calibration tone appears 
in double sideband form on the cosine channel at 1900 Hz. 

The polarization tone is inserted similarly at 2100 Hz. 

The pointers are 

PCC = [1900 (^) = 0.5]j (8.39) 

s 

and 

PCV = [2100 (^) + 0.5]j (8.40) 

for calibration and polarization respectively. When no 
tone at PCV appears (compared to noise level) then hori- 
zontal polarization has been transmitted. 


Noise floor pointers may be calculated for C-band or 
L-band respectively by 

PNZ = [1800 (^) + 0.5], (8.41) 

C % ^ 

PNZ, = [1200 (^) + 0.5] j (8.42) 

Step 7 - It is quite possible that Interrupt-! will occur 
prior to completing Step 6.1 of SIGMA! . Therefore, the 
module should be coded to handle such an event. The flow 
charts in Appendix F indicate a way to process the interrupt 
and alert the operator to that f^ct. This capability is 
designed to be used in the program development and testing 
only. 

SIGMA2 (Additional considerations): 

Step 2 - This step should accomplish the following; 

a. Evaluate power at the pointer PCH and PCV if C-band 
(or PCC and PCV if L-band). 

b. Evaluate power at the noise floor, PNZ, Equation (8.41) 
or (8.4^. Reset calibration alarm. 

c. For C-band scatterometers: 

1) If PT|^ (power at PCH pointer) » PT^ (power at PCV 
pointer) set ^(PZFLG) = 2 to indicate horizontal 
polarized transmitter. 

2) If PT^«PTy set ^(PZFLG) = 0 to indicate vertical 
polarized transmitter. 

3) If ‘^(PZFLG) = 2 compare PT^ to PT,^ (power at PNZ 

pointer) and if 9° to C.(5), else set 

calibration alarm and go to C.(5). 

4) If ^(PZFLG) = 0 compare PT^ to PTj^ and if PT^»PTj^ 
go to C.(5), else set calibration alarm. 

5) Continue, transmitter type and quality have been found. 


d. For L-band scatterometers : 

1) If PT (power at PCC pointer)»PTf. go to d. 

(2), else set calibration alarm. 

2) If PTy (power at PCV pointer)»PT|^ set ^(PZFLG) 

= 2 to indicate horizontal polarization and go 
to d (3), else set ’^(PZFLG) = 0 to indicate 
vertical polarization. 

3) Continue, transmitter type and quality have been 
found. 

Step 2 - Total received power within each band associated 
with is calculated by 

PRL.+N.-l 

PR. = E DSL. (8.43) 

k = PRL. 

where DSL|^ is the magnitude squared of kth spectral line coefficient of 
the 512 point DFT. In the CZT implementation of the DFT the Fourier 
components are squared and integrated over sub-records in a dedi- 
cated hardware processor, so that only the summation over the 

NF. = PRR. - PRL. + 1 (8.44) 

set of spectral lines is required to estimate the returned power. 

Actual transmitted power is found by correcting for the fact that 
the calibration tone may not appear in the center of the filter band- 
width. The measured transmitted power is adjusted by the shape function 
in Equation (8.21). 


When using g?in and beamwidth look-up tables the values of 9'^^^ 
will be used to index into the tables, rather than the "desired" view- 
ing angles, 6^^. . 

Step 4 - All data except flag bytes going into the output 
buffer lines must first be converted to binary coded 
decimal (BCD). 

Appendix F shows the line format. 

Step 5 - As each output buffer line, indexed by KO, Equation 
(8.34), is written to the Bi-(J) L port it should be backfilled 
with BCD '9' in each byte. 

8.2.6 Output to the Bi-<f) L Interface 

Output to the Bi-(J) L interface is accomplished by transferring the 
contents of the block of memory, one byte at a time, to the Bi-t() L output 
port. The software used in the 13.3 GHz processor built at TAMU in 1977 
may be used with appropriate modifications to account for the differences 
in the BCD data frame. 

8.2.7 Estimated Core Memory Usage 

Expected core usage may be divided into five principal blocks 
as illustrated in Figure 8.11. This breakdown does not imply specific 
address spaces for the modules involved. The size of each usage space 
has been conservatively estimated and is expected to be somewhat less, 
but of course will only be determined when actual coding is accomplished. 

Block 1 is a IK segment set aside to handle all hardware related 
address ports; e.g.. Front Panel, Bi-(|) L, NERDAS, and 8259 (Interrupt 
Processor). 


BLOCK 


0 

IK 


9K 


UK 


25K 


TOTAL: 64K 


1 

HARDWARE RESERVE AREA (RAM) 

2 

CONTINGENCY ROM/RAM 
AND 

CZT BUFFER AREA 

3 

WORKING STORAGE AREA, 
< 2K bytes (RAM) 

4 

RESET, RUN, AND HALT 
Program Modules 
< 14K bytes (ROM) 

5 

Time Aligned Output 
BCD Data 


^ 200 lines x 200 bytes/ line 
(RAM) 


FIGURE 8.n ESTIMATE CORE USAGE 



Block 2 is contingency space containing both RAM and ROM, up to 
8K bytes. 

Block 3 is a set of RAM reserved for working storage; i.e., space 
required for data manipulation (sigma-zero calculations), NERDAS decod- 
ing, Bi-(|) L encoding, input CZT data, etc. 

Block 4 is a 14K byte set of ROM reserved for the main driver 
program; i.e., RESET, RUN, and HALT module. 

Block 5 is a 39K byte region of RAM required for the BCD Output 
data set containing sigma-zero values aligned with NADIR time tags. 

This buffer is large enough to allow aircraft operations up to altitudes 
of 6500 feet above the terrain being viewed. 

8.3 Operating Procedures 

8.3.1 Set-Up 

Consistent with the design goals of the system, minimum set-up 
operations will be required. Operator interactions consist of only 
two functions: 

a. Power - Turn-On: assure that the system has proper electrical 

power and that In-Put/Output ports are properly connected and 
powered . 

b. Enter any NERDAS over-ride parameters. 

Over-ride parameter inputs, as discussed in Section 8.2.2, must be entered 
before placing the system in the "RUN" state. 

8.3.2 Operator Inputs 

The operator may communicate with the system through the keyboard 
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only. The keyboard will have three state keys (RESET, RUN, HALT), 
eleven digit keys (0,1,..., 9'.' CHS) and three function keys (CLR, 

DISP, ENTER). The state keys are used in place the system in the 
desired mode (or state) of operation. When power turn-on is accom- 
plished the system is automatically placed in the RESET state. While 
in this state the operator should enter those aircraft data values, if 
any, that will be used in place of NERDAS values. After the over-rides 
are entered, they may be verified by the display function (DISP key). 

If during entry of value on the keyboard a key error (0,1 ,. . . ,9, ' . ' , 

CHS) is made, that entry may be cleared from the board by pressing the 
CLR key. The appropriate sign may be given to the value by pressing the 
CHS key. Once the value desired is placed in the keyboard display, the 
operator may implement the over-ride by pressing ENTER and then the 
destination key; i.e., VEL, ALT, DRIFT, ROLL or PITCH. Procedural errors 
during over-ride entries and during display requests are signaled by 
flashing the front panel display. Return to the RESET mode is then 
made by pressing the CLR key. Table 8.2 summarizes the operator input 
procedures in their correct sequence. 

Only the one error message (flashing display) is provided to the 
operator while in the RESET state. This always signifies an invalid 
key sequence during data entry on display functions. 

System defaults are needed for those values of the aircraft data 
set which have not been over-ridden, but fail to pass a reasonable value 
test in the NERDX module. It is suggested that these default values be 
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TABLE 8.2 OPERATOR INPUT SEQUENCE 


Step 1. Power— on and I/o connections made. 

Step 2. Over-ride Values: 

a. If no over-rides desired go to Step 3. 
u. Key in value and sign, (floating point number) 
Example: (key order is left to right) 

"-0.19", displayed shows - 0.19 

press "ENTER", display blinks 
press "ROLL", display shows 

R 0 L = - 0.19 

c. Press CLR to blank front panel display. 

d. Key in next value; as in 6 above, until all desired 
over-rides entered. 

e. Sequentially display all values previously entered. 

Example: press "DISP", display blinks 

press "ROLL", display shows 

R 0 L = - 0.19 

Step 3. Change to "RUN" state by depressing RUN key. 

Step 4. Change to "HALT" state by depressing HALT key. 

Step 5. Change to "Reset" state by depressing RESET key. 


selected based upon the most likely operating envelope for the system; 
e.g., ALT=1500', VEL=77.1 mps, PITCH=R0LL=DRIFT=0.0. 

8.3.3 System Outputs 

8. 3. 3.1 The Output Data Frame 

Each cycle through the SIGMA1/SI6MA2 modules will cause an output buf- 
fer line, or 'frame', to be written to the Bi-(j) L interface. Each 
frame will contain up to 200 bytes of BCD data and flags. As the frame is 
now designed, pages 47-50 of Appendix F, it will consist of 128 (80H) bytes. 
All data will be coded into BCD except the flag words 78H through 7FH. 

8. 3. 3. 2 Data Frame Annotation 

Output flags are place in the frame to provide the user with suf- 
ficient data to analyze sigma- zero results. These are information bits 
which indicate specifically the conditions and over-rides, if any, used 
in calculating the sigma-zero estimates. 

Each output frame will contain the following trailer bytes ('>' 
indicates binary notation). 

Polarization: >11 = HH,>10 = HV,>00 = VV,>01 = VH 

H/V Alarm: >00 = off, >01 = on 

Time Alarm: >00 = off, >01 = on 

Calib. Time Alarm: >00 = off, >01 = on 

Over-ride Flag: OOH = No over-rides used 

Bit 6 on = alt. over-ride used 
Bit 5 on = vel . over-ride used 
Bit 4 on = pitch over-ride used 
Bit 3 on = roll over-ride used 
Bit 2 on = drift over-ride used 

Std-value Flag: (same format as over-ride) 
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9.0 CONCLUSIONS AND RECOMMENDATIONS 


9.1 Executive View 

The CZT approach represents a novel and powerful technique for 
filtering the fan beam scatterometer return into many parallel Doppler 
channels. When the radar return is regarded as a complex valued signal, 
the CZT, when appropriately modified, is capable of filtering the fore 
and aft spectra simultaneously. When amalgamated with a dedicated 
digital post processor, the PSD of the radar return can be estimated 
with precision over a large dynamic range as required by the character 
of the radar return. Since many parallel spectral estimates are avail- 
able simultaneously the filtered data can be selected adaptively with 
changes in aircraft velocity to track the desired viewing angles. 

Since the selection of spectral data is under software control, 
the mode of operation can also be altered to satify the needs of 
specific experiments. As a consequence, the spectral data could be 
selected to provide constant precision estimates, constant ground 
resolution estimates, constant angular resolution estimates or combi- 
nations of these modes. The operating system developed within this 
effort is based on a constant angular resolution mode. An angular 
resolution of three degrees is provided for the C band system and six 
degrees for the L band system. Higher or lower angular resolutions 
are obtainable for both scatterometers. 

The cha i coupled device and associated circuitry were demonstrat- 
ed as a good approach to filtering the radar data. A dynamic range of 
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more than 50 dB was observed in the evaluation tests. Excellent linearity 
was demonstrated over the 50 dB of dynamic range. The quadrature process- 
ing technique yielded a fore to aft signal rejection ratio of more than 
30 dB. 

A comprehensive software design was developed to control the process- 
ing system and to convert the spectral estimates into calibrated estimates 
of a°. The software was designed into flexible modular units having cohe- 
sive objectives. Flags were included in the design to assist in testing^ 
the program once the coding is developed. An identical software 
architecture was utilized for the C-band and L-band processors. The same 
architecture is extendable to other NASA fan beam scatterometers. 

Various interface units were also designed to provide communication 
in and out of the system. Among these, an operator console and display 
unit was designed to permit operator control of the system. The console 
features a keyboard entry and an alpha-numeric display. The display 
incorporates many intelligent functions such as clear display, shift left 
or right, blink, read/write, etc. It consequently represents a powerful 
technique to monitor the operation of the processing system when used 
appropriately. For example, it will be helpful in debugging the soft- 
ware during the future development stages. 

An aircraft interface unit was also designed to transform NERDAS 
data to a binary form for use by the scatterometer processing system. 

A similar interface unit was actually fabricated and proven under 
another NASA contract [1] and features a level control to condition 
the Bi-(j) L modulation for decoding and a special digital decoder 






174 


to properly identify the binary transitions. 

The afore described processing system consequently meets or exceeds 
the requirements for a flexible operator and experimentor oriented 
system. The system architecture is sufficiently universal to handle the 
data characteristics unique to each of the NASA scatterometers. It is 
conceivable that a single processor could adapt to each of the scattero- 
meters, if desired. However, it is recommended that a processor be 
developed for each of the scatterometers to provide multi -sensor data 
simultaneously. In the latter case, since a common design approach 
applies to all the systems, the development and fabrication costs will 
be considerably cheaper. 

It is important to point out that the efforts reported here repre- 
sent a comprehensive paper design with certain elements of the design 
partially proven through breadboard testing. This work is consequently 
a sizable first step in a design and development program. However, a 
major effort is yet required to fully develop the design since a design 
and development program is an iterative process including design, fabri- 
cation and evaluation. 

The results of the current effort have demonstrated that it is 
necessary to develop and evaluate an engineering model of the processor 
system before fabricating the target system (see Section 6.0). An 
evaluation of the engineering model is required to assess the actual 
speed of the processor and the interaction among all elements of the 
system. In approaching the development of the engineering model, it 
is recommended that certain key problem areas, identified within this 
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effort, be investigated before implementing the developmental system. 

First, it is necessary to perform a comprehensive study to verify 
the current a® inversion technique or to replace it with a better one. 

The current inversion technique is based on a intuitive interpretation 
of the scatterometer equation and has not been verified through simula- 
tions. It is essential that this be performed to assure accurate 
estimates accounting for the effects of the antenna pattern, Doppler 
filter and aircraft attitude. Secondly, the software design should be 
encoded to verify timing, optimize interrupts and develop certain inter- 
active coding for debugging purposes' and system check-out. Lastly the 
breadboarded PSD estimation subsystem should be employed as a "test 
bed" to aid in the development if the engineering model. In this regard 
the source of the noise floor problem described in Section 6.5.3 should 
be identified and remedied. Then the software and hardware should be 
united to demonstrate system capability. 

9-2 Engineering View 

Experience in debugging and evaluating the PSD estimation breadboards 
has led to a number of insights which could lead to improvements or 
simplifications in the system. They are simply listed below: 

1) Consider digitization to 8 bits rather than 12. The additional 
resolution may be available in averaging records and spectral 
lines. The accumulator word size could conceivably be reduced 
to 24 bits. 

2) Consider reducing memory size to store only aft data. 
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3) Use doubly clad PC board to improve noise immunity of analog 
and A/D electronics. 

4) Consider methods of shielding the filtering unit from the 
digital and controller boards. 

5) Replace discrete sample and hold circuits with integrated 
circuits. 

6) Establish the methodology to balance the common mode 
rejection in the differential current integrators. 

7) Evaluate the effect of op amp offset voltages on pre-chirp 
multiplication accuracy. 

8) Reduce the bandwidth of the input amplifiers by a factor of 
two and compensate for the loss of calibration tow amplitude 
in the software. 

9) Optimize CZT gating signals to reduce probable noise sources. 
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APPENDIX .A 

OPTIMIZING THE PRECISION 
FOR A GIVEN COVERAGE 

Theorem; 

It is shown here that for a given coverage interval 
Lc = Pb + vT 

that the precision factor BT can be maximized if the scan time is 
chosen so that 

Pb = vT 

ie., the beam resolution is equal to the scan length. As a consequence, 
BT can be maximized if the beam resolution occupies half the coverage 
length during a single integration period. The reader will note that 
the result is independent of view angle. 

Proof ; 

The Doppler frequency on the ground track is given by 

X _ 2v Sine 

A . 

When the x coordinate is chosen in the up-track direction as illustrated 
in Figure A-1, the Doppler gradient along the x coordinate is given by 


where x = h tane has been used. For a fixed bandwidth the precision 
factor may be related to the gradient in the following way 


HT _ df ^'“c " 

- 3x Pb V 


The factor (L^. - Pq)/v is clearly the integration time. Now the 
precision can be maximized with respect to beam resolution by 
determining when 



This yields the result 
_ *"c 

Pb ■ T 


which is equivalent to 
Pb = vT 

It is easily shown that this criterion results in a maximum BT rather 
than a minimum. 
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APPENDIX B 


THE THEORY AND APPLICATION OF 
THE ALIGNMENT GENERATOR 

The alignment generator produces quadrature periodic square wave 


-T/g < t < -T/^ 

-T/i < t < T/^ 

T/4 < t < T/2 


0 < t < T/g 

T/g < t < T 

over a single period. As can be easily shown, the spectra associated 

with these periodic signals are given by 
00 

Sj = E i. ^ ~ ^ n) 

n = -00 

n odd 


signals given by 
-1 

Vj = +1 

-1 

and by 




1 

-1 


and 


Op 

1 j TTn /) 

Sn = ^ 1 6 e^T 6(0) - ^ n) 

y “ n y 

n = -«> 
n odd 
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t 






respectively, where 1/n is negative on the set ne -7, -3, +3, 

+7,...} and positive on ne -5, -1, +1, +5, Sj and Sq are 

illustrated in Figure B-1. If Vj is inserted into the I channel and 
Vq is inserted in the Q channel , then the CZT processor will yield 
a spectrum given by 

So ' 5, + J Sp 

The spectrum of is shown in Figure B-2. It is noted there that the 
spectrum is not symmetric and that alternate lines have disappeared. 

The complementary spectrum can be formed by interchanging the inputs to 
the channels. Then the complex valued linear combination of spectra 
becomes 

So' ■ Sq * j s, 

which is illustrated in Figure B-3. It is noted that 

Sq' (w) = j Sq(-o)) 

i.e., Sq'(w) is the convolved and phase rotated spectrum of Sq(w). 

The spectra of Sj, Sq, Sq, and Sq' are helpful aligning the CZT 
filtering unit. The gains through the four chirping channels must be 
matched to assure operation of the CZT transversal filter consistent 
with the theory cited in Section 3.4. Before using the alignment gener- 
ator, the dc offsets in the channels prior to the chirping multipliers 
must be removed. An offset adjustment is provided in each of the four 
channels. 
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Ga in Adjustments 

Inject Vj into the I channel (x) with no input on the Q channel (y). 
Adjust and Xg channel gains so that Re and Im inputs to the trans- 
versal filter have identical peak to peak amplitudes. An oscilloscope 
may be employed to observe the Re and Im inputs to the transversal filter. 
Also adjust gains in the y^ and ^ 2 , channels in the same way by injecting 
Vj into the Q channel. 

After balancing the gains it is advisable to fine adjust the dc 
offsets in all four channels. 

DC Pedestal Adjust 

The chirped signals must ride on a positive bias when entering the 
CCD transversal filter. A bias somewhere between 4 and 9 volts is 
required to achieve the full dynamic range inherent to the transversal 
filter. 

To adjust the dc pedestal inject Vj into the I or Q channel and 
monitor the delayed output VOR with an oscilloscope. Adjust the dc 
pedestal on the Im channel so that the delayed output is neither com- 
pressed or clipped at the top or bottom of the chirped waveform. For 
maximum undistorted dynamic range, the input to the transversal filter 
should be near 1.5 volts peak to peak. 

Similarly, adjust the dc pedestal on the Re input by monitoring 
the delayed output VOI. 

The envelopes of VOR and VOI will be amplitude modulated over the 
CZT cycle. This is normal for the CCD device and is induced by a high 
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frequency roll-off inherent to the filter. The filter taps have been 
compensated for this characteristic. 

Alignment Verification 

Before verifying the alignment of the CZT filtering unit, the 
dc offsets in the signal extraction channels should be removed. This 
may be accomplished by adjusting VOl and V03 to zero volts dc when no 
inputs to the I and Q channels are applied. 

With the dc offsets removed, inject Vj and Vq into the I and Q 
channels, respectively. To stabilize the spectral lines appearing at 
VOl and VOS it is advisable to synchronize Vj and Vq with the CZT clock. 
This may be accomplished by driving the alignment generator at a sub- 
harmonic of the sampling clock. The outputs VOl or VOS or both should 
contain a spectrum identical to that of Figure B-2 when monitored with 
an oscilloscope. The oscilloscope gain may be increased to examine the 
intensities of the suppressed lines. If they are present, then the 
CZT filtering unit is not properly aligned. The location of the 
suppressed lines may be ascertained by removing the Q or I input. If 
the output spectrum matches that of Figure B-S, then the input signals 
have been interchanged. 
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APPENDIX C 


THE DESIGN AND OPERATION OF 
THE ALIGNMENT GENERATOR 

In order to emmulate quadrature radar return signals, a signal 
alignment generator was designed and constructed to produce quadrature 
periodic square wave signals. The circuit in Figure C-1 was utilized for 
that purpose. 

The incoming square wave is sent directly to the clock input of one 
D type flip flop (FF) and, at the same time, is complimented by use of a 
NAND gate and sent to the clock input of the other D FF. Due to the fairly 
low frequency of the system input, the delay associated with the NAND 
gate can be neglected. The clear and preset inputs to both FF's are 
disabled by typing them high. Since a D FF is triggered by a rising edge 
on the clock cycle, the effective output frequency is one half of the 
input frequency. The clocking signals to the D FF's are 180 degrees out 
of phase, which yields signals out of the divide-by-two circuitry that 
are 90 degrees out of phase. This is demonstrated in the timing diagrams 
in Figure C-2. The top signal is the input square wave (CLl) from the 
controller board. The second waveform is the complementary clock signal 
(CL2) from the output of the NAND gate. The third and fourth waveforms 
are the outputs (Q^. and Q^) of the in-phase and quadrature FF's respectively. 
Note that the outputs only change on the rising edges of the clock inputs. 

By connecting the complementary output (Q^.) to the input of the in- 
phase FF, the in-phase FF will operate in a toggle mode. The complementary 
output of the in-phase signal (Q^) and the quadrature signal (Q^) are 


NANDed together and the resulting output is used as the input to the 
quadrature FF. The output of the NAND gate only goes low when both 
inputs are high. For that case to be true, the quadrature signal (Q^) 
must be high and the in-phase signal (Q^) must be low (state A), With 
the output of the NAND gate in the low state, the next CL2 rising edge will 
cause the quadrature FF to toggle into the low state (state B). This 
change in the quadrature output to the low state will cause the output 
of the NAND gate to go high. One half cycle late, CLl will produce a 
rising edge which toggles the in-phase outputs (state C). The change will 
not affect the output of the feedback NAND gate therefore the input to 
the quadrature FF will remain in the high input (state D) and the resulting 
output will cause the output of the NAND gate to go low again (state E) 
after the next CLl rising edge and the cycle begins again. The input 
to the quadrature FF can be seen in the fifth waveform displayed in Figure 
C-2. This arrangement will assure that the quadrature (SIN) channel refer- 
ence signal will always lag the in-phase (COS)' channel by a quarter period. 

The output signals from theFF's are buffered by the diode, transistor, 
and resistor circuits before being converted to analog signals by the 
LM301 operational amplifiers (op amps). By careful selection of the re- 
sistors used in the op amp circuits, the output analog signals can be 
balanced. 

In order to provide stable spectral lines at the output of the CZT 
filter, the alignment generator should be synchronous with a sub-harmonic 
of the CZT clock signal. This signal is obtained by tapping outputs from 
the CZT clock signal counter (ZGl) on the controller board. For ease in 
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presentation on the oscilloscope, the fifth sub-harmonic (1/16 of the 
CZT clock rate) was utilized (pin 7). 
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1.2K* 



C-1 SCHEMATIC OF ALIGNMENT GENERATOR 
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APPENDIX D 

PARTS LIST FOR SCAHEROMETER 
SIGNAL PROCESSOR 


I. Controller Board 

Number 

8 

7 

15 

2 

2 

3 

2 

3 

7 

1 

1 

2 

3 

1 

1 


II. Digital Board 

Number 

16 

5 

8 

8 

8 

4 

14 

5 
2 
2 
4 
2 
2 
3 


Device 

SN7474N Dual D flipflop, TI 
SN7408N Quad 2-input NAND Gate, TI 
SN74193N Binary up-down counter, TI 
316A102 Allen Bradley resistor package 
C270B EPROM INTEL 

SN74174N 4 bit bidirectional shift register, TI 

SN7417N Hex Buffer/Driver, TI 

SN7432N Quad 2-input OR gate, TI 

SN7404N Hex Inverter, TI 

SN7430N 8-input NAND gate, TI 

P8205 1 of 8 binary decoder, INTEL 

SN74273N Octal D flipflop w/clear, TI 

P8212 8 bit input/output port, INTEL 

SN7411N 3 input positive AND gate, TI 

SN74LS112N Dual J-K flipflop (neg. edge), TI 


Devic e 

P2111 RAM, INTEL 

SN7474N Dual D flipflop, TI 

P8216 4 bit bidirectional bus transceivers, INTEL 

SN74273N Octal D flipflop w/clear, TI 

SN7408N Quad 2-input AND gate, TI 

SN74LS352N Dual 4-line to 1-line data selectors, TI 

SN74283N 4-bit binary full address, TI 

SN7404N Hex inverter, TI 

ADC-HS12BGC Datel A/D Convertor 

MPY-12AD TRW 12 x 12 Multiplier 

64x503 Spectrol Trim Pot 20 turn 50K 

18M Resistor 

IM Resistor 

IK Resistor 


* Unless otherwise specified: 

All resistors h watt 5"/ 

All Capacitors CTC mono cap 10% 
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III. Analog Board 


Number Capacitors 


27 

18 

17 

12 

5 

4 

8 

4 

4 

4 

2 

2 

2 

2 

4 


0.1 mf 
33 pf 
0.01 mf 
0.001 mf 
5 pf 
560 pf 
470 pf 

470 pf silver mica BOV 5% 

420 pf silver mica BOV 5% 

39 pf 

1.0 pf electrolytic hermatically sealed 
330 pf 
1500 pf 
10 pf 
82 pf 


Number Resistors 


4 

4 

6 

4 

4 

8 

6 

4 

4 

3 
2 

15 

1 

4 
2 

15 

4 

4 

4 

2 

1 

2 

1 

1 

1 

2 

2 

2 


56K 
39K 
21 K 1% 
12. 4K 1% 
18K 
2.7K 
14. 7K V/o 
8.66K n 
39 K 
200 K 

1500 Q n 

IK 1% 
680fi 
5. IK 
lOM 

lOK 1% 
6.8K 
1780^2 1% 

820 n 

1.2K 

470fi 

150Q 

1000 

470 

2000 ]% 

5K 

IM 

2.32K ^% 
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V • ‘f - '.jf - 


',*rr 


Analog Board (continued) 


Number 

Resistors 

2 

7.5K n 

2 

8.2K 

4 

lOOK 

2 

4.7K 

1 

536^2 

10 

.51fl 

Number 

Trim Pots 


6 

64x502 Spectrol 20 turn 5K 

4 

64x103 Spectrol 20 turn lOK 

2 

64x102 Spectrol 20 turn IK 

Number 

Integrated Circuits 

18 

LM 308A Operational Amplifvar, TI 

15 

LF 356P Low Offset Monolithic JFET Input, TI 
Operation Amplifier 

2 

IM 301A Operational Amplifier, TI 

3 

P2708 or P2704 EPROM INTEL 

4 

DAC-HA12BC Datel D/A Convertor 

4 

DS0026CN 5MHz 2-Phase MOS clock driver, NATIONAL 

1 

R5601-1 Reticon Quad Chirp Transversal filter 


Number 

Transistors and Diodes 

8 

SD210, SIGNETICS 

5 

IN4448, TI 

2 

IN703A, TI 

1 

IN753A, TI 

1 

IN756A, TI 


Display Interface 

Number Device 

1 SN7404N Hex Inverter, TI 

1 SN7432N Quad 2-input OR Gate, TI 

1 SN7408N-Quad 2-input AND Gate, TI 

2 SN7417N Hex Buffer/Driver, TI 

1 SN7474N Dual D Flip Flop, TI 

20 AB316A102 Allen Bradley Resistor Package 

20 SN74LS273N Octual D Flip Flop w/clear, TI 

20 TIL305 Alphanumeric Display, TI 

2 P8205 3-to-8 Decoder, INTEL 

2 MTXAl Matrox Alphanumeric Display/ Keyboard Controller 

1 MCI 4490 Hex Contact Decoder, MOTOROLLA 

1 Hex Keyboard CHERRY B65 Series 

4 SN74366AN Tri -State Hex Buffer, TI 
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DRAWING LIST FOR SCATTEROMETER 
SIGNAL PROCESSOR 


Identification Description 


3556.01.01 

.02 

.03 

.04 

.05 

.06 

.07 

3556.05.Q11 

.012 

.021 

.022 


.03 

.041 

.042 

.043 

.044 

.05 

.06 

.07 

3556.06.01 

3556.07.01 


Display, Section 1 
Display, Section 2 
Display, Section 3 

Display Driver, Right Half and Keyboard Control 1 

Display Driver Left Side 

Keyboard 

Display Data Buffer 

Signal Conditioner for the I Channel 
(L~Band) 

Signal Conditioner for the I Channel 
(C-Band) 

Signal Conditioner for the Q Channel 
(L-Band) 

Signal Conditioner for the Q Channel 
(C-Band) • 

S/H and Signal Splitters 
Sine and Cosine ROMS 
Cosine Chirp Premultiplier 
Sine Chirp Premultiplier 
Summing Circuitry 
CZT Filter 

0^ Signal Extraction Circuity 
0g Signal Extraction Circuitry 
A/D Converters and Squares 
I and Q Summer 


Identification 


Description 


3556.08. 

.01 

.02 

.03 

.041 

.042 

.043 

.044 

.051 

.052 

.053 

.054 

.06 

.07 

3556.09.01 
.02 
.053 

3556.10.01 

3556.10.02 


Accunnilator and Memory 

Preprocessor .MDR and Byte Selector 

Preprocessor Memory Data Bus Multiplexer 

Preprocessor Memory Address Register 

Memory 

Memory 

Memory 

Memory 

Multibus Interface Logic and Buffers 

DMA Address Registers 

OD Gate, R/W Gate and E0-E7 Buffer 

CZT Synchronization and Gating 

Accumulator Adder 

Reset Circuitry 

DMA Transfer Controller 

DMA State Controller 

Clock Generator and Clock EPROM 
Address Register 
NERDAS Interface 
Bi-(J) L Interface 


APPENDIX F 


FLOW CHARTS FOR THE PROCESSOR SOFTWARE 


RESET MODULE 


ciPS" 




RESET MODULE, C- and L-Band Scatterometer Data Processor 




Set to accept Interrupts 

For keyboarij entries; RUN & HALT 

LOAD DMA Address Registers 


Accept ; 



Vectored j 
Interrupt^ I 


Interrupt^ 

Processor 
i = 0,2... 
7 



X -j 

DECFG I 
DISFG f 
Enter J 
NNBR = 1 
‘^(Staek) 


= 'b' (blank) 


A 
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RESET Module: 

INI59, Initialize 8259 Interrupt Controller, set Addr. Registers 


^ Entry ^ 


INI59: 


N 


ICWIZ 

N 

L 

OCWl 

z' 


L 

CZTI 



f 


QEI) 


1 . 2.1 


Issue ICWI and ICWZ to DAH & DBH, 

ICWl - A7A6A5 10110 

ICW2 - A15A14...A8 « 1st VECTOR ADDRESS 


Issue OCWl to Loc. DBH, 

OCWl * 02H, Masking the End-of -Convert 
(EOC) or Integration Timer interrupt 


Initialize CZT Address Registers 

Send least significant part of 
’SPWR' ADDRESS to port OAH, and 
most significant part of ADDRESS 
to port OEH 
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1 1 




RESET Module: 

INIBU, Initialize Buffers for E 


c 


1 . 3 . 


(NPASS)-<-o> 
J - o 
k = o 
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RESET Module: 

RFPD, Reset Front Panel Display 


Entry 


RFPD; 


l.A 

Clear Puts ’O' in each led 
of the Display 

Load Display Length Registers 
Addr ■ 12. n 


pg. 17 
RFPDl ; 



STATUS DISPL^, 1-1, 13 


1 

2 

3 

4 

5 

6 
7 
B 


NERDAS Status 
Buffer Length 
Time Cycle 
Altitude Over-ride 
Velocity Over-ride 
Pitch Over-ride 
Roll Over-ride 
Drift Over-ride 


1( 

1 ] 

12 

13 


3 character 
ID Field 
character 


'b' or 


'b' (blank) 
Digit Display 




DISPL 

12 

y 


L 


Addr = 14 


10 


DG^, i=l, 5 (Digit Displays) 


c!d 


1 

2 

3 

4 

5 


• 1 I lot In 

or ’ ’ 


iQf 
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RESET Module 


RESIP , Reset Interrupt Processor 


RESIP: 



f 

I 


^ Entry ^ 


\ 

/ 1 

X 

DECFG 1 
DISFG r=o 
ENTER J 

NNBR+-1 

\ 


(stack^)-^'b ' 

^ 

i 


^ Exit ^ 


Reset Digit Counter, Decimal 
Display, Enter, and Negative 
Flags 


Fill ASCII Character 
Save Area with Blanks 


211 
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INTERRUPT PROCESSOR . Interrupt-Q 

RDKEY; Created by a digit key on Front Panel 



Load Acc. A 
with COMMAND 


Put *^(A) on 
Data Buss 


Acc. A gets 
^(Data) » Key Value 


Save Key Value 
at Loc. INDIG. 
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INTERRUPT PROCESSOR . Interrupt-0 
DECIM; Digit Key Value > 11 








INTERRUPT PROCESSOR . Interrupt-0 
CHS; Digit Key Value»10 







upt-0 

Key Value < 9 
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INTERRUPT PROCESSOR . FARM Module 

DISPL: Digit Key Value < 9, c (ENTER) i 0, c (ENTER) i 1 


0.16 



PCAR: 


'DRF = 
'ROL = 
• PTC = 
’VEL - 
'ALT - 






Character 
Strings for 
Front Panel 
ID Field 
Displays 


218 



INTERRUPT PROCESSOR . PARM Module 

ERRl; Digit Key Value _< 9, ®(ENTER)-0, ®(X)^5 


0.15 


^ Entry ^ 

^ 

f 

Disable 

Interrupts 

/ 

^ 

f 

CHG Mask 
to ACCEPT 
’CLR 'or 
’RESET’ onlj 


f 

N 

FLA 

SH 


Reset Interrupt 
Controller to Respond 
only to 'CLR' function 


Flash Front Panel 

@ 500 ms Intervals 

until ’CLR’ or ’RESET’ Keyed 



INTERRUPT PROCESSOR. PARM Module 

ACPRM; c (ENTER) j^O, ^ (ENTER) -1, 3< Digit Key Value < 8 


0.21 



Calculate Address 
Pointers for Farm 
Value, ID String, 
and Value String 
Save Location: 

VADDR- PARV+ ( INDIG) -4 ) *2 
ID = PCAR + (c(INDIG)- 4)*5) 
STRNG * STGSV + (c (INDIG)-4)*5) 

Display Appropriate 
Character String on 
Front Panel ID Field 



PARV: 

0 

; Values Store 


0 

• 

; In FLTG Point 


• 

• 

0 

; Form 

STGSV: 

0 

; Character String 


0 

; Saved from Front 


PCAR: 


0 ; 

Panel Entries 

DRF - 

; Character 

ROL = 

; String for 


; F/P ID 

• 

; Field 

ALT=’ 

9 
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INTERRUPT PROCESSOR . Interrupt-0 

DGSTK; Digit Key Value < 9, ‘^(ENTER)-O, *^(X)< 5 






INTERRUPT PROCESSOR . 

RFPDP; Reset Part of Front Panel, Interrupt-0, ‘^(INDIG)<10 
c(ENTER)-0, °(X)-0 


RFPDP 



pg. 4.0 
RFPDl : 







INTERRUPT PROCESSOR. DISPL Module 


WRID; Write ID to Front Panel, Interrupt-0, c(ENTER) i 0, c (ENTER) 
‘^(DISFG) ^ 0, (INDIG) < 9 


0.18.2 




N 

L 

REG H, L-^- 
C(DSPAD), 
CURSOR 
ID Field 

N 

/Q.18.2 

FPRNT 



L 

Reg H,L -t- 
‘^’(PARR) , 
Cursor 
Digit Field 

\ 

f 0.18. 

FPRNT 

i 

L 


Reg. H, L » ADDR of ID string 
to be displayed. Set Cursor 
to Point to field on front 
Panel. Set Reg B ■ NER of 
Characters (4) . 


Write to front Panel 
a Sequence of 1 characters 
from the location » ^(REG H,L) 


Reg. Hj,L * ADDR of Character 
String for Parameter Value 
to be Displayed. Set Curs r 
to Point to Front Panel Digit 
Field. Set Reg. B ■ NBR of 
Characters 


C ) 
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INTERRUPT PROCESSOR, DISPL Module, WRID Module 

FRPNT, Output Character String to Front Panel, 

Reg J,L - ADDR of ASCII String to Display 
Reg B " NBR of Characters 




INTERRUPT PROCESSOR 


INTERRUPT-6, Reset Key Response 



6.0 (- 1 . 1 ) 


INTERRUPT PROCESSOR 
Interrupt-7, CLEAR Key Response. 


GrED 


. N 

( 

Blank 

Digit 

LEDla / 

_ \ 


X 

DECFG 

ENTER 

DISFLG, 

NNRR.«-1 






( RETURN ) 
to 1 .9 V 


7.1 


7.2 (=1.6) 
Reset 

Interrtupt 

Processor 




225 



INTERRUPT PROCESSOR 


DISPLAY, INTERRUPT-4, Function Key *DISPL' Depressed 


4i^Q 

Entry J 


\ 

f 4.1 

^(DISFLG)-*-! 

‘^(ENTER)'^-2 

N 

^ 5.2 

FLSHl 

i 

^ 5.3 


4.2 


Flash Front 
Panel LED's 
Once 


( RETURN ) 
^ to 0.5 y 





INTERRUPT PROCESSOR. PARM Module . ACPRM Module 

SBOR; Set Bit of Over-ride Word, Interrupt-0, c (ENTER) - 1, 
3 < Digit Key Value < 8 


0.21.1.1 


^ Entry ^ 


N 

/ 

Compute 'OR* 
Pattern 
Address = 
PADDR 

^ 

^ 

c(ORFWD) -t- 
‘^(ORFED) .OR. 
‘^(PADDR) 

i 




CjEO 


*0R' Pattern Storage: 


ORFWD: 

ORWRD: 


OOH 

20H 

lOH 

OSH 

04H 

02H 


Over-ride Flag Word 

Drift Mask 

Roll Mask 

Pitch Mask 

Vel Mask 

ALT Mask 


PADDR « ORWRIH-c(INDIG)-4 
Bit Posltlon-BPOS“^(INDIG)-1 


Over-ride Flag Word Bit Assignments 


(MSB) I 0| 




e 

(0 

s. 


as 

o 


c 

CO 


^ 5 


S’ 


H- 

n 

n 

p‘ 


< 

M 

tr* 
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INTERRUPT PROCESSOR 

ENTR, INTERRUPT-5, Function Key, 


5.0 


^ Entry ^ 


N 


c(ENTE 

R)-<-l 

N 


FLSHl 


, 5. 


C Return ^ 

to 0-5 J 


'ENTER', Depressed 


Flash Front Pantl 
LED's Once 

* 0 . 5 ) 



INTERRUPT PROCESSOR 

HALT . Interrupt-3. Front Panel State Key, 'HALT', Deressed 



3.0 


3.1 

Set up Interrupt 
Mask for Interrupt-6 
only Recognition 


3.2 

Build a Fill Frame 
for End of File 
Identification 


3.3 

Write 50 fill 
frames to the 
Output 


3.4 

Reset Front 
Panel 


3.5 


3.6 

Enter a Walt 
Loop until 'Reset' 
Is Depressed 









INTERRUPT-2, Key Labelled ’RUN’ 
RUN Module 



2.0 


Reset Interrupt Mask 
Read A/C NERDAS and 
Set Over-ride 


Set and Start CZT 
Evaluate 

ao' (0L^), i-1, 8 


Wait for end of 
Integration Cycle and 
Interrupt from CZT 


Acquire CZT Data 
and Evaluate 
do (0L^), 1-1,8 


Write time Correlated 
Data to Output 
in Bi (}!-L Format 
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lNTERRUPT-2 PROCESSOR . RUN Module 
INITAL, Key Labelled 'RUN' depressed. 


r 

!• 



Reset Interrupt Mask to 
Accept 'HALT' and 
Integration Timer (EOC) only 


Read NERDAS Into TEMPl/2 
Perform Valldlty/Over-rlde tests 
Set appropriate alarms and 
flags. 





INTERRUPT-2 PROCESSOR. RUN Module. INITIAL Module 


RMSK; Run Mask Set, InterrupC-2, 'RUN' Key Depressed 


2 . 1.1 


^ Entry 

a 

N 

^ 


OCWI 

1 

^ 



^ Exit 

i 


Issue OCWi to the 8259, Location 
OODBH, for Mask; 

OCWI = 00F5H, Masking all but 
'HALT' and 'EOC Interrupts 
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INTERRUPT-2 PROCESSOR . Run Module . Inltal Module 

NERDX; Interrupt-2, SIGMl Module 

Reg(H,L) •> Destination Starting Address 
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INTERRUPT-2 PROCESSOR. RUN Module. INITIAL Modu.le, NERDl Module 
NERDZ: NERDX Routine • SIGMl Module, Interrupt-2 

Decode & Block A/C Data into TEMPI or l'EMP2 Buffer 
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Storage Arrangement; 

TEMPI: 0 ; BCD Values-of 

0 ; A/C Data, Each 
0 ; WORD = 4 LSB's 

■ 

• 

0 ; Of 107 WORDS 
0 ; Of ACBLK, Thru 
0 ; Run Number 

TEMP2: 0 ; BCD Value of 

0 ; A/C Data, each 
0 ; W0RD-4LSB's 

• 

• 

0 

0 

0 ; 




'i. 

Hi: 



INTERRUPT-2 PROCESSOR . RUN Module . INITIAL Module . NERDl Module 
VALID; NERDl Routine, SIGMl Module, Interrupt-2 
Converts Time & Other A/C Parameters 

to FLTG Point. Where Appropriate, Over-ride Values are 
Substituted & Alarms & Flags are Set. 


^ Entry ^ 


2 . 1 . 1. 2 


^ 

^ 

TCONV 

^ 

/ 

HFLTG 

^ 

f 

DFLTG 

N 

L 

FROLL 

\ 


FPITC 

\ 

f 

FGSPD 

\ 



Time Is Converted 
from BCD to FLTG 
Point Value <SECS*10) 
Saved at FLTME(<10°) 

Altitude Is checked & 
Converted to FLTG Point 
& Saved at Alt In Feet, 
Unless Over-ride Set 


If Not Over-ridden 
Drift Is Checked & Saved 
as Radians, FLTG Point. 


If not Over-ridden 
Validate Roll, CHG to 
FLTG Point & Save (Radians) 


If Not Over-ridden, Check 
Convert Pitch to FLTG Point 
(Radians) & Save 


If not Over-ridden, Check & 
Convert GND Speed to FLTG 
Point (Knots) & Save 0 VEL 


C 


FLTMB: 

ALT; 

DRFT; 

ROLL; 

pitch: 

VEL: 




;Cell Time 
;Altitude 
;Drift (RAD) 
;Roll (RAD) 
;Pitch (RAD) 
;GND Speed 
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SIGMl ; Key Labelled 'RUN* depressed 


C 


Entry 






Pg 38 
SIGM2 


2 . 1.2 



Evaluate Cell Sizes 


L^, 1 - 1, 8 


2.1.3 Compute Integration time, 

^ " T1 Compute Number of Sub-Recoras, NS 



Set and Start the 
CZT board Integration 
and timer. 


2.1.4, 2.1.5 

Compute Buffer Length, NL 
Initial or Adjust 
Buffer Poiners for Loading 
And for Write-Outs. 


Evaluate Sub-Parameters 
And Then The Term 


ao'^, i * 1, 8 


Set c(WAIT) = 1, Loop Until 
CZT Integration Timer 
Interrupt-1. 



SIGM2 
Pg 36 
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INTERRUPT-2 PROCESSOR . Run Module . SIGMl Module 


STINT; Start Integration 


2.1.3 


^ Entry 

3 

N 

L 


STINT 

1 



V 


^ Exit 

3 


Evaluate 

NS=[(Tk/512)f8+0.5]j 

Send Value of NS to 
CZT Board, starting 
Integration 
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T.NTERRUPT-1 PROCESSOR . RUN MODULE 

SIGM2 ; End of Integration Time (EOG) , Or 
Exit From WLOOP With lEOC Set. 


1 . 2.1 
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Transfer CZT Spectral 
Data Into p-P Memory 


Evaluate Calibration 
And Polarization Power 
Set Flag Bits, And 
Alarms If Appropriate 

















INTERRUPT-1 PROCESSOR . RUN Module 
SIGM2 ; (Cont. B) 



Write Appropriate Output 
Line To The Bi-Phase-L 
Interface, Set ^(NPASS)“1 


Read A/C Data 8. 
Clock Time, T , Into 
TEMP2 ^ 


Calculate Cycle Time, TC 
Move Ti, 


(TFLG) ^0 
Clear From 
Front Panel Time Alarm 


Pg31 


SIGMl 


'tW/\ 























INTERRUPT-1 PROCESSOR . RUN Module . SIGM2 Module 
I»1A; Initiates DMA Transfer of CZT FREQ Components 


^ Entry ^ 


■ 

■ 

OUT CZT 2 


■ 



4.3.3 


^ Exit ^ 


CZTl EQU 08H 
CZT2 EQU 09H 
CZT3 EQU OAH 
CZT4 EQU OBH 


SPWR:0 

0 

0 

• 

• 

0 

0 
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;512 X 32 
;BIT 

; BUFFER FOR 

• 

• 

; SPECTRAL 
; POWER 







jEE 


liEC- 




Entry 


Using Constant Filter Index 
MC-[FCAL * TR + 0.5] , 
FILTL II Evaluate PT-SPWRC^^ (MC) *32) 


Where CR*correct Ion factor 


MN-{0.6*MC}^, 

Evaluate PN-SPWRC*^ (MN) *32) 


PTR-PT/PN 


PTR;STD 


CALRM 


FILTP 


Evaluate Filter Index 
For Polarization Tone 
R+0 
S 

Evaluate PTP=PP/PN 


Set Polarization 


Flags and 


Gain Flags 


PTP : STD 


HPOLZ 


VPOLZ 

















lNTERRUPT-1 PROCESSOR . RUN Module . SIGM2 Module 
POWER; Computer Spectral Power for Each Doppler Band 



SPWR:0 ; Spectral 
0 ; Power 

: ;Buffer 


; 512x4 Words 


PR:0 ;Received Power 
0 ;For Each 
0 ;Doppler Band 


0 ;Up to 8 Values 










INTERRUPT-1 PROCESSOR . RUN Module. SICM2 Module 
GAIN ; Determine Gain Values From Table (L-Band) 


^ Entry ^ 


\ 

L 

TRCVT 



L 

GLX 





CZED 


.3.7 


Determine Type of 
Transmit/Receive 
Configuration; I. E. 
'Like' or 'Cross' 


Calculate Gain 
Table Index As 
Function of Polz 
And ICON (like, cros 






INTERRUPT- 1 PROCESSOR . RUN Module . SIGM2 Module 

A A 

SIGB ; Computes or ao^^^ For Each Doppler Band 
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LOAD ; SIGM2 ROUTINE, INTERRUPT-1 

Computes Pointers Into Out-Buffer, Moves 
Nerdas, 5o(0|^) Into Appropriate Lines 
And Columns. Moves Balances of Other 
Out-Put Parameters Into Out-Buffer. 



.3.10 


Evaluate Line Pointer For NADIR & 

For Each Look Angle 

LP =(DF/l / 256)*N,k=l, 8 

lire. -M 0D( (LDFTR-LP. ) , (N+10) ) ,k-l , 8 

LINE -MOD (LDPTR, (N+10) ) 
n 


Move a (0 ) to LINE , Col Of The 
Output Buffer. Also, Move Values 
of PT, Polz, Alarm Words, Over-Ride 
And Flags To Output Buffer LINE 


Move Nerdas Data Into 
Approrlate Segment of Line 
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LOAD ; SIGWa INTERRUPT-1 

Computes Pointers into Out-Buffer, Moves 
Nerdas, do(a|^) Into Appropriate Lines 
And Columns. Moves Balances of Other 
Out-Put Parameters Into Out-Buffer. 



4.3.10 

Evaluate Line Pointer For NADIR & 

For Each Look Angle 

LP =(DFA. /256)*N,k=l, 8 

LINE. »M0D( (LDPTR-LP. ) , (N+10) ) ,k=l , 8 

LINE =M0D(LDPTR, (N+10)) 


Move a (0 ) to LINE^, Colj^ Of The 
Output Buffer. Also, Move Values 
of PT, Polz, Alarm Words, Over-Ride 
And Flags To Output Buffer LINE 


Move Nerdas Data Into 
Approriate Segment of Line 
















OUTPUT BUFFER LINE CONTENTS 



BYTE 

(BCD) 


BYTE (BCD) 

ADDRESS 

CONTENTS 

ADDRESS 

CONTENTS 

XX 00 

*0.1 


. XX 13 

*0.1 


01 

*1 

SEC 

14 

*1 

LONG, 

02 

*10 


15 

*10 


03 

*1 

MIN 

16 

*1 


04 

*10 


17 

*10 

LONG, 

05 

*1 

HRS 

18 

*100 


06 

*10 


19 

SIGN 


07 

*1 

DAY 

lA 

*1 


08 

*10 


IB 

*10 


09 

*1 

MONTH 

1C 

*100 

ALT 

OA 

*10 


ID 

*1000 


OB 

*1 

YR 

IE 

*10000 


OC 

*10 


IF 

*0.1 


OD 

*0.1 


20 

*1 

HDG 

OE 

*1 

LAT, MIN 

21 

*10 


OF 

*10 


22 

*100 


10 

*1 


23 

*0.1 


11 

*10 

LAT, DEG 

24 

*1 

DRIFT 

12 

SIGN 


25 

*10 





26 

SIGN 



MIN 


DEG 


OUTPUT BUFFER LINE CONTENTS (CONT.) 



BYTE 

(BCD) 


BYTE 

(BCD) 

ADDRESS 

CONTENTS 

ADDRESS 

CONTENTS 

XX 27 

*0.1 


XX 3C 

*0.1 


28 

*1 

ROLL 

3D 

*1 

THETA-A(l) 

29 

*10 


3E 

*10 


2A 

SIGN 


3F 

*0.1 


2B 

2C 

*0.1 

*1 


40 

41 

*1 

*10 

THETA-A(2) 

PITCH 


2D 

*10 


42 

*0.1 


2E 

SIGN 


43 

*1 

TAETA-A(3) 

2F 

*1 


44 

*10 


30 

*10 

GND 

SPEED 

45 

*0.1 


31 

*100 


46 

*1 

THETA-A(4) 

32 

*1 


47 

*10 


33 

*10 

MISSION 

48 

*0.1 


34 

*100 


49 

*1 

THETA-A(5) 

35 

*1 


4A 

*10 


36 

*10 

SEQ 

48 

*0.1 


37 

*100 


4C 

*1 

THETA-A(6) 

38 

*1 

FLT ■ 

4D 

*10 


39 

*10 


4E 

*0.1 


3A 

*1 

LINE 

4F 

*1 

THETA-A(7) 

38 

*10 


50 

*10 
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OUTPUT BUFFER LINE CONTENTS (CONT.) 



BYTE 

(BCD) 


BYTE 

(BCD) 

ADDRESS 

CONTENTS 

ADDRESS 

CONTENTS 

XX 51 

*0.1 


XX 64 

*0.1 


52 

*1 

THETA-A(8) 

65 

*1 

SIGMA(5) 

53 

*10 


66 

*10 


54 

*0.1 


67 

SIGN 


55 

*1 

SIGMA(l) 

68 

*0.1 

SIGMA(6) 

56 

*10 


69 

*1 


57 

SIGN 


6A 

*10 


58 

*0.1 


6B 

SIGN 


59 

*1 

SIGMA(2) 

6C 

*0.1 


5A 

*10 


6D 

*1 

SIGMA(7) 

5B 

SIGN 


6E 

*10 


5C 

*0.1 

• 

6F 

SIGN 


5D 

*1 

SIGMA(3) 

70 

*0.1 


5E 

*10 


71 

*1 

SIGMA (8) 

5F 

SIGN 


72 

*10 


60 

*0.1 


73 

SIGN 


61 

*1 

SIGMA(4) 

74 

*0.1 


62 

*10 


75 

*1 

PT 

63 

SIGN 


76 

*10 

(CALIB. PWR) 




77 

SIGN 





256 


OUTPUT BUFFER LINE CONTENTS (CONT.) 


BYTE (BCD) 

ADDRESS CONTENTS 


XX 78 
79 
7A 
7B 
7C 
7D 
7E 
7F 


Polarization (HH-llH, HV=10H, VV-OOH, VH«01H) 
BUF Alarm (0=0FF, l=ON) 

Time Alarm (0=OFF, l»ON) 

Calib. Alarm (0»OFF, 1-ON) 

Over-Ride Bits 
Last-Value Bits 
FFH 

End of Frame 

FFH 
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APPENDIX G 




APPENDIX G 

INTEL SBC 80/20 SOFTWARE 


•09» M4CM0 ASSSMBUF** VCR 


ERRORS ■ O'PACF I 


REMOTE SENSING CENT«» 
TEXAS *m UNIVERSITY 

PROJECT 3S5fi 


OESCRIPTION 


THIS PROGRAM CONTAINS ALL SOFTWARE UFVElOPED 
FOR the IN’fFL SBC SO/20 DURING PRnjFC- 3S66. 
THERE ARE 5 PROGRAMS WITHIN THE MAIN 
PROGRAM. UPON EXFCU^EION OF THE PROGRAM, 

A MENU IS PRINTED. THE USFR MAY SELECT THE 
DESIRED program crcm thf MENU. THE SELECTED 
PROGRAM WILL EXECUTE AUTOMATICALLY, PRf'MPTiNG 
THE OPERATOR FtlR ANY AODITIONAL REQUIRED 

information. 


THE S PROGRAMS ARE: 

II 'SAMPLE ACCUMULA'^ION AND DMA TRANSFER* 
this program initialises THC CZT 
HAPOWARR, performs SAMPLE ACCUMULATION, 
AND ODES TMR 0«A tp4m«jfep tO THE S3C 
BO/20 AT M*-MCRY location 6000 HEX. TO 
terminate the program, hit *ESC* WHEN 

prompted for number or records, 

21 'SAMPLE accumulation. OMA TRANSFER AND 
transfer to ti« 

THIS OROGRAM deRRORMS EXACTLY AS RROGRA« 

I, EXCEPT, AT the end of THE OMA TRANSFER 
TH*' data is TRANSFEREO TO THE TI OAO MINI- 
COMPUTER. the program 'TX90S0' MUST 
BE RUNNING ON THE TI RQR A TRANSFER TO 
TAKE PLACE. 

SI 'ADC ALIGN' 

THIS PROGRAM IS USED t<3 ALIGN THE ANALOG 
TO DIGITAL CCNVSRTERS IAOO. IT CAUSES 
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•9t) N4C90 AtSCNBLC** VER >•« 


ERRORS ■ 0 PACE 2 


the processor to T*KE 2S« records • DO 

A DMA transfer »N0 INMEOJATFLY START 
TAKING MORF qrCOROS. SINCE THE DMA 
REQUIRES ONLY ABOUT SO MSEC TO COMOLETE 
IT gives the impression THAT THE PROCESSOR 
IS IN A MODE OF CON^’INUOUS SAMPLE 

accumulation. terminate the 

PROGRAM, PRESS •ESC*. 

41 ‘CASSETTE load* 

▼HIS PROGRAM mill LOAD PROGRAMS OR 
DATA STORFO ON CASSET-^E TARES FROM THE 

silfnt -*00 Terminal -o the ssc so/zo 
MEMORY. cassette drive I MUST BE USED. 

AND THE UNIT must RE I N THE PLAYBACK 
MODE. THE OPPRATOR MiJST POSITION THE 
TAR* TO the beginning OF THE DESIRED 
PROGRAM, the program PPOMOTS the OPERATOR 
FOR •'HE MEMCRY ADDRESS MHERE THE PROGRAM 
ON tape is to be STOPED. UPON COmo(_eTION 
CONTROL IS returned •O ThE MONITOR. 

SI aCASSETtc puKCH* 

this program mill PFCORD programs from 

COMPUTER MEMORY TO CASSETTE TAPE. 
cassette drive I must be USED. AND THE 
UNI^ MUST 0E IN pe-CORO MODE. THE 
QPFraTOP IS prompted for STARTING AND 
ENDING addresses OF DATA OR PROGRAM TO 
BE RECORDED. THE •OECOED* REFEREO TO 

IN the operator instructions is the 

ON/OFF SWITCH IN the RECORD CONTROL 
SECTION OF THE SIL«NT TOO. UPON COM- 
PLETION, control is returned to the 
MON|TOP» the RECORO control MUST PE 
TURNED OFF BEFORE THE TAPE CAN BE MOVED, 


PROGRAM equates 


©OFF 

ICWt 

EOU 

OFFM 

090B 

ICM2 

equ 

OBH 

0999 

DMA 

EOU 

09H 

0020 

OCH2 

EOU 

20H 

0000 

LOAON 

EOU 

OSH 

OOOA 

LSB 

EOU 

OAH 
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•OBO M^C^O ASSCMBUe** VER 7.4 


(O' 


Cl 


009B 

MSB 

EOU 

9BH 

9727 

CETMX 

EOU 

02»7M 

0914 

MON 

FOU 

OOIAM 

92 SC 

Gf TNM 

FOU 

9ESCM 

0220 

6ETCH 

EOU 

9220H 

oica 

CO 

EOU 

OIERM 

OIOF 

CNVBN 

EOU 

OlOFM 

0IC2 

BREAK 

EOU 

OICPH 

02AI 

HILO 

EOU 

92A1M 

92CS 

NMOUT 

EOU' 

92C4H 

0IS5 

STMLF 

FOU 

0355H 

0)00 

CR 

EOU 

OOH 

OOOA 

LF 

EOU 

OAH 

0097 

BELL 

EOU 

9TH 

907F 

RUB 

EOU 

7FH 

3FF9 

TEMP 

FOU 

3FF9H 

OS 90 

t 

ORC 

0S09M 


INITHLIZATtON OF THE SZE? INTERRUPT PROCESSOR 


OSOO 

0300 

0301 

'3 

35FF 

I 

Int: 

01 

MVt 

A.ICRt 

lOISABLE the 8930 INT«^PRUPTS 
;SETS N-^STED MODE INTEORUPTS. 

0303 

0 333 


nuT 

903H 

SARTTE ICYl TO 8239 

0305 

3£0B 


MVI 

AtICwa 

IPUTS CALL TABLE AT 03E0 

0307 

333R 


OUT 

0O9H 

JARITE ICW2 TO 8259 

0009 

FB 


El 


;ENA8LE 8930 INTERRUPTS 


4 BYTES APART 


PROGRAM SELECTION ROUTINE 


03 OA 


sel: 




030 A 

06C8 


MVI 

QfMSILG 

: PR int menu 

03 OC 

218B09 


LXI 

MtMRSG 1 

: 

93)F 

COSI 09 

. 

CALL OUTHS 

• 

0312 

C02792 

enter: 

CALL GETmX 

5GF.T PROGRAM NUMBFP 

0315 

79 


MOV 

AtC 

;mcve to accumulator 

0316 

32730S 


STA 

LOC 

;s»ORF FOR later use 

0319 

FEOI 


CPI 

01H 

Its IT PROGRAM 17 

03IB 

CA3003 


JZ 

INIT 

IIP Se, JUMP TO CZT 

031E 

FE02 


CPI 

02H 

IIS IT PROGRAM 27 


INITIALIZATION ROU''INE 




oS 




•0(19 MAC^O ASSENBLCHt VeM <*A 


o«zo 

CA3D0a 

JZ INIT 

0923 

FE03 

CPI 0 3H 

•425 

CASA >>9 

JZ ATOO 

9920 

FEOA 

CB| 94H 

992A 

CACA09 

JZ CL 

9920 

FEOS 

CPI OSH 

092F 

CAITOO 

JZ CSOUN 

0932 

061 C 

MVI B.MSFLC 

0934 

2IS30A 

LXI H.MESC2 

0937 

COS 109 

CALL OIJTms 

093A 

C31209 

JMP ENTER 


;|F SOt JJMP TO CET INITtALUZATION OOU^INF 
;IS I'*’ OOOGPAM 37 
S IF V*St JUMP TO ADC F.OUTINE 
SIS XT PP3GRAM A7 

;|F SO JUMP TO CASSETTE LOAOEP ROUTINE 

SIS IT PROGRAM S7 

SIF SO JUMP TO cassette PUNCH ROUTINE 


sre-entfo program selection 


CZT INITI ALIZAticn routine 


0430 


INITI 

0930 

C05B09 

CALL LOOMA 

044 0 

0623 

MVI B.MS3LC 

9942 

2I6F0A 

LXI H.MESG3 

0945 

CD6I09 

CALL OUTMS 

994 9 

C02T02 

CALL GETHX 

094B 

79 

MOV A.C 

094C 

0300 

OUT LOAON 

044E 

06IB 

MVI B.MS4LG 

0950 

2I920A 

LXI H.MFSCA 

0993 

CD6I 09 

CALL OUTMS 

0956 

76 

MLT 

0457 

00 

NOP 

0959 

C33099 

JMP INIT 


SLOAO DMA DESTINATION AORESS 


;i»CVF TO accumulator 

SST-'PE ON CZT controller BOARD 


S*AIT FOR INTERRUPT 


SUPON RETURN FROM INTERRUPT. CO BACK INIT 


LOAD DESTINATION ADPESS ROUTINE 


0953 


looma: 


0958 

3E00 

MVI 

A. 00 

0950 

D30A 

OUT 

LSB 

095F 

3E60 

MVI 

A. 6 OH 

0961 

O30B 

OUT 

MSB 

0063 

C9 

RET 



5SETS LS9 DF AORESS TO 00 
SAND MSB TO 60. THEREFORE 
5TME DMA OESTINA'riiON ADRESS 
SIS 6900H. 

S 




1 


•Ot> MAC40 ASSCMBtFR* VKR 2.4 


CKROaS n 0 PAGE 5 


AOC POOTINE 


0464 

1 

ATOOt 

9064 

CD5404 

CALL LDOMA 

9067 

3EFF 

MVI A,8FFH 

0969 

0308 

OUT LOADN 

096B 

COC20I 

CALL break 

096E 

OAI404 

JC MON 

0471 

76 

t 

HLT 


;SETS NO. OF PECOPOS To 256 
ZtOAOS NUMBER OF P='COROS 
5SFF IF fJDPPATOR WANTS OUT 
;IF so CO TO MONITOR 


SUCCESSFUL INTEORUPT PCUTINE 


0472 


good: 



0472 

0309 


OUT 

DMA 

0S74 

3A730B 


LOA 

LOC 

0877 

FE03 


CPI 

03H 

0S79 

C26A0S 


JNZ 

SAMP 

0S7C 

3E20 


MVI 

A,OCW2 

947E 

0336 


OUT 

9D9H 

0680 

FB 


El 


99SI 

C36A09 


JMP 

ATOD 

0984 

FEOI 

samp: 

CPI 

OIH 

0996 

C291 09 


JNZ 

tpans 

0689 

3E20 

back: 

MVI 

A ,OCW2 

099B 

0336 


OUT 

008H 

0960 

FB 


El 


066E 

C33O09 


JMP 

INI T 

0691 

C0A209 

TRANS: 

CALL 

, TX 

0694 

C3690S 


JMP 

BACK 


istart oma transfer 

SRFCALL which PROGRAM WAS LOADED 

!ie IT 37 

! IF NOT GO TO SAMP 

i IF SO, S=’ND 6259 A NON- 

SSPECIFIC end of INTFpdupT ( 0CW2I 

IRFENAfiLE 4080 INTERRUPTS 

ICO BACK TO APC routine 

IWAS IT program 17 

SIF NOT, GO TO TRANS 

;IF SO, SEND 8259 A NON- 

ISPECIFIC END OF INTERRUPT (OCW21 

IREENabue 8080 INTERRUPTS 

;go back tq init 

iPQOGRAM selected was 2, OO A 
5TRANSFER TO Tl, THEN CO TO BACK 


FALSE INTERRUPT RCUTINE 


0697 

0997 

0612 

BAO: 

MVI B,MS5LG 

0699 

2IA00A 


LXt HtMESr.K 

069C 

C06I09 


CALL OUTMS 

099F 

C31409 

t 

JMP MON 


;CO BACK TO MONITOR 


i" ■It; _ 


■aao M4CRO asscmble** veo z.a 


•nnms > o pace a 


> 


TRANSFER DATA TO Tl ROUTINE 


0SA2 


Tx: 




BSA2 

3ECE 


MVt 

AtOCEH 

;SET UP SERIAL PORT 

OBAA 

0 3E0 


OUT 

OEDH 

;F0R P690 BAUD 

99A6 

3E2T 


MVI 

A .2Th 

t 

9BA8 

03ED 


OUT 

OEDH 

8 

• 

OBAA 

2I70IT 


txt 

M.S903 

:loao start address of data 

OBAO 

OBEO 

place; 

IN 

9EDH 

;VAIT until TRANSMITTER READY 

OBAF 

EiiOl 


an: 

91 

i 

3981 

CAAOOS 


32 

PLACE 

t 

OSSA 

77 


MOV 

M.A 

;SENO DATA BYTE TO TI 

093S 

u 

111 

o 


OUT 

9ECM 

• 

• 

0937 

23 


INX 

M 

ICO TO NEXT BYTE 

0938 

3E6B 


MVI 

A.99H 

IREPEAT IF not LAST ONE 

093A 

BC 


CMP 

H 

• 

• 

090B 

C2A008 


JNZ 

PLACE 

I 

OBOE 

OBEO 

wait: 

IN 

OEDH 

IWAIT till TRANSMIT BUFFER 

09C0 

e»OA 


ANI 

OA 

SIS FMPTy 

09C2 

C ABE 08 


32 

BAIT 

• 

OBCS 

3EA0 


MVI 

AtAQH 

SRESET SERIAL PORT 

99C7 

03E0 


OUT 

OEDH 

• 

a 

0BC9 

C9 


RET 


IRE TURN 


CASSETTE LOADER ROUTINE 


09CA 

CL 

• 

• 

OSCA 

3636 

MVI B.MSBLG 

OSCC 

2IBF3A 

LXI H.MFSG6 

09CF 

COB 109 

CALL OUTMS 

0902 

3S«F 

MVI aoffm 

090A 

32F93F 

STA TEMP 

0907 

C02702 

CALL CEtmK 

030A 

SO 

MOV n«B 

09OB 

59 

MOV E.C 

090C 

OEIO 

MVI C«I0H 

OBOE 

C0E801 

CALL CO 

09EI 

3E36 

MVI Ci36H 

09E3 

C0E90I 

CALL CO 

09£6 

OEM 

MVI CtllH 

9BEB 

C0E90I 

CALL CO 


;PRINT INSTRUCTIONS FOR PROGRAM 
• 

jsET for UPPER NtORLE 

icET LOAD AORESS 
JPUT IN OE 

icSI TO PLAY 
( 

; 

t 

I 

:flavback on/kevboaro off 


•Otr MAC»0 ASSEMSLRX* Ve<9 2*4 


lOi 


CF r 


€' 


OSEB 

CO2902 

starti 

CALU GETCH 

IVAIT FOB START SYMBOL 

OSEE 

FE3F 

CPI •>• 

I 

OSFO 

C2EnOS 


JHZ START 

• 

t 


C02002 

GET I 

call GETCH 

;CET A CHARACTER 


F£TF 


CPI tch 

t IGNORE RUB OUTS 

9tfB 

CAF309 


JZ OFT 

1 

9tFB 

FE3C 


C»l •<• 

I STOP |r end symbol 


CAI *90 


Jt MON 

IGC B*CK TO monitor 


COBF9I 


CALL CNVPN 

ICCNVFRT TO BINARY 


*F 


MOV C.A 

1 

T r 

CD5503 


CALL STHLF 

SSTCoe IN MFMORY 

T "* T ‘^1 

3AFP3F 


LOA TEMP 

IINCOCMENT TO NEXT NIBBLE 

. 

d7 


OPA A 

I 

^ r ?fl 

C20F09 


JNZ SKIP 

t 


13 


INK 0 

I 

' '^1 

EEFF 

skip: 

KOI OFFH 

I 


32F93F 


5TA TEMP 

• 

s 

ml 

C3F306 

t 

UMP GET 

IREPEAT 


cassette pumch <recopoi routine 


0917 
091 T 

067E 

C SPUN I 

MVI B.MSTLG 

IPRINT OPERATOR INSTRUCTIONS 

0919 

2IFS0* 


LXI H.MPSGT 

9 

09IC 

COO 109 


CALL OUTMS 

9 

09IF 

OE02 


MVI C.02 

ICFT start and STOP AORESSES 

0921 

C0SC02 


CALL GETNM 

0 

092* 

OEI 0 


MVI Ct lOH 

;CS! TO RECORD 

0926 

COcBOt 


CALL CO 

9 

0929 

0E35 


MVI Ct3SH 

I 

092U 

COEflOl 


CALL CO 

• 

092E 

0EI2 


MVI C.I2H 

IRFCCRO CN/PRINTFR OFR 

0930 

CDEBOl 


CALL CO 

• 

t 

9933 

0E3F 


MVI €.•>• 

; MR ITS start symbol 

0935 

COE90I 


CALL CO 

: 

0939 

01 


POP 0 

;oUT STOP AOPESS IN 0 e 

0939 

Et 


POP H 

;PUT start adress in hl 

093A 

7E 

again: 

MOV A.M 

;get a data byte 

093B 

COCS02 


CALL NMOUT 

;aoite it out 

093E 

COA102 


CALL HILO 

5 REPEAT TILL DONE 

09*1 

23 


ink H 

t 

09*2 

C23A09 


JNZ AGAIN 

• 

0 

09*S 

0E3C 


MVI Ct*<* 

I WRITE END SYMBOL 

09*7 

COEBOl 


call CO 

I 

09*A 

0613 


MVI Ctl3H 

IWRITC DC! (PLAYBACK OFF 1 

094C 

CDEBOl 


CALL CO 

t 










MAC^O ASSFMSLRtt* Wn F«A 


CMOffS « 0 PACr • 


09*f 

OEI A 


HVI C* I4H 

tPECOBO r*r/PK INTER ON 

0951 

COE SOI 


CALL CO 


0954 

OEIO 


HVI C«IOH 

ICSI TO PLAV 

0956 

coEsni 


CALL CO 

I 

0959 

0E96 


HVI C«36H 

1 

095S 

COE SCI 


CALL CO 

i 

0956 

Clt409 

t 

1 

JMP MON 

ICC TO MCNITOR 



I 

t MES5406 PqtNTING 
I 

routine 

0961 

0961 

46 

1 

OUTMS: 

HOV C«M 

IMCVF character to C 

0962 

79 


MOV 4.C 

lANO TO A 

0963 

FEOO 


CPI C9 

S5«E IF CHARACTFR IS 4 RETURN 

0965 

C47I09 


32 C9P0H 

I IF SO. JUMP TO COPUB 

0966 

COE90I 


CALL CO 

I IF HOT OUTPUT character 

0963 

23 


INX H 

IINCRFycHT POINTER 

096C 

05 


OCP B 

ICECPpMENT counter 

0960 

C26I09 


3HZ OUTHS 

SJUMP BACK IF NOT FINISHED 

09T0 

C9 


OF’’ 

I IF FINISHED . return 

0971 

C0ES9I 

C99UB: 

CALL CO 

ICUTPUT carriage rfturn 

09T4 

0E7F 


HVI CtPUB 

IFUT ASCII RUB OUT IN C 

9976 

COE90I 


CALL CO 

; OUTPUT 5 OUOOUTS. 

097 9 

COE90I 


CALL CO 

;iT PREVFNTS THE SILENT 700 

097C 

C0E90I 


CALL CO 

;frch MISSING ANV characters 

J97F 

CDECOI 


CALL CO 

;OUE TO ITS SLOV 

0952 

COES9I 


CALL CO 

ICAORIAGE RETURN 

0985 

23 


INX H 

IINCREMENT POINTER 

0956 

05 


OCP 9 

lOECOFMENT counter 

0997 

C26109 


JNZ OUTHS 

IJUMB BACK IF NOT FINISHED 

099A 

C9 

t 

SET 

I IF finished RETURN 



0 

• 



0993 

0993 

00040453 

: 

HE SGI : 

03 CP.LFtLF, 

•SELECT program: • .CR.LF 

099F 

454C4543 




0993 

54205052 




0997 

4F47524I 




0993 

40340004 




09 9F 

31203020 


09 *1 ■ sample ACCUMULATION ANO DMA TPANSFER 

0943 

53414050 




0947 

4C4S204I 




0940 

43435540 
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•0*0 MAC40 ASSEMBUen. VE» t«4 


VBROBS ■ 0 RACE « 


0«AE 
0983 
0BB7 
0988 
09BF 
09C3 
09C7 
09C9 
09C0 
0901 
0995 
9909 
0900 
09EI 
09CS 
0**9 
99E0 
0971 
0975 
097 9 
9970 
OAOl 
0A95 
9A99 
9A90 
OAI I 
0AI2 
OAIB 
OAI A 
9AIE 
9A22 
9A25 
0A29 
0A20 
9A3I 
0A35 
9A3 9 
0A3O 
OAAI 
0AA5 
0A49 
9A40 
OASI 
09C0 


9A53 

9A53 

0A57 

0A5B 

0A57 


55404154 

49474*20 

41474420 

44404120 

545241 4C 
53444552 
0004 

32201020 
S34I40S9 
4C4S204I 
43415540 
554C4I54 
49474E2C 
44404120 
S4S24 145 
53464552 
2C204I 47 
4426S452 
4I4F5146 
45522754 
47295449 
0O0A1320 
30204144 
43204I4C 
49474590 
OA 

34203020 

43415353 

45545445 

204C474I 

440DOA 

35293020 

43415353 

45545445 

20595547 

4 345000 A 

45534320 

30205245 

54555245 

205*4720 

4047474Q 

5**75290 

0A07 

MSILC 

I 

i 

MESC2: 

00*94756 

4I4C4944 

204S47S4 

S2S92C5* 


OB *t a S4M7LE ACCU9UUAT1 ONrOMA TRANS77^, aNO TRANSFFo *3 


08 CR.L7,*3 a AOC AL t GN* t CR 1 1.7 

08 »A a cassette LOAO*«C»tLF 

08 *9 a CASSETTE 0UNCH*,CRtL7 

08 ‘ESC a return TQ MONITOR •• CR tL7 .CELL 

’20U fMESCt 

08 CRt*INVAL10 EMTRVtrov AGAIN. • «CR «LF .OEM. 




oo^ 
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•••9 M«C^O ASSCMBLERt VE* 


ERRORS B O.PACr 10 


0463 

0A67 

0460 

03IC 


0A6F 

946F 

0473 

0477 

0470 

047F 

04S3 

0497 

9490 

049F 

0003 


0402 
0492 
0496 
94 94 
94 9E 
9442 
0446 
0444 
9910 


OOSO^OAI 

474I494E 

2cOOAA37 

MS2LG 

t 

t 

ME scat 

0O4C4F4I 

44204E55 

4 ^ 424^(92 

OO'iFASO') 

5249434F 

52449320 

29293223 

4S4SS929 

20200 "» 

MS3LG 


EOU ••MESG2 


OB CRt *1.040 AUMBER OF RECORDS 


EOU S-MESC3 


(2 HEA> 


I 

MESG4: 

90414343 

5540994C 

41S4494F 

4E20494C 

2050924F 

47&245S3 

S3000A 

MS4LC 


00 CRt'ACCUMULAVtON IN RROCRESS* tCRiLF 


EOU ««MESG4 


• lOELL 


0440 

0440 

0401 

9405 

9489 

0400 

0012 


04&F 

040F 

04C3 

0AC7 

OACB 

OACF 

0403 

0407 

0400 

040F 

9AE3 

0AE7 


MESCS: 

004641 4C OB CRt*F4LSE INTERRUPT • ,CR .LF 

53452049 

4E5A4552 

52555054 

009A 

MSSLC EOU S-MESC5 
t 


MESG6: 

0O454E5A OB CR<*ENTER ADDRESS WHERE PROGRAM IS TO BE LOADED ( 

45522041 

44445245 

53532057 

43455245 

2050524F 

47524 1 40 

20495320 

S44F2042 

45204C4F 

41444544 


4 HEXI 
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•0*9 NAC40 ASSEMBLER • VEA »«A 




CE 11 


OAEO 

OAEF 

OAES 

0936 


OAES 

OAFS 

0AF9 

OAFD 

OB)I 

2 995 
309 
0300 
031 1 

oat S 

0819 

0810 

0321 

9325 

0329 

0323 

032F 

0333 

0337 

9833 

083F 

0343 

0847 

0343 

0B4F 

0353 

9357 

9853 

035F 

0363 

0367 

OBOB 

03<>F 

OOVE 


087V 


20283420 
48455829 
2007 


00454554 

4SS22053 

54415254 

494E4720 

41444452 

45635320 

23342048 

45552920 

54484545 

20535041 

43452041 

4E44P045 

45544552 

OOOA 

45454449 
45472041 
44445245 
5353204F 
46204441 
54412054 
4F204245 
20524543 
4F524445 
44250004 
54555245 
204F4E20 
52454 34F 
52442041 
4E442048 
49542052 
45545552 
4E000407 


MS6L6 

1 

1 

MESCTt 


E0U>«>ME$«6 


03 CO«*ENTFM STAOTINC ADDRESS (A HFXl THEN SFACE AND ENTE«* 


03 •ENOING A00R5SS OF DATA TO BE RFCOROEO* • tCO ,LF 


03 'TURN ON RECORD AN9 HIT RETURN* «CR tLF.BELL 


NS7LC EOU 4-M5S67 


Loc: 


OS 01 


INTERRUPT JUMP TABLE 




riTrm>*iT^ 


•a«0 MACRO ASSEMBLER* VER 2* A CRMORS ■ O.PACE 12 


OREO 


QR6 

0BE9N 

OSEO 


▼ABLE! 


easo 

C3970S 

JMP 

BAO 

0SE3 

00 

NOR 


OBEA 

C3720R 

JRP 

COOO 

OBE7 

00 

NOB 


oaEs 

C3970B 

JMP 

BAO 

OBEB 

00 

NQR 


03EC 

C3970S 

JMR 

BAO 

OSEF 

00 

NOR 


OBRd 

C3970S 

jmb 

BAO 

0SR3 

00 

NOR 


OBFA 

C39709 

JMR 

BAO 

03F7 

09 

NOR 


OBRB 

C3970B 

JMP 

BAO 

03FB 

00 

NOR 


OBFC 

C3970S 

JNR 

BAD 

OBRF 

00 

NOR 




END 



NO PROGRAM ERRORS 
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•0*0 MACRO ASSEMBLFR* VCR 2.A 


ERRORS ■ O.RACF 13 


SYMBOL TABLE 


* 

A 

0307 

AGAIN 

0R3A 

ATCO 

0B64 

R 

0000 

BACK 

03B9 

BAD 

0«9T 

BELL 

eooT 

BREAK 

01C2 

C 

0001 

CL 

OACA 

CNVBN 

einp 

CO 

OlEB 

CR 

OOOO 

CRRU9 

ORTI 

CSPUN 

091-» 

O 

OC02 

OMA 

0 009 

E 

0003 

rN»EP 

0912 

GET 

06F3 

GETCH 

0220 

GFTMX 

0??T 

GFTNM 

029C 

GOOD 

0672 

M 

0004 

HILO 

02AI 

ICMt 

OOFP 

ICW? 

OOOB 

INIT 

0930 

INT 

0900 • 

L 

0009 

LOOP A 

OBSB 

LR 

OOOA 

LOAPN 

0009 

LOC^ 

OBTl 

LSO 

OOOA 

M 

000& 

MESGI 

0 99B 

PPSG2 

0A93 

MESC3 

0A6F 

MESC4 

0A92 

MH SGS 

OAAO 

MVSG6 

OABF 

MESG7 

OAFS 

MO'a 

0014 

MSILC 

OOC9 

MS2LC 

OOlC 

MS3L C 

0023 

msalc 

OOtB 

MS5LG 

0012 

M?*LG 

0036 

MSTLC 

007F. 

MSB 

OOOB 

NMOUT 

0*C9 

nc «2 

002 0 

OOTMS 

0961 

PLACE 

09AO 

PSW 

0006 

PUB 

OO’F 

SAPP 

0664 

SSL 

090 A 4 

SKIP 

OOOF 

SP 

006* 

START 

06E6 

STMLF 

0355 

TABLE 

09EO * 

TERR 

3FP9 

TRANS 

0B9I 

TX 

05A2 

WAIT 

09BE 
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“iff ■ 


APPENDIX H 


CONTROLLER ROM CONTENTS 




< 

^ Q _l _l 

q; o 3: a: 

HEX 

, Equivalent 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 4 1 0 0 0 

8 6 1 0 0 0 

8 4 1 0 0 0 

8 5 1 0 0 0 

8 4 1 0 0 0 

9 4 1 0 0 1 

9 4 1 0 0 1 

8 4 1 0 0 0 

8 4 1 0 0 0 


gii: ;:v.Q t-ooo 
oooo_i_j oH-H-f- pi-PS 


0 1 0 0 0 1 i 0 1 o"n f 
0100 olio 1011 
0100 0110 1011 
OlOO 0110 1011 

0100 0110 1011 
0100 Olio 1011 
0100 0 1 10 1011 
0100 0 110 1011 
0100 0010 0111 
0100 0010 0111 
0100 0010 0111 
0100 0010 0111 
0100 0011 0111 
0100 0011 0111 
0110 0011 0111 
0100 0011 0111 
0 10 1 0 0 11 1111 
0100 0011 1111 
0100 0011 1111 
0100 0011 1111 
0100 0011 1111 
0100 0011 1111 







I ' «i^- 


13 

"V. a 

of o 


Equivalent 

8 

4 

8 

4 

A 

4 

A 

4 

8 

4 

8 

4 

C 

4 

C 

4 

4 

4 

4 

4 

C 

4 

C 

4 

8 

4 

8 

4 

8 

4 

8 

4 

8 

4 

8 

4 

8 

4 

8 

4 

8 

4 


C <C 


0 0 
0 0 
1 0 
1 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


oc a: oi 
<13: Q 

H- ^ Q Q 
CO CO _J _l 


0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 
0 10 0 


_j 

(_} 

H* 

_J 

O 

_J 

O 

_l 

o 

.J 

u 


Q£ 

z 

liJ 

M 

f— 

«\J 

CO 

<• 

t— • 

CO 


O 

t- 

»- 

1- 



1- 

a 


0 0 11 
0 0 11 
0 0 11 
0 0 11 
0 0 11 
0 0 11 


10 11 
10 11 
10 11 
10 11 
0 0 0 0 
0 0 0 0 
0 0 0 0 


1110 

1110 

1110 

1110 

1110 

1110 


0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 


1111 

1111 

1111 

1111 

1111 

1111 

1111 


HEX 

Equivalent 
3 f E 1 




' x,>< 




i 

i 

i 


h 

[■ 

r 

\ 

j 

I 

t 

I 


fc 


ae 



1 3 

Q 

oi o 

HEX 

Equivalent 

WLAT 

RLAT 

START 

1 3= 

to 

LD/R 

LDDR 

43 

8 

4 

1 0 

0 

0 

0 

"T 

0 

0 

44 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

45 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

46 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

47 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

48 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

49 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

50 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

51 

8 

0 

1 0 

0 

0 

0 

0 

0 

0 

52 

8 

0 

1 0 

0 

0 

0 

0 

0 

0 

53 

8 

0 

1 0 

0 

0 

0 

0 

0 

0 

54 

8 

0 

1 0 

0 

0 

0 

0 

0 

0 

55 

8 

0. 

1 0 

0 

0 

0 

0 

0 

0 

56 

8 

0 

1 0 

0 

0 

0 

0 

0 

0 

57 

8 

0 

1 0 

0 

0 

0 

0 

0 

0 

58 

8 

0 

1 0 

0 

0 

0 

0 

0 

0 

59 

8 

C 

1 0 

0 

0 

1 

1 

0 

0 

60 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

61 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

62 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 

63 

8 

4 

1 0 

0 

0 

0 

1 

0 

0 


27 






CZT CLK 

TICK 

T2CL 

T3CL 

T4CL 

TIK 

T3R 

DIEN 

HEX 

Equivalent 

0 

0 

0 

0 

1 

1 

1 

1 

0 

F 

0 

1 

0 

0 

1 

1 

1 

1 

4 

F 

0 

1 

0 

0 

1 

1 

1 

1 

4 

F 

0 

1 

0 

0 

1 

1 

1 

1 

4 

F 

0 

1 

0 

0 

1 

1 

1 

1 

4 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 ' 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 

0 

1 

1 

0 

1 

1 

1 

1 

6 

F 


6 F 




The REMOTE SENSING CENTER was established by authority of the Board of Directors of 
the Texas A&M University System on February 27, 1968. The CENTER is a consortium of four 
colleges of the University; Agriculture, Engineering, Geosciences, and Science. This unique 
organization concentrates on the development and utilization of remote sensing techniques and 
technology for a broad range of applications to the betterment of mankind. 










